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V  > 

'  The  Spot  of  Arago  is  the  bright  diffraction  spot  in  the  center  of 
the  shadow  behind  a  circular  obstacle.  Because  this  pattern  is  the 
result  of  diffraction  at  the  edge  of  the  obstacle,  the  spot  contains 
information  about  the  incident  beam.  If  the  incident  beam  contains 
aberrations,  the  shape  and  intensity  of  the  diffraction  pattern  should 
change.  1 

'For  points  on  and  near  the  optical  axis,  analytical  solutions  to 
the  Rayleigh-Sommerfeld  diffraction  integral  are  possible  for  a  circu¬ 
lar  obscuration  and  an  annular  aperture  illuminated  by  a  plane  wave  of 
uniform  or  gaussian  intensity  distribution  and  containing  defocus 
and/or  spherical  aberration.  Computer  studies  of  these  cases  show 
excellent  agreement  with  experiment.  The  amount  of  defocus  and  spher¬ 
ical  aberration  can  be  deduced  by  intensity  measurements  and  shift  in 
positions  of  on  axis  intensity  extrema  behind  an  annular  aperture.  An 
empirical  study  of  astigmatism  revealed  a  predictable  change  in  the 
diffraction  pattern  allowing  for  the  verification  of  the  aberration  to 
within  .05  wavelengths. 

The  Spot  of  Arago  is  inherent  in  many  optical  systems  due  to  the 


geometry  of  components.  The  technique  of  aberration  determination  by 
changes  in  the  diffraction  pattern  should  be  useful  in  such  systems. 
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aperture,  e  *  .625*  Viewing  distance:  2  -  2.05m. 
Magnification  6X.  Exposure  time  1  second. 

10°  countertilt.  Minimum  on  axis . 

6.70  5/16"  diameter  obstacle  with  1/2"  diameter 

aperture,  e  =  .625.  Viewing  distance:  2  -  2.05m. 
Magnification  6X.  Exposure  time  1  second. 

11°  countertilt.  Maximum  on  axis . 

6.71  5/16"  diameter  obstacle  with  1/2"  diameter 

aperture,  c  =  .625.  Viewing  distance:  2  -  2.05m. 
Magnification  6X.  Exposure  time  1  second. 

12°  countertilt.  Maximum  on  axis . 

6.72  5/16"  diameter  obstacle  with  1/2"  diameter 

aperture,  e  =  .625.  Viewing  distance:  2-2. 

05m.  Magnification  6X.  Exposure  time  1  second. 

13°  countertilt.  Minimum  on  axis . 


6.73  5/16"  diameter  obstacle  with  1/2"  diameter 

aperture,  c  =  .625-  Viewing  distance:  2  -  2.05m. 
Magnification  6X.  Exposure  time  1  second. 

13°  countertilt.  Maximum  on  axis . 
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6.74 

6.75 

6.76 

6.77 

6.78 

6.79 

6. 80 

6.81 

6.82 

6.85 
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5/16"  diameter  obstacle  with  1/2"  diameter 

aperture.  £  =  .625*  Viewing  distance:  2-2. 
05m.  Magnification  6X.  Exposure  time  1  second. 


14°  countertilt.  Minimum  on  axis . 228 

5/16"  diameter  oostacle  with  1/2"  diameter 

aperture.  £  =  .625*  Viewing  distance:  2  -  2.05m. 
Magnification  6X.  Exposure  time  1  second. 

14°  countertilt.  Maximum  on  axis . 225 

5/16"  diameter  obstacle  with  1/2"  diameter 

aperture.  £  -  .625*  Viewing  distance:  2  -  2.05m. 

Magnification  6X.  Exposure  time  1  second. 

15°  countertilt.  Minimum  on  axis . 229 

5/16"  diameter  oostacle  with  1/2"  diameter 

aperture.  £  =  .625-  Viewing  distance:  2  -  2.05m. 

Magnification  6X.  Exposure  time  1  second. 

15°  countertilt.  Maximum  on  axis . 229 

3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  e  =  .68.  Viewing  distance:  .6  -  ,65m. 

Magnification  4X.  Exposure  time  1  second. 

22°  countertilt.  Minimum  on  axis . 230 

3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  c  =  .68.  Viewing  distance,:  .6  -  .65m. 

Magnification  4X.  Exposure  time  1  second. 

22°  countertilt.  Maximum  on  axis . 230 

3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  c  =  .68.  Viewing  distance:  .6  -  .65m. 
Magnification  4X.  Exposure  time  1  second. 

24°  countertilt.  Minimum  on  axis . 231 

3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  e  =  .68.  Viewing  distance:  .6  -  ,65m. 
Magnification  4X.  Exposure  time  1  second. 

24°  countertilt.  Maximum  on  axis . 231 

3/16"  diameter  obstacle  with  7mm  oiameter 

aperture,  e  =  .63.  Viewing  distance:  .6  -  .65m. 
Magnification  4X.  Exposure  time  1  second. 

26°  countertilt.  Minimum  on  axis . 232 

3/1 6"  diameter  obstacle  with  7mm  diameter 

aperture,  c  =  .68.  Viewing  distance:  .6  -  .65m. 
Magnification  4X.  Exposure  time  1  second. 

26°  countertilt.  Maximum  on  axis . 232 
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6.84  3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  e  =  .68.  Viewing  distance:  .6  -  .65m. 
Magnification  4X.  Exposure  time  1  second. 

28°  countertilt.  Maximum  on  axis . 253 

6.85  3/16"  diameter  obstacle  with  7mm  diameter 

aperture,  e  =  .68.  Viewing  distance:  .6  -  .65m. 
Magnification  4X.  Exposure  time  1  second. 

30°  countertilt.  Maximum  on  axis . 233 

6.86  1/2"  diameter  obstacle  with  5/8"  diameter 

aperture.  e  =  .8.  Viewing  distance:  2.05  -  2.15m. 
Magnification  6X.  Exposure  time  1  second. 

9°  countertilt.  Minimum  on  axis . 234 

6.87  1/2"  diameter  obstacle  with  5/8"  diameter 

aperture,  e  =  .8.  Viewing  distance:  2.05  -  2.15m. 
Magnification  6X.  Exposure  time  1  second. 

10°  countertilt.  Minimum  on  axis . 234 

6.88  1/2"  diameter  obstacle  with  5/8"  diameter 

aperture.  r.  =  .8.  Viewing  distance:  2.05  -  2.15m. 
Magnification  6X.  Exposure  time  1  second. 

10°  countertilt.  Maximum  on  axis . 235 

6.89  1/2"  diameter  obstacle  with  5/8"  diameter 

aperture,  c  =  .8.  Viewing  distance:  2.05  -  2.15m. 
Magnification  6X.  Exposure  time  1  second. 

12°  countertilt.  Maximum  on  axis . 235 

6.90  3/16"  diameter  obstacle  with  5/8"  diameter 

aperture.  Viewing  distance  2m.  8°  counter¬ 
tilt . 241 

6.91  3/16"  diameter  obstacle  with  5/8"  diameter 

aperture.  Viewing  distance  2.05m.  8°  counter¬ 
tilt . 242 
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Chapter  1 


INTRODUCTION 


Annular  apertures  are  commonplace  in  optical  systems.  For  exam¬ 
ple,  Cassegrain  telescopes  (reflecting)  have  a  secondary  mirror  which 
acts  as  a  central  obscuration  in  an  otherwise  circular  field.  A 
circular  stop  deliberately  placed  on  the  objective  of  a  refracting 


Figure  1.1 

Cassegrain  Mirror  System  and  Cross-section  of  Field^ 


telescope  can  increase  the  depth  of  focus  and  the  resolving  power  of 

2 

the  system. 
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Statement  of  the  Problem 

Laser  systems  often  have  annular  apertures  as  part  of  an  optical 
resonator  configuration.  The  diagram  below  depicts  an  elliptical  out- 
coupling  mirror  set-up  that  acts  as  the  optical  analog  of  an  annular 
aperture. 


incident  beam 

Figure  1.2 

Elliptical  Scraper  Mirror  and  Its  Optical  Analog 


This  last  configuration  appears  in  the  design  of  many  High  Energy 
Laser  (HEL)  resonators.  It  has  created  a  particular  nuisance  for  its 
developers.  In  the  center  of  the  shadow  produced  by  this  effective 
annular  aperture  is  a  hot  spot.  (The  frequencies  of  the  HEL  output 
are  in  the  infrared  range.)  The  spot  is  a  consequence  of  the  con¬ 
structive  interference  of  radiation  diffracted  from  the  outer  and 
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inner  edges  of  the  "aperture."  This  spot  is  related  to  the  spot  of 
Arago  —  the  diffraction  spot  formed  by  light  incident  upon  a  circular 
obstacle.-^  The  intensity  of  the  spot  can  be  as  high  as  four  times  the 
intensity  of  the  incident  beam.  This  makes  it  an  unwanted  addition  to 
a  designated  "shadow"  area.  Left  unchecked  the  spot  can  burn  up  vital 
optical  components  or  produce  thermal  deformations  in  components. 
Designs  must  include  provisions  to  diffract  or  reflect  the  aggravating 
manifestation  into  an  optical  dump.  The  attitude  of  the  HEL  workers 
is  reminiscent  of  Lady  Macbeth's  plea  "Out,  damned  spot!  out,  I  say!"^ 

It  is  frequently  necessary  to  perform  real-time  HEL  diagnostics 
to  test  beam  quality.  This  requires  additional  gratings,  beam  split¬ 
ters,  or  other  beam-sampling  components  which  add  to  the  overall 
weight,  size  and  complexity  of  the  design.  If  the  hot  spot  contains 
information  about  beam  quality,  namely,  if  its  shape  and  size  is 
characteristic  of  aberrations  in  the  incident  wave  front,  it  could  be 
used  as  an  integral  part  of  the  required  diagnostic  process.  Such  a 
scheme  would  eliminate  the  need  for  a  lot  of  hardware. 

James  E.  Harvey  and  James  L.  Forgham  have  demonstrated  that  the 
shape  of  the  spot  in  the  shadow  of  a  circular  obstacle  illuminated  by 
a  plane  wave  does  change  with  increased  aberration  (astigmatism).'1 

The  feasibility  of  a  diagnostic  package  using  the  spot  as  input 
depends  on  the  answers  to  the  following  questions: 

1.  What  types  of  aberrations  produce  measurable  changes  in  the 
spot? 

2.  How  sensitive  are  changes  in  the  spot  to  changes  in  amount  of 
aberration? 


3*  Given  specific  spot  characteristics,  can  the  aberrations 
present  be  determined? 

Thesis  Content 

This  thesis  will  address  itself  to  the  solutions  of  the  preceding 
questions  for  the  following  Seidel  aberrations:  defocus,  spherical 
aberration,  coma  and  astigmatism.  First  the  remarkable  history  of  the 
spot  of  Arago  will  be  traced  from  the  eighteenth  century  to  current 
research.  A  detailed  analysis  for  the  unaberrated  cases  of  the  spots 
produced  by  circular  obscurations  and  annular  apertures  will  serve  as 
standards  against  which  to  compare  the  aberrated  cases.  It  turns  out 
that  analytic  solutions  to  the  Rayleigh-Sommerfeld  diffraction  inte¬ 
gral  are  possible  for  the  rotationally  symmetric  aberrations  (defocus, 
spherical  aberration)  for  certain  restrictive,  but  physically  impor¬ 
tant  cases.  The  non- rotational ly  symmetric  aberrations  (coma,  astig¬ 
matism)  require  numerical  or  empirical  analysis.  Where  possible, 
theory  is  tested  against  experimental  results.  Each  aberration  is 
handled  separately  but  various  strategies  to  separate  combinations  of 
aberrations  are  discussed.  The  main  body  of  the  thesis  culminates  in 
a  discussion  of  experimental  and  theoretical  results,  suggestions  for 
future  research,  and  overall  conclusions.  Several  appendices  are 
included  which  contain  material  of  historical  interest,  lengthy  deri¬ 
vations  and  computer  program  listings. 

Potential  applications  exist  for  a  variety  of  optical  systems; 
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Warren  J.  Smith,  Modern  Optical  Engineering  (New  York:  McGraw 
Hill,  1966),  p.  38. 

*  W.  T.  Welford,  "Use  of  Annular  Apertures  to  Increase  Focal 
Depth,"  Journal  of  the  Optical  Society  of  America,  50  (August  I960), 
749;  E.  H.  Linfoot  and  E.  Wolf,  "Diffraction  Images  in  Systems  with  an 
Annular  Aperture,"  Procedings  of  the  Physical  Society  B,  66  (1953), 
145. 


^  Max  Born  and  Emil  Wolf,  Principles  of  Optics,  6th  ed.  (Oxford: 
Pergamon  Press,  1980),  p.  375. 

^  Macbeth  V.i.38. 

^  James  E.  Harvey  and  James  L.  Forgha,.,  "The  Spot  of  Arago:  New 
Relevance  for  an  Old  Phenomenon,"  American  Journal  of  Physics,  (Janu¬ 
ary  1984). 


One  of  the  great  ideological  clashes  in  the  history  of  physics  is 
the  controversy  between  the  corpuscular  and  wave  theories  of  light. 

The  center  of  the  conflict  rested  on  the  very  nature  of  light.  Advo¬ 
cates  of  the  corpuscular  theory  agreed  with  Sir  Isaac  Newton  that 
light  was  composed  of  minute  particles.  The  wave  theorists  sided  with 
Rene  Descartes,  Robert  Hooke  and  Christian  Kuyghens  who  felt  light  was 
the  result  of  waves  undulating  in  a  medium  called  the  ether.1 

By  the  beginning  of  the  nineteenth  century  the  corpuscular  theory 
was  firmly  ensconced  as  the  preferred  theory.  Peter  Anton  Pav  pro¬ 
poses  the  popularity  of  the  particle  theory  originated  in  the  philo¬ 
sophical  battle  between  inductivism  and  deductivism: 


What  occupied  scientists  was  not  a  battle  of  waves  vs.  corpuscles 
.  .  .  nor  blind  allegiance  to  Newton's  authority.  Rather,  it  was 
a  basic  and  deeply  felt  difference  of  opinion  as  to  how  and  where 
scientific  truth  was  to  be  found  and  whether  the  ultimate  canons 

p 

of  scientific  justification  were  to  be  experimental  or  deductive. 


Newton's  approach  in  his  work  Opticks  was  experimental  and  this  had 

great  appeal  for  the  scientists  of  the  eighteenth  century.  The  cor- 

•3 

puscular  theory  was  adopted  almost  by  consequence. 


To  win  over  the  scientific  community  the  wave  theorists  needed  to 
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show  experimentally  the  conclusive  advantages  of  their  theory.  The 
first  advance  was  made  by  Thomas  Young.  He  conducted  a  series  of 
experiments  exploring  a  variety  of  diffraction  effects  (including  the 
famous  double  slit  experiment).  Despite  publication  of  his  papers  and 
the  delivery  of  a  lecture  and  demonstration  before  the  Royal  Society 
of  London  in  1801,  Young's  work  won  little  immediate  support  for  the 
wave  theory.^  However  he  did  gain  a  very  important  convert  —  the 
Frenchman  Dominique  Francois  Jean  Arago. 

In  1811  Arago  studied  the  Newton's  rings  phenomena  in  calcite  in 
conjunction  with  his  work  on  polarization.  The  failure  of  the  corpus¬ 
cular  theory  to  explain  many  of  his  results  encouraged  him  to  consider 
the  wave  theory.  He  began  a  correspondence  with  Young  and  was  soon 
convinced  of  the  truth  of  Young's  theory. ^ 

The  next  advance  in  the  cause  of  the  wave  theory  came  from  Augus¬ 
tin  Jean  Fresnel.  He  met  briefly  with  Arago  in  1815-  Arago  intro¬ 
duced  Fresnel  to  Young's  treatises  and  persuaded  him  to  come  to  Paris 
in  1816  to  begin  his  researches  in  diffraction.  The  timing  of  the 
research  was  ripe,  for  in  January  1817  a  commission  composed  of  Pierre 
Simon  Laplace,  Jean-Baptiste  Biot,  Claude  Louis  Berthollet  and  Jacques 
Alexandre  Cesar  Charles  chose  diffraction  as  the  subject  for  a  prize 
sponsored  by  the  Academie  des  Sciences.0 

Fresnel  entered  his  memoir  on  diffraction  in  the  contest.  In 
July  1818  a  panel  of  five  was  chosen  to  judge  the  entries.  The  panel 
members  were  Laplace,  Biot,  Arago,  Simeon  Denis  Poisson  and  Joseph 
Louis  Gay-Lussac.  Laplace,  Biot  and  Poisson  were  staunch  adherents  of 
the  particle  theory  of  light  while  Arago  advocated  Fresnel's  view. 
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Gay-Lussac,  a  chemist,  was  impartial. 

In  reviewing  Fresnel's  theory,  Poisson  saw  that  the  simple  geome¬ 
try  of  light  diffracted  by  an  opaque  circular  disc  suggested  the 
existence  of  a  bright  spot  of  light  in  the  center  of  the  disc's 
shadow.  This,  he  felt,  was  contrary  to  experience  and  thus  Fresnel's 

O 

theory  must  be  wrong. 

Arago  tested  Poisson's  prediction  in  the  laboratory.  He  affixed 
a  disc  2mm  in  diameter  to  a  plate  of  glass  and  using  sunlight  observed 
the  spot  in  the  center  of  the  shadow.^  (This  is  an  irony  Pav  should 
appreciate  —  Poisson  attempted  to  discredit  Fresnel's  wave  theory 
based  on  deductive  reasoning;  Arago  vindicated  Fresnel  through  experi¬ 
ment.)  Since  this  time  the  phenomena  has  been  known  as  either  Pois¬ 
son's  spot  or  Arago's  spot.  (This  paper  adopts  the  terminology  of 
crediting  Arago,  who  trusted  in  the  spot's  existence.)  Arago's  dis¬ 
covery  was  powerful  evidence  in  favor  of  the  wave  theory  of  light. 
Fresnel  won  the  prize  from  the  Academie  and  his  entry  became  known  as 
his  Crown  Memoir.^ 

Although  the  demonstration  of  the  spot  did  much  to  usurp  the 
corpuscular  theory,  the  controversy  between  waves  and  particles  con¬ 
tinued  through  the  nineteenth  century  in  a  series  of  incarnations. 


Earlier  Observations  of  the  Spot 


Many  sources  mention  that  Arago  was  not  the  first  to  observe 
the  bright  point  of  light.  John  William  Strutt  (Lord  Rayleigh)  and 
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E.  Verdet  in  his  Lecons  d'Optique  Phisique  note  that  Joseph  Nicholas 
Delisle,  a  French  astronomer,  discovered  a  bright  spot  in  the  center 
of  a  circular  shadow  "in  the  earlier  half  of  the  18th  century,"  a 
discovery  which  subsequently  "passed  into  oblivion."1^  John  Strong, 
Max  Born  and  Emil  Wolf  credit  Maraldi,  another  French  astronomer,  with 
the  finding.  Born  and  Wolf  say  that  Maraldi’s  observation  occured  "a 
century  earlier"  than  Arago's  experiment  and  "had  bee’"  forgotten."^ 
Strong  claims  that  Maraldi  found  the  bright  spot  "over  half  a  century 
earlier. 

Joseph  Nicholas  Delisle  (1688-1768)  was  one  of  a  family  of  famous 
French  scientists.  One  brother,  Guillame,  was  a  geographer  and  histor¬ 
ian,  another  brother,  Louis,  was  also  an  astronomer.  Catherine  the 
Great  invited  Joseph  and  Louis  to  come  to  Russia  in  1726  to  run  the 
observatory  at  St.  Petersburg.  Louis  died  while  in  Russia  and  Joseph 
Nicholas  came  back  to  Paris.  Delisle  became  a  member  of  the  Academie 
des  Sciences  in  1714  and  in  1724  was  elected  to  the  Royal  Society  of 
London. ^  ^ 

As  a  member  of  the  Academie  des  Sciences,  Deslisle  could  have  met 
and  worked  with  two  Maraldis.  Neither  Born  and  Wolf's  text  nor 
Strong's  book  mentions  Maraldi's  first  name. 

Jacques-Phi 1 lipe  Maraldi  (1665-1729)  was  the  nephew  of  Dominique 
Cassini  —  the  discoverer  of  the  famous  division  in  the  rings  of 
Saturn.  Jacques-Phil lipe  also  became  an  astronomer  and  joined  the 
Academie  des  Sciences  in  1694. 

Jean-Dominique  Maraldi  (1709-1788)  was  the  nephew  of  Jacques- 
Phil  lip  Maraldi.  He  was  inducted  into  the  Academie  des  Sciences  in 
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The  vague  references  from  Rayleigh,  Verdet,  Strong  and  Born  and 
Wolf  do  not  allow  one  to  pinpoint  the  date  of  the  first  observation, 
the  circumstances  surrounding  the  discovery,  who  made  it  or  even  which 
Maraldi  was  involved.  Because  Delisle  and  both  Maraldis  lived  and 
worked  in  Paris  at  overlapping  times  and  belonged  to  the  same  scien¬ 
tific  organization  it  is  possible  that  the  "discovery"  by  the  second 
party  was  actually  a  verification  of  a  phenomenon  observed  by  the 
first. 

The  sources  also  suggest  the  observations  of  Delisle  and  Maraldi 
about  the  spot  faded  into  the  limbo  of  untimely  scientific  discov¬ 
eries.  That  Poisson  was  unaware  of  these  findings  is  evident  from  the 
challenge  he  made  against  Fresnel’s  theory.  Arago,  however,  did  know 
about  the  earlier  sightings.  On  February  26,  1816  (two  years  before 
he  sat  on  the  judging  commission  with  Poisson)  he  presented  to  the 
National  Institute  observations  about  diffraction  and  an  early  draft 
of  Fresnel's  memoir.  Here  he  mentions  that  besides  Grimaldi,  Maraldi 
and  Delisle  showed  that  light  could  penetrate  the  geometric  shadow 
behind  an  opaque  bodyJ^  This  knowledge  may  have  provided  additional 
impetus  for  Arago  to  verify  Poisson's  skeptical  prediction. 

Rings  Around  the  Spot 

The  ease  with  which  one  could  produce  the  spot  made  it  an  ideal 
subject  for  study  by  the  new  legions  of  wave  theorists.  The  following 


is  a  description  of  an  experiment  by  Richard  Potter  to  test  the  inten¬ 
sity  of  the  spot. 


I  prepared  .  .  .  discs  of  brass  of  1/10,  2/10,  3/10,  4/10  and 
7/10  and  with  watch  makers  implements  of  1/20  inch.  The  discs 
were  attached  by  thin  films  of  cement  to  a  plate  of  glass.  The 
sun's  light  was  reflected  horizontally  through  a  window  shutter 
of  a  darkened  room,  and  the  sun's  image,  formed  by  a  lens  of  1/6 
inch  focal  length  was  at  60  inches  distance  from  the  disc  at  the 
same  time  the  focus  of  the  eye  lens,  by  which  the  phenomena  of 
diffraction  were  examined,  was  at  60  inches  distance  on  the  other 
side.  When  the  whole  was  adjusted,  on  looking  through  an  eyelens 
of  about  1  inch  focal  length,  at  the  centre  of  the  shadow  cast  by 
the  disc,  there  was  seen  a  bright  central  spot,  of  a  white  color 
slightly  tinged  with  brown,  surrounded  by  a  greater  or  lesser 
number  of  coloured  rings,  according  to  the  size  of  the  discJ"^ 


The  colors  observed  by  Potter  stem  from  the  fact  that  light  from  the 
sun  is  not  monochromatic. 

The  concentric  ring  structure  was  given  a  mathematical  foundation 

by  the  British  Astronomer  Royal,  George  Biddell  Airy.  In  a  paper 

entitled  "On  the  Diffraction  of  an  Annular  Aperture,"  Airy  constructed 

a  table  by  which  one  could  deduce  the  spacing  of  the  ring  structure 

produced  in  three  cases:  monochromatic  light  falling  on  a  circular 
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aperture,  a  circular  obscuration  or  an  annular  aperture. 


Lommel's  Work 


The  definitive  study  of  diffraction  by  a  circular  obstacle  and  a 
circular  aperture  was  made  by  Emil  Lommel  in  the  1880s.  Lommel  solved 
the  diffraction  integral  by  using  Bessel  functions.  He  found  the 


intensity  distribution  of  light  in  the  shadow  of  a  circular  obstacle 


satisfies  the  expression: 


IT  =  C2  +  S2  (II-1  ) 

2 

where  M  is  the  intensity  and  C  and  S  are  integrals  whose  solutions 
contain  the  Loramel  functions: 

Born  and  Wolf  list  the  Lommel  functions  as 


Un(u,  v)  =  £  (-1)  W+2S  ,-n+2s(v) 


s=0 


,  \  ,  vS.uNn+2s  ,  x 

V  <n,  v)  -  £  (-»  (7)  J„+2s<v) 

s=0 


(II-2) 


(H-3: 


where  u  and  v  are  dimensionless  parameters  based  on  the  geometry  of 
the  set-up  and  Jn+2S  the  Besse2  function  of  the  first  kind  and  n  + 
2sth  order.1^  The  Lommel  functions  have  the  disadvantage  of  being 
slowly  convergent. 

Lommel  performed  the  experiment  to  verify  his  theory.  He  used 

sunlight  and  a  prism  to  obtain  fairly  monochromatic  light.  With  an 

eyelens  etched  with  a  micrometer  he  measured  the  ring  structure  of 

both  the  obstacle  and  the  aperture  and  found  it  to  be  in  excellent 

on 

agreement  with  theory. 


Experiments  at  the  Turn  of  the  Century 


It  is  interesting  that  Lommel  relied  on  the  somewhat  crude  method 


of  sunlight  and  prism  to  make  a  light  beam  with  narrow  band  width. 


The  nineteenth  century  saw  major  advances  in  the  production  of  high 

intensity  light  sources.  In  1825  Thomas  Drummond  invented  lime- 
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light.  The  carbon  arc-light  was  developed  by  a  number  of  people. 

Frederick  Holmes  built  a  cumbersome  working  model  in  1857.  By  the  mid 

22 

1870s  practical  commercial  arc-lamps  were  being  sold. 

The  widespread  use  of  arc-lamps  in  laboratories  at  the  turn  of 
the  century  renewed  interest  in  the  spot  of  Arago.  Using  a  spec¬ 
trometer  one  could  isolate  an  intense  monochromatic  signal  from  the 
carbon  spectrum.  This  was  important  because  photographs  of  the  phe¬ 
nomenon  were  possible. 

W.  Arkadiew  published  a  famous  series  of  photographs  in  1913- 
Using  a  variety  of  lens  combinations  he  altered  the  virtual  distance 
between  object  and  image  planes  to  observe  the  change  in  spot  size.2^ 
In  April  1914  two  papers  appeared  which  demonstrated  that  a 
circular  obstacle  could  be  used  in  place  of  a  lens  to  produce  an 
image.  Alfred  W.  Porter  noticed  that  "the  bright  spot  is  in  reality 
an  image  of  the  source  of  light. "2^  Porter  used  a  set-up  similar  to 
the  one  described  by  Mason  E.  Hufford  and  depicted  below. 


box 


arc  collimating  opening  suspended  photographic 

lamp  lens  with  sphere  plate 

transparency 

Figure  2.1 

Experimental  Set-up  to  Demonstrate  Imaging  by  a  Disc 
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Hufford  reported  that  "each  point  in  the  aperture  (opening  to  the  box) 
produces  a  corresponding  effect  in  the  shadow.  This  suggests,  there¬ 
fore,  that  the  aperture  may  be  any  sort  of  figure  and  in  the  shadow  of 
the  disc  or  sphere  there  should  be  an  inverted  image  of  it."^  Porter 
used  a  triangular  opening  to  produce  a  three-sided  image.  Hufford, 
with  a  bit  more  flair,  made  a  monogram  out  of  the  letters  I  ana  U  (he 
was  a  professor  at  Indiana  University)  as  one  opening  and  a  transpar¬ 
ency  of  the  profile  of  Woodrow  Wilson  as  another.  The  resultant  image 
of  the  latter  appears  in  Jenkins  and  White. 

Hufford  also  made  photographs  of  the  ring  structure  in  the  shadow 
of  a  circular  obstacle.  Lommel  had  only  measured  the  radii  of  a  few 
rings  immediately  surrounding  the  central  spot.  By  exposing  a  photo¬ 
graphic  plate  for  several  hours  Hufford  got  photographs  showing  a 
large  number  of  rings.  These  rings  could  be  measured  to  check  Lom- 
mel's  theory.  There  was  excellent  agreement  with  theory  (less  than  2 % 
difference  in  radius  size  for  the  outer  rings). ^ 

At  a  meeting  of  the  Physical  Society  of  London  in  1925,  A.  0. 

Rankine  projected  several  images  using  a  7mm  ball  bearing.  He  noted 

that  each  point  in  the  light  source  produces  an  image  in  the  shadow. 

Since  only  one  source  point  is  on  the  optical  axis,  a  circular  obsta- 

28 

cle  creates  some  aberration.  The  aberrations  can  be  understood  by 
noticing  that  points  in  the  source  off  the  optical  axis  do  not  enjoy 
the  special  symmetry  that  the  axial  point  does. 


source 


obstacle 


source 


obstacle 


Figure  2.2 

Disk  Illuminated  by 
On-axis  Source  Point: 
Equal  Optical  Path  Distances 


Figure  2.3 

Disk  Illuminated  by 
Off-axis  Source  Point: 
Different  Optical  Path  Distances 


Optical  path  differences  in  Figure  2.3  lead  to  phase  differences  upon 
diffraction.  If  one  substitutes  a  spherical  obstacle,  the  circular 
cross  section  restores  the  symmetry  for  every  source  point. 


source  obstacle 

Figure  2.4 

Sphere  Illuminated  by  Off-axis  Source  Points: 
Equal  Optical  Path  Distances  for  Each  Point 


Dissimilarities  in  the  intensities  of  the  Arago  spots  for  circular 
obstacles  and  spherical  obstacles  was  the  subject  of  a  1926  paper  by 
C.  V.  Raman  and  K.  S.  Krishnan.  They  mounted  a  sphere  and  a  disc  to  a 
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glass  plate  and  measured  the  central  spot  intensities  for  both  cases. 
If  the  observation  distance  behind  the  obstacles  was  small,  the  circu¬ 
lar  obstacle  produced  a  brighter  central  spot  than  the  spherical 
obstacle.  Raman  and  Krishnan  explained  it  this  way.  Rays  diffracted 
by  the  edge  of  a  circular  obstacle  can  reach  the  observation  point 
directly.  For  a  spherical  obstacle  one  must  construct  a  cone  tangent 

x  x 

p 

Circular  obstacle  Spherical  obstacle 

Figure  2.5 

Diffraction  Edge  Sources  for  Rays  Converging  at  Point  P 

to  the  back  side  of  the  sphere  with  the  observation  point  as  the  apex. 
"The  disturbance  incident  on  the  surface  of  the  sphere  has  to  creep 
round  it,  as  it  were,  over  the  arc  XY  before  the  rays  diffracted  out 
by  the  sphere  can  reach  the  point  of  observation,  and  must  suffer  a 

very  considerable  dimunation  in  the  process."  As  the  distance  from 

the  obstacle  to  P  increases,  the  arc  XY  decreases  in  length  and  the 
difference  in  intensity  between  the  spherical  and  circular  obstacle 
diffraction  spots  becomes  smaller.^ 

The  promise  of  imaging  by  a  disc  or  sphere  was  never  fully  real¬ 
ized.  The  ring  structure  associated  with  each  image  point  led  to 
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smearing  and  poor  image  quality.  With  no  forseeable  practical  appli¬ 
cation,  interest  in  the  spot  of  Arago  waned  until  the  next  major 
advance  in  production  of  monochromatic  light. 

The  Advent  of  the  Laser 

A  laser  is  a  source  of  coherent,  intense,  monochromatic  light. 

It  provides  a  simple  means  of  observing  the  spot  of  Arago.  A  circular 
or  spherical  obstacle  placed  in  an  expanded  beam  from  a  laser  creates 
a  bright  spot  and  ring  structure  that  in  many  cases  is  visible  with 
the  naked  eye.-^ 

The  importance  of  the  laser  prompted  innovative  uses  for  the 
spot.  One  area  with  great  potential  is  the  use  of  the  spot  in  the 
detection  of  aberrations. 

In  an  interdepartmental  correspondence  from  Hughes  Aircraft, 

J.  A.  Jenny  showed  a  way  to  detect  tilt  in  a  wavefront.  A  plane  wave 

propagating  in  the  direction  of  an  optical  axis  normal  to  a  circular 

obstacle  will  image  a  bright  spot  at  every  point  on  axis  behind  the 

obstacle.  If  the  incident  wave  direction  vector  makes  a  small  angle 

with  the  optical  axis  the  spot  is  offset  from  the  axis  by  a  calculable 

distance.  Thus  a  circular  obstacle  placed  in  a  wavefront  to  be  tested 

■z  1 

can  be  used  to  measure  the  amount  of  tilt  present. 

Harvey  and  Forgham  solved  the  Fresnel  region  Rayleigh-Sommerfeld 
diffraction  integral  for  a  spherical  obstacle  with  an  incident  plane 
wave  of  uniform  or  gaussian  distribution  of  intensity.  The  intensity 


pattern  in  the  image  plane  has  a  JQ  dependence,  where  is  the 
zeroth  order  Bessel  function.  They  also  derived  a  formula  for  the 
intensity  on  axis  behind  an  annular  aperture.  Introducing  known 
amounts  of  aberration  (astigmatism)  in  the  incident  wavefront  they 
demonstrated  the  shape  and  pattern  of  the  spot  and  ring  structure 
changes. 

This  last  experiment  is  the  embarkation  point  for  this  thesis 
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Chapter  5 


THEORY  I:  UNABERRATED  CASES 


Consider  a  monochromatic  plane  wave  incident  upon  an  aperture  in 
an  otherwise  opaque  plane  screen.  The  disturbance  at  a  point  P>  be¬ 
hind  the  screen  is  given  by  the  Rayleigh-Sommerfeld  diffraction  inte¬ 
gral:  ^ 

U(P- )  =  7X  f  f  U(Pj)  —  cos  (n,  1)  r1dr1d01  (III-l) 

plane 


Figure  5.1 

Diffraction  by  an  Aperture  in  a  Plane  Screen 
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X  is  the  wavelength  of  the  incident  wave.  U(Pj)  is  the  disturbance 
function  in  the  plane  of  screen  where: 


=  A(rx,  0j )  eikW(rl’  T(rj  ,  0X) 


(HI-2) 


A(rlt  0j)  =  amplitude  function 
W(r,,  0, )  =  aberration  function 
T(r,,  0 , )  =  aperture  function 


=  1  if  ( r-.  ,  0X)  is  in  the  aperture 


=  0  otherwise 


n  is  a  unit  vector  normal  to  the  plane  of  the  screen 
cos  (n,  l)  is  the  cosine  obliquity  factor 
1  is  the  vector  from  point  P,  to  point  P? 

It  is  possible  to  obtain  'll  in  terms  of  the  factors  r  ,  r  .  0,0 

~  12  12 


and  z,  the  normal  distance  from  point  P2  to  the  screen  plane. 


screen  plane 


observation  plane 


Figure  3.2 

General  Geometrical  Layout 


As  shown  in  the  above  diagram: 


1  =  z  +  r 


(III— 3 ) 


25 


But  by  the  law  of  cosines: 

r32  =  Tj2  +  t2'  -  2r1r2cos  (02  -  0j)  (III-4) 

.*.  I2  =  z2  +  r  2  +  r  2  -  2r  r  cos  (0  -  (5  )  (III-5) 

1  2  12  2  1 

Case  1 :  Circular  Obstacle  with  Uniform  Illumination 

The  first  case  to  consider  is  that  of  a  circular  obstacle  in  the 
screen  plane.  The  diameter  of  the  obstacle  is  d.  The  optical  axis 
will  be  defined  as  the  line  normal  to  the  obstacle  and  passing  through 
the  center  of  the  obstacle,  which  is  the  origin  of  a  cylindrical 
coordinate  system. 


If  the  incident  plane  wave  is  of  uniform  illumination  (constant 
amplitude  A)  and  contains  no  aberrations  then  the  disturbance  function 
U  in  the  obstacle  plane  is  given  by: 

U ( r  ,  0  )  =  AT  ( r  ,  0  ) 

11  11 


(II 1—6 ) 


The  aperture  function  T  can  be  transferred  to  the  limits  of  integra' 
tion  and  the  Rayleigh-Sommerfeld  diffraction  integral  becomes: 


Intensity  on  Axis 

The  cosine  obliquity  factor  can  be  replaced  by  the  term 
is  convenient  to  change  the  variable  of  integration  to  1. 


Figure  3.4 

Determination  of  1  for  On-axis  Points 


dr 


1 


Thus  the  integral  becomes: 


U(0,  0,  z) 


ikl 
e _ 

1 


dl 


(HI-7) 


=-  .  It 

1 


(III-8) 


(II 1-9) 


To  evaluate  this  integral,  integrate  by  parts: 


f(l) 

f’(l) 


1 

1 


g’(l)  =  e 


ikl 


...  1  ikl 

8d)  "  ^  * 


U(0,  0,  2) 


CD 

+ 

^ z2+d2/4 


1 

i 


f  - 


.  ikl 
1  e _ 

kl  1 


(III-10) 


This  second  integral  contains  a  —  term,  otherwise  it  is  identical  to 
the  original  integral.  Unless  one  looks  immediately  behind  the  obsta¬ 
cle  (a  distance  on  the  order  of  a  few  wavelengths  or  less)  then  kl  is 
much  greater  than  1.  Thus  the  second  integral  is  negligibly  small. 
The  expression  for  the  disturbance  at  P2  is: 


U(0,  0,  z) 


2  Az  ik  \/ z2+d2/4 

z^+d2/ 4 


but  since  k 


2H 

X 


U(0,  0,  z)  =  A 


ik  V  z  2  +d 2  /  4 


V 


z-+d^/4 


(III-ll) 


¥hat  does  this  expression  mean  physically?  The  factor  A  is  the 
amplitude  for  the  incident  plane  wave.  z  / +  d  ~/4  is  the  cosine 
obliquity  factor  for  rays  diffracted  by  the  edge  of  the  obstacle  to 


P  .  e^  ^  z  +d  /4  is  phase  difference  from  the  edge  points  to  the 
point  P  .  This  indicates  only  the  rays  diffracted  from  the  edge  of 
the  obstacle  cause  a  disturbance  on  axis.  Since  the  optical  path 
length  is  identical  for  each  edge-diffracted  ray  converging  at  P2  the 
overall  interference  effect  will  be  constructive. 

The  on-axis  intensity  distribution  is  given  by  the  expression: 

1(0,  0,  z)  =  (U(0,  0,  z)|2 

-  A?  z2 

‘  A  ~zr+d'^b  (111-12) 

Substituting  the  incident  intensity  I0  =  Az  in  the  expression  yields: 


1(0,  0,  z)  = 


(111-13) 


Thus  the  intensity  on  axis  behind  a  circular  obscuration  is  equal  to 
the  incident  intensity  times  the  square  of  the  cosine  obliquity  factor 
for  edge  diffracted  rays.  For  large  values  of  z,  d2/4z2  is  much  less 
than  1  (the  cosine  obliquity  factor  in  essence  becomes  1)  and  the 
intensity  on  axis  is  the  same  as  if  no  obstacle  were  present. 

Below  is  a  graph  of  l(0,  0,  z)  /  IQ  as  a  function  of  z  for  a  1cm 


diameter  obstacle. 


■BMC 
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(III— 16) 


These  two  restrictions  insure  that  the  cosine  obliquity  factor  is 
approximately  equal  to  1 . 

With  these  restrictions  it  is  possible  to  expand  1.  By  III-5 


l2=z2+r2+r2-2rr  cos  (0  -  0 

1  2  12  2  r 


For  the  region  where  ~  <<  1  then 
r  2  r  2  r  r 

1  =  2  *  2T  *  IT  ‘  "P  C°S 


(111-17) 


If  1  is  multiplied  by  a  large  number  (like  k)  then  all  the  above 
expansion  terms  must  be  kept,  otherwise  1  is  approximately  equal  to  z. 
The  Rayleigh-Sommerfeld  integral  of  interest  is: 

.  C  00  f  2tt  ikl 


f  00  /"  2tt  ikl 

r2’  02 ’  z)  =  lx  J  ,  J  cos  (£’  ridrid0i 


'd/2  ^0 


(III-7) 


Inserting  the  cosine  obliquity  factor  of  1  and  the  approximate  expan¬ 


sions  of  1  yields: 


2  -  -  -  2tt 


U(r2,  02,  z) 


r,  C  [  r 

A  ikz  ik  — I  I  ik  — 

‘  iXz  6  *  2Z  J  J  62 


d/2  0 


r  r 

-ik  - —  cos  (0  -  0  )  r  dr  d0 

x  e  z  2  llil 


(III-18) 


To  solve  the  integral  over  0,  one  can  use  the  identity  for  the  Bessel 
•5 

function: ^ 


(Ill— 19) 


J 

n 


(x) 


ixcosa  ina , 
e  e  da 


When  n  =  0  and  a  *  (02  -  0j)  then 


4-fi) 

< 


rlr2 


ik 
e  z 


cos  (0,  -  0, ) 


d0 


(III-20) 


Substituting  the  above  in  III-18  gives: 


U(r2,  02,  z) 


2ttA  ikz 
iXz  6 


ik 


e 


2 


r2 

2z 


x 


rle 


ikr. 


(III-21) 


(Up  to  this  point  the  treatment  is  identical  to  Lommel’s  formulation. 
He  proceeded  to  solve  the  integral  by  substituting  in  a  Bessel  func¬ 
tion  identity  and  integrating  by  parts  to  obtain  a  series  of  cascading 
higher  order  Bessel  functions.) 

The  integral  to  be  evaluated  is: 


Both  the  exponential  term  and  the  zero  order  Bessel  function  are 
oscillatory  of  similar  period.  (The  period  of  the  Bessel  function 


tines  the 


approaches  2m.)  The  argument  of  the  Bessel  function  is 


ro 


argument  of  the  exponential.  When  —  <<  1  the  Bessel  function  will 


1 


oscillate  much  mere  slowly  than  the  exponential  term.  One  can  divide 
the  interval  [d/T ,  «■]  into  subintervals  sufficiently  small  that  the 
Bessel  function  term  is  essentially  constant  over  the  subinterval. 

One  can  also  imp:se  the  restriction  that  for  neighboring  subintervals 
the  Bessel  function  terms  are  approximately  equal.  That  is,  if 


m 


represents  the  Bessel  term  for  the  i^  subinterval  then: 


/krir2\  /kri+1r2\ 

>v— )  =  — H 


(111-22) 


Since  this  term  is  taken  as  constant  it  can  be  pulled  out  of  the 
integral.  The  integral  over  the  i**1  subinterval  is  thus: 


x  v. 

It  is  clear  by  III-22  that  the  upper  integration  limit  for  the  itn 
subinterval  will  negate  the  lower  integration  limit  for  the  i  =  1 


subinterval.  This  pairwise  cancellation  will  continue  over  the 
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subintervals  and  leave  two  unpaired  terms  corresponding  to  r  =  d/2 
and  r  =  00 .  Thus  the  original  integral  becomes: 


d/2  (111-24) 

The  final,  radially  symmetric  form  for  the  disturbance  at  P0 
is: 


U(r2,  z) 


2-tA  ikz  ikr„2/2z  z  „  2^  ik 

*  ui e  e  2  ik  ]«\~2r)  e  8z 

,  ikz  ikr^  llz  ikd2/8z  ,  (^X2\ 

=  -A  e  e  2  e 


(111-25) 


(Notice  the  disturbance  function  is  independent  of  angle 

The  intensity  at  P  is: 

2 


I(r„»  z) 

C. 


(III-26) 


Since  J^(0)  =1,  it  is  readily  apparent  the  intensity  on  axis 
equals  the  incident  intensity  (A2),  which  is  expected  for  the  region 

where  —  <<  1  • 

2z 

An  intensity  contour  plot  and  related  radial  intensity  profile 
for  a  3/16  inch  diameter  obstacle  viewed  50cm  away  is  given  below. 
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As  is  characteristic  of  an  Airy  pattern,  the  spot  of  Arago  is 
surrounded  by  a  series  of  increasingly  dimmer  rings.  (The  reference 
to  the  Airy  pattern  is  intentional.  Both  the  pattern  and  the  spot  are 
images  of  a  light  source  produced  by  a  circular  aperture  and  obstacle 
respectively.)  The  ratios  of  the  maximum  intensities  to  the  intensity 
on  axis  as  well  as  the  positions  of  the  maxima  and  minima  can  be 
determined  readily  from  a  table  of  zeroes  for  and  J^. 

Extrema  occur  when: 


-A2  M  j 
2z  0 


hr) M  - 0 


(III-27) 


/kdr,\  /Kdr  \ 

Minima  occur  when  JQ  I  1  =  0  and  maxima  when  l  j  =  0.  The 

table  below  shows  the  radial  distances  (in  terms  of  ~)  and  normalized 

Kd 

intensities  (normalized  to  l(0,  0,  Z))  for  the  first  fourteen  intensi¬ 
ty  extrema  based  on  the  table  of  values  for  JQ  and  found  in  the  CRC 
Standard  Math  Tables,  13th  student  edition,  pages  326-329* 


Table  3> 1 


Extrema  of  Diffraction  Ring  Structure 


Radial  distance 

kd 

X  27 

I (r  /z)  /  1(0, z)  = 

v(“?) 

Extremum  type 

0 

1.0 

Maximum 

2.4C5 

0.0 

Minimum 

3.832 

0.162 

Maximum 

5.502 

0.0 

Minimum 

7.016 

0.09 

Maximum 

8.654 

0.0 

Minimum 

10.174 

0.062 

Maximum 

1 1.792 

0.0 

Minimum 

13.324 

0.048 

Maximum 

14.931 

0.0 

Minimum 

16.471 

0.039 

Maximum 

18.071 

0.0 

Minimum 

19.616 

0.032 

Maximum 

21.212 

0.0 

Minimum 

An  examination  of  the  factors  k,  r  ,  d  and  z  in  the  argument  of 
the  Bessel  function  provides  insight  into  characteristics  of  the  spot 
and  its  ring  structure.  The  argument  is: 

kdr  irdr 

_ 2  _ 2 

2z  Xz  (III-28) 

If  the  argument  is  set  equal  to  a  constant,  then  r,  is  proportional  to 
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X  and  z  and  inversely  proportional  to  d.  The  larger  r is  for  a  given 
argument  the  more  "spread  out"  the  ring  pattern  is  in  the  observation 
plane.  The  pattern  will  increase  in  size  the  longer  the  wavelength  of 
incident  light,  the  farther  along  the  optical  axis  the  observation  is 
made  or  the  smaller  the  diameter  of  the  obstacle.  Likewise  the  pat¬ 
tern  shrinks  if  one  has  a  shorter  wavelength  light,  moves  closer  to 
the  obstacle  or  uses  a  larger  diameter  obstruction. 

For  a  given  obstacle  diameter  and  incident  wavelength,  the  spot 
of  Arago  (the  central  maximum)  can  be  thought  of  as  a  slowly  widening 
sheath  surrounding  the  optical  axis  for  increasing  distance  z.  This 
concept  has  an  intriguing  consequence.  If  one  observes  the  spot  at  a 
sufficient  distance,  the  diameter  of  the  spot  will  be  equal  to  or 
greater  than  the  diameter  of  the  shadow.  In  short,  the  shadow  will 
disappear.  If  the  preceding  theory  is  valid  in  this  range  then  this 
will  occur  when: 

kr2d  . 

— 2z  =  2.405  with  r2  _>  ~ 


z  > 


(III— 29) 


For  a  1cm  diameter  obstacle  illuminated  by  He-Ne  laser  light  of  6328a 
wavelength,  then  theory  predicts  shadow  saturation  when: 

z  _>  103m  (III-30) 

One  must  view  the  above  numbers  with  some  caution.  One  restriction 

2r  „ 

the  theory  was  based  on  was  — -  <<  1.  The  theory  may  not  be  valid  for 

d 

this  large  a  value  of  r„ . 


The  edge  diffracted  rays  define  the  disturbance  at  points  in  the 


shadow.  One  expects  a  reproduction  of  the  previous  results  multiplied 

c 

by  the  gaussian  factor  e  F  for  the  diffraction  of  a  wave  with  a 
gaussian  distribution. 

The  Ray leigh-Sommerfeld  integral  is: 


“  2n 


2ttA 

iAz 


ikz 

e 


(III— 3 1 ) 

The  evaluation  of  the  integral  follows  the  method  presented  in  the 
preceding  section  with  the  result  for  the  disturbance  function: 


I 


(III-32) 


The  intensity  distribution  is: 


k2d4  (III-33) 

2  2 

For  most  cases  of  interest  L  z  c -  <<  1  and  can  be  neglected.  (For  z 

k2d4 

on  the  order  of  1  meter,  0=10  and  d  on  the  order  of  .01m  then 
z2c2  -4 

-  — 7— r-  '  10  Then  the  intensity  distribution  reduces  to  the  antic- 

k  d 

ipated  form: 

_  £  /kdrz\ 

I(r2,  z)  =  A2  e  2  JQ2  \^~)  (III-34) 

The  resulting  intensity  on  axis  and  radial  intensity  distributions  are 

£ 

reduced  by  the  gaussian  factor  e“  2* 

The  following  graphs  demonstrate  the  intensity  distribution  on  and 
off  the  optical  axis  for  different  values  of  the  gaussian  constant  c. 


OH  AXIS  NORMALIZED  INTENSIIV  BEH I ND  A  1  CM  OBSTACLE:  INCIDENT  GAUSSIAN  NAMES 


DISTANCE  BEHIND  OBSTACLE  (CM) 

GAUSSIAN  CONSTANTS:  SOUL  r  0  DASHED  =  .25  DOTTED  :  1  SQUARE  r  2.5 


Graph  J. 1 1 


Graph  ;>.13 


Case  3:  Annular  Aperture  with  Uniform  Illumination 


The  geometrical  layout  for  diffraction  by  an  annular  aperture  is 
depicted  below.  The  Rayleigh-Sommerfeld  integral  is  given  by  III-7 


Figure  3-5 

Geometrical  Layout  for  Diffraction  by  an  Annular  Aperture 


substituting  —  in  for  the  upper  limit  for  the  integral  over  r^. 


D/2  2 


U(r2,  02,  z)  = 


*/..  /  4 


cos  (n,  1)  r  ,dr,d0, 

^  ^  III 


d/2  0 


(III-35) 


On  Axis  Intensitj 


The  arguments  follow  those  for  the  circular  obscuration  with  the 


final  result  for  the  disturbance  at  point  P  (from  I I I — 1 0 )  being: 


on  axis  intensity  distribution  is  thus: 


2  ?  z  z 

1(0,  0,  z)  «  U (0,  0,  z)  =  A  2  %  ■  +  - 

z  +D  /4  z  +d  /4 


x  (cos(k(  */z2+D2 /4  - 


z2+d2/4))  +  cos  (k(*/z2+d2/4  - 
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*  is  the  amplitude  for  the  disturbance  of  the  incident  wave  in  the 
aperture  plane.  The  coefficients  in  front  of  the  exponentials  are  the 
cosine  obliquity  factors  for  points  on  the  edges  of  the  aperture.  The 
exponential  terms  are  the  phase  factors  for  rays  reaching  from  the 
inside  and  outside  boundaries  (through  distances  ^  z2+d2/4  and 
z2  +D2 /4  respectively).  As  with  the  case  for  the  circular  obstacle 
the  edge  diffracted  rays  cause  the  disturbance  at  P2*  Those  for  the 
circular  obstacle  are  all  in  phase  for  axial  points  and  thus  interfere 
constructively.  For  the  annular  aperture  the  rays  from  the  inside  and 


outside  borders  will  interfere  to  varying  degrees  based  on  the  rela¬ 


tive  phases.  This  phase  difference  is  given  by 


For  large  z  the  phase  difference  will  vary  slowly  with  changes  in  z. 


Thus  one  has  a  disturbance  function  that  oscillates  between  -2A  and  2k, 
rapidly  close  behind  the  aperture.  The  intensity  oscillates  between  0 

and  4  times  the  incident  intensity. 

E2  d2 

In  the  region  where  — r  <<  1  and  * — ~  <<  1,  maxima  will  occur  on 

4zl 

axis  when  the  cosine  term  equals  -1.  This  happens  when  the  argument 
is  an  odd  multiple  of n. 


(D2  _  d2)  =  (2n  +  1 ) fr 


n  =  0,  1,  2, 


P2-d2 

4A(2n+l) 


(III-40) 


Minima  will  occur  when: 


P2-d2 

8mA 


m 


1,  2, 


(111-41) 


The  last  maximum  on  axis  occurs  when  n  =  0  or  when: 

D2-d2 

2  "  4A  (III-42) 

The  last  miminum  ’before  the  intensity  is  driven  to  zero  at  infinite 
z)  occurs  when  m  =  1  or: 

_  P2-d2 

2  8A  (III-43) 


The  last  maximum  defines  the  boundary  for  the  oscillatory  behavior  of 
the  intensity  on  axis.  The  distance  to  the  last  maximum  will  decrease 
for  increased  wavelength  or  decreased  factor  (D  -  d  ).  (D  -  cT  )  can 
be  written  as  D“ (1  -  e2)  where  e  is  the  obscuration  ratio,  the  ratio 
of  the  inner  to  outer  diameters.  Thus  a  smaller  obscuration  ratio 
(for  fixed  D)  will  decrease  the  distance  to  the  last  maximum  and  a 
smaller  outer  diameter  (for  fixed  e)  will  accomplish  the  same  thing. 

Below  are  graphs  depicting  the  intensity  on  axis  for  variations 
in  the  parameters  A,  D  and  e.  (Unless  specified,  the  wavelength  is 


assumed  to  be  6 323$.) 


Graph  J>.  1  5 
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a  circular  obstacle  results  from  interference  between  rays  diffracted 
from  the  edge  and  from  the  center  of  the  aperture. 

The  intensity  on  axis  is  given  by: 


For  the  circular  aperture  the  last  maximum  on  axis  occurs  when 

z  =  4A(2n+l)  (111-47) 

A  graph  of  the  intensity  on  axis  behind  a  1cm  aperture  appears  below. 


AXIAL  DISTANCE  (METERS) 

DIAMETER  OF  OBSCURATION  -  0.0000  M.  DIAMETER  OF  APERTURE  =  0,0108  K, 


Graph  5-21 


If  r2  =0  then  the  above  reduces  to  equation  III-3S*  If  D  becomes 


infinite  then  J 


(^) 


goes  to  zero  and  the  above  reduces  to  the 


result  for  the  circular  obscuration.  When  d  equals  zero  one  obtains 
the  intensity  function  for  diffraction  by  a  circular  aperture. 


I(r2,  z)  =  A 


[jo2(tt) + 1  -  2Jo  (ir) cos  (fe  ”2) 


(111-51) 


The  intensity  pattern  produced  by  an  annular  aperture  is  a  com¬ 
plicated  series  of  bright,  dim  and  dark  rings  of  various  radii.  The 
following  intensity  contour  plots  and  radial  intensity  profiles  show 
the  diffraction  patterns  for  a  .3555cm  diameter  circular  obstacle,  a 
1  cm  diameter  circular  aperture  and  their  resultant  annular  aperture 
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combination  at  .5  meters. 
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Graph  3.24 

Circular  Aperture  Intensity  Distribution 


LEGEND:  |  I  <  .01  lo  B  I  <  .05  lo  «  I  <  .5  la  I  <  !"B  !o 


0.495  K. 


I  >  .8  lo 


Graph  3.25 

Circular  Aperture  Pattern 


Beats  and  the  ring  structure .  The  radial  intensity  plot  for  the 
annular  aperture  of  obscuration  ratio  1/3  above  yields  zeroes  for 
intensity  at  the  radial  distances  depicted  in  the  table  below. 

Table  3-2 

Radial  Positions  of  Intensity  Zeroes  for  Annular  Aperture 
(D  »  .01,  c  -  1/3,  z  =  .495m) 

Radial  distance 
from  center  (mm) 

.031 

.175 

.269 

.362 

.456 

average  .09375mm  +  .0005 


Distance  to 
next  zero  (mm) 

.094 

.094 

.093 

.094 


This  very  regular  pattern  also  appears  in  the  spacing  of  minima  for 
the  circular  aperture.  A  table  of  the  radial  positions  of  these 


appear  below. 


Table  3-3 


Radial  Position  of  Intensity 

Minima  for  Circular  Aperture 

0 

11 

0 

z  =  .495m) 

Radial  distance 
from  center  (mm) 

Distance  to 
next  minima  (run) 

.009 

.061 

.070 

.064 

.134 

.061 

.195 

.064 

.259 

.062 

.321 

.063 

V>J 

CO 

.062 

.446 

.063 

.509 

average  .0625mm  .0012 


The  frequency  of  zeroes  for  the  annular  case  is: 

u  =  - - -  =  10666.7m'1 

ann.  9.375x10  5m 

The  frequency  of  minima  for  the  circular  aperture  case  is: 

v  .  =  - - - r-  =  16000m'1 

circ.  ,  oc  ,  -=> 

6.25x10  m 


(III-52) 


(II 1—5  3 ) 
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To  understand  the  above  ratio,  one  can  approximate  the  Bessel  function 
as  a  sinusoidal  function.  The  period  for  JQ  for  large  arguments 
approaches  2~.  The  amplitude  for  large  arguments  is  essentially 
constant  (very  slowly  decreasing).  The  Bessel  functions  in  equation 
I I 1—51  are  of  the  form: 


Jo(^~vir')  =  csin(2~v^r. '  i  =  1,  2 

n  a 

The  sinusoidal  approximations  are  thus: 


where  v,  =  — —  and  v  ~  -rr 

1  2 A z  2  2X 


J  (2irv.r  )  =  c  sin  (2cv.r  A  i  =  1,  2 
0  i  2  x2 


(111-55) 


v1  is  different  than  v,  and  the  combination  of  the  two  periodic  pat¬ 
terns  will  lead  to  beating.  This  beat  frequency  is  given  by  the 
, .  4 


expression: 


beat 


V1  V2 


(III-56) 


2\z 


(D  -  d) 


=  2X7(1-£) 


=  Vj  (1  -  e) 


(111-57) 


For  the  case  where  £=1/3  then: 


v 

beat 


v 

1 


1 


(III — 58 ) 


The  frequency,  v  ,  corresponds  to  the  frequency  of  intensity 

oea  o 

zeroes  radially  for  the  annular  aperture  case.  The  frequency  ; 
corresponds  to  the  frequency  of  minima  for  the  circular  aperture  case. 
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The  frequency  o  applies  to  the  circular  obscuration  case. 


v  . 
circ . 


v 

1 


D 

2\z 


_ .01m _ 

2  x  6.328  x  10~7m  x  .495m 


15962m-1 


(III-59) 


v 

ann. 


v 


beat 


D_-_d 

2Az 


(.01  -  .003333)m 
2  x  6.328  x  10  7m  x  .495m 


10642m-1 


(III-60) 

The  concept  of  beats  provides  insight  into  the  outlying  ring 

structure  of  the  diffraction  pattern  of  an  annular  aperture.  The 

ring  structure  will  oscillate  within  a  beat  envelope  whose  radial 

width  is  given  by  1  /  v  or  Ql..  The  intensity  pattern  is  the 

drm*  D— d 

square  of  the  beat  structure  from  the  disturbance  function.  The 
result  is  a  series  of  dark  bands  (corresponding  to  the  regions  around 
the  zeroes  of  the  beat  envelopes)  and  bright  bands  (corresponding  to 
the  regions  around  the  maxima  of  the  beat  envelopes)  which  contain  an 
underlying  dark  and  bright  ring  structure.  These  bands  will  be  wide 
for  large  obscuration  ratios.  (The  beat  envelope  width  is  inversely 
proportional  to  (l  -  c) ) •  For  e  very  close  to  1  the  beat  envelope  is 
very  broad  and  one  will  observe  a  ring  structure  whose  frequency  of 
intensity  zeroes  approaches  that  of  the  frequency  of  minima  for  the 
circular  aperture  and  circular  obstacle.  For  very  small  u  the  beat 
envelope  is  very  narrow.  (The  study  of  aberration  effects  or.  the 
diffraction  pattern  will  involve  changes  in  the  ring  structure.  A 
wide  beat  envelope  will  provide  a  better  opportunity  to  .bserve  devia¬ 
tions.  Thus  an  aperture  with  a  high  obscuration  ratio  is  better  for 


aberration  testing. 


Graph  3.29 


Annular  Aperture  Pattern 
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Graph  3*34 

Intensity  Distribution  for  e  =  .9 


Graph  3.35 

Annular  Aperture  Pattern  for  e  =  .9 
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Table  3.4 

Comparison  of  Measured  Envelope  Widths  to  Theory 


Graph 

which 

ments 

from 

measure- 

taken 

Average  measured 
distance  between 
beat  zeroes  (mm) 

Theoretical 

distance 

2  z 

D-d 

Percent 

Difference 

Measured 

V 

ann. 

(1-0 

V  . 

circ. 

Graph 

3-30 

.0784  +_  .0019 

.0779 

.6  % 

.795 

.8 

Graph 

3-32 

.1270  +.  .001 1 

.1271 

.08% 

.500 

.5 

Graph 

3.34 

.5330  +_  .0005 

.6333 

.05% 

.100 

.1 

Graph 

3-36 

•  9053  _+  .0009 

•  9052 

.11% 

.666 

.  667 

Graph 

3.33 

.665  +  .0005 

.679 

2.1% 

.102 

.1 

Tolerances  in  patterns  with  on-axis  shifts.  The  preceding  charts 

and  graphs  were  made  for  planes  at  positions  of  maxima  along  the 

optical  axis  behind  the  annular  apertures.  The  intensity  of  the 

central  spot  changes  with  distance  z.  Mathematically  this  stems  from 

the  cosine  interference  term  in  equation  III-50.  One  can  view  the 

argument  for  the  cosine,  4—  (D2  -  d2),  as  a  phase  term.  For  constant 

oz 

z,  D  and  d  this  phase  term  will  be  fixed.  This  phase  term  will  affect 
the  intensity  of  the  rings  and  their  positions.  A  null  on  axis  will 
push  ring  structure  out  from  the  center. 

Below  are  graphs  which  demonstrate  the  change  in  the  diffraction 
patterns  for  shifts  to  the  next  minima  on  axis  for  two  previously 
calculated  cases  for  c  =  1/3  and  e  *  .9. 
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The  diffraction  pattern  makes  its  most  dramatic  change  as  the 
central  spot  shifts  from  on-axis  maximum  to  on-axis  minimum.  This 
relates  to  a  change  in  the  argument  of  the  cosine  term  by  * .  Desig¬ 
nating  the  argument  by  a  one  has: 


Aa  = - — -  '.''2  -  d2)  Az 

8z2 


(keeping  k,  D  and  d  constant) 


Let  Aa  =  air 
Az  =  - 


.  Then : 

AaX  ? 
D2-d2 


(III-61) 


(111-62) 


Thus  the  necessary  shift  on  axis  by  a  phase  factor  a  n  is  quadratic 
in  the  initial  distance  z.  As  z  increases  a  greater  shift  on  axis  is 
necessary  to  produce  a  given  change  in  the  diffraction  pattern. 


Case  4_:  Annular  Aperture  Illuminated  with 
Gaussian  Intensity  Distribution 


The  treatment  of  an  annular  aperture  illuminated  by  a  plane  wave 

r£ 

with  a  gaussian  intensity  of  the  form  e  Dz  follows  that  already  done 

p 

for  the  circular  obscuration  with  the  substitution  of  —  for  the  upper 
limit  of  the  integral  in  equation  III— 31 .  The  disturbance  function 


becomes: 
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U(r2,  z) 


A  ikz  ik 

e  e  2z 


1-i 


2zc 

kD2 


£  ..  d£ 

e"  4  e1  8z  J 


BO- 


£  d£ 

e  4 


ikfi 

D2  6  82  J0 


© 


(III-63) 


The  intensity  function  is: 


I(r2,  z)  « 


1+4 


z2c2 


k2D4 


e  2  J, 


/kDr,\  £  d2  /kdr,\ 

far)  +e  2  d2  jo2(it)  • 


£ 

2e  4 


(III-64) 


2  2 

Z"C 


As  with  the  circular  obscuration  case,  the  term  4  1  for  small 

k^D 

z  and  can  be  neglected. 

On  axis  =  0  and  the  Bessel  functions  go  to  1.  The  intensity 
function  becomes: 


[■ 


£  £  d2 

1(0,  z)  =  A2  I  e"  2  +  e'  2  7 


-  2e  4 


[1+  7]  cos  (jjj  (D2  -  d2))  j 


(III-65) 


The  positions  of  maxima  on  axis  are  given  by  equation  III-40.  At  the 
maxima,  the  intensity  is: 


£  £  ck 

I  (0,  z  )  =  A2  I  e  4  +  e  4^' 
max  I  D 


[■ 


J 


(III— 66) 


The  position  of  minima  on  axis  are  given  by  equation  III-41.  At  the 
minima,  the  intensity  is: 


1(0, 


z  . 
mm 


)  =  A2 


( II 1—67 ) 


Thus  for  a  gaussi  .n  intensity  distribution  the  intensity  on  axis 
oscillates  between  the  two  limits  given  by  the  equations  above.  If  c 
is  not  zero  then  the  intensity  on  axis  never  reaches  zero.  Likewise 
the  intensity  at  the  maxima  is  always  less  than  four  times  the  inci¬ 
dent  intensity. 

Physically  the  above  result  stems  from  difference  in  amplitude  of 
the  diffracted  rays  from  the  inner  and  outer  edges.  The  amplitude  of 


light  rays  reaching  an  aperture  point  a  distance  r  from  the  center  are 

r2 
-c- 


diminished  by  the  gaussian  term  e  d2  .  The  amplitudes  of  rays  at  the 
outer  edge  are  decreased  more  than  amplitudes  of  rays  at  the  inner 
edge.  Therefore  the  destructive  interference  of  diffracted  rays  is 
not  total.  The  constructive  interference  of  diffracted  rays  is  re¬ 
duced  by  a  combination  of  the  gaussian  terms. 


The  following  graphs  show  the  intensity  on  axis  behind  an  annular 
aperture  of  obscuration  ratio  .5  for  gaussian  constants  c  of  .5.  2  and 


4 

l- 

Graph  5.42 


y  ^  «  4 

For  large  z  or  for  large  gaussian  constants  c  the  term  b  -  in 

k^D4 

the  denominator  in  equation  III-64  is  no  longer  negligible.  This 
drives  the  intensity  on  axis  to  zero  at  large  observing  distances. 

As  the  gaussian  constant  c  increases  the  maximum  intensity  on 
axis  decreases.  The  minimum  intensity  starts  at  zero  for  c  equal  to 
zero  and  increases  as  c  does.  But  it  reaches  a  maximum  value  beyond 
which  the  exponential  decrease  starts  to  drive  the  minimum  intensity 
back  to  zero.  For  fixed  obscuration  ratios  one  can  maximize  the 
minimum  intensity  on  axis  by  changing  the  gaussian  constant.  The 
following  graphs  depict  the  maximum  and  minimum  on  axis  intensities  as 
functions  of  gaussian  constants  for  obscuration  ratios  of  .1,  .2,  .4, 
.6  and  .8. 
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Solid  line:  e  =  .1 
Heavy  dashed  line:  e  =  .2 
Dotted  line:  e  =  .4 


Squared  line: 
Thin  dashed  line: 


Graph  3.43 


MINIMUM  ON  AXIS  INTENSITIES  FOR  OBSCURATION  RATIOS  .1  .2  .4 


Case  5_:  Displaced  Obstacle  and  Aperture 
Illuminated  with  Uniform  Intensity 

If  the  obstacle  and  aperture  do  not  lie  in  the  same  plane,  it  is 
still  possible  to  derive  expressions  for  the  on  and  off  axis  intensity 
distribution  functions.  The  layout  depicted  below  cor  ponds  to  the 
eventual  experimental  set-up.  As  with  the  previous  tr  ment  the 
origin  of  the  cylindrical  coordinate  system  is  taken  t*.  h  che  center 
of  the  aperture. 

aperture 


Geometrical  Layout  for  Diffraction  by 
Displaced  Aperture  and  Obstacle 


The  observation  plane  lies  a  distance  z  behind  the  aperture  and  a 

distance  z  -  z  *  Lz  behind  the  obstacle.  The  rays  diffracted  by  the 
2  1 

obstacle  which  converge  at  point  P  project  an  effective  obstacle  in 
the  aperture  plane.  By  similar  triangles,  the  diameter  of  this  phan- 
tom  obstacle  is  d'  =  — d.  The  light  rays  striking  the  aperture  are 


parallel  to  the  optical  axis,  but  those  striking  the  effective 


To  evaluate  the  Rayleigh-Sommerfeld  integral  for  the  effective 
annular  aperture,  the  following  assumptions  are  made: 


«  1, 


2.  The  effect  from  rays  diffracted  by  the  obstacle  striking  the 
aperture  is  negligible.  The  rays  striking  the  aperture  are  considered 
parallel  to  the  optical  axis. 

3.  As  with  previous  treatments  only  those  rays  diffracted  by  the 
outer  aperture  edge  and  the  inner  effective  edge  contribute  to  the 
ultimate  diffraction  pattern. 

4.  The  rays  hitting  the  inside  and  outside  edges  of  the  effec¬ 
tive  annular  aperture  have  a  phase  difference  given  by  equation  III- 


The  Rayleigh-Sommerfeld  integral  for  the  disturbance  at  a  point  P 
in  the  observation  plane  is: 


,  2 it  D/2 


U(r2,  Zj,  z2)  iX 


W  - 


cos  (n,  1)  6  (rx)  r^drj 


where : 


(HI-71) 


B(rx)  =  1 


for  =  y 


sik  a  2  u2  -  v  £or  ri.  jii 


(III— 72) 


The  result  follows  from  previous  derivations: 
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Thus  the  disturbance  function  has  the  reassuring  form  of  having  the 
effect  from  the  aperture  edge  at  distance  Z} ,  coupled  with  the  effect 
from  the  obstacle  edge  at  distance  The  intensity  distribution 

function  is: 


81 


Notes 


^  Joseph  W.  Goodman,  Introduction  to  Fourier  Optics  (San  Fran¬ 
cisco:  McGraw-Hill ,  1968),  p.  44* 

p 

James  E.  Harvey  and  James  S*  Forghan,  "The  Spot  of  Arago:  New 
Relevance  for  an  Old  Phenomenon,"  American  Jou ~nal  of  Physics,  (Janu¬ 
ary  1984). 

^  Max  Born  and  Emil  Wolf,  Principles  of  - ptics,  6th  ed.  (Oxford: 
Pergamon,  1980),  p.  395* 

^  Harvey  E.  White,  Modern  Co  1  lege  Physics,  6th  e:.  (New  York:  Van 
Nostrand  Reinhold,  1972),  p.  304* 


hapter  4 


THEORY  I.  THE  ABERRATED  CASES 

General  Aberration  Considerations 

Wavefront  aberrations  are  deviations  from  spherical  wavefront 
behavior.  Consider  an  arbitrary  wavefront  in  space.  (A  wavefront  is 
defined  as  a  surface  of  constant  phase.)  It  is  convenient  to  project 
this  wavefront  onto  a  spherical  wavefront.  This  reference  wavefront 
shares  the  phase  and  optical  axis  point  with  the  arbitrary  wavefront. 
The  radius  of  curvature  is  the  distance  to  the  source  if  a  diverging 
wavefront  or  the  distance  to  the  focal  point  if  a  converging  wave- 
front.  The  projections  onto  the  reference  wavefront  reflect  phase 


Figure  4.1 


Aberrated  Wavefront 
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differences  due  to  the  optical  path  differences  between  corresponding 
points  on  the  two  wavefronts.  The  phase  differences  can  be  charac¬ 
terized  by  an  aberration  function  W  where  W(r,  r,  0)  is  the  optical 
path  difference  between  the  point  (r,  '*)  on  the  reference  wave- 

front  and  its  complement  on  the  arbitrary  wavefront.  Thus  the  phase 
differences  appear  in  the  form: 

ikW 

e 


First  Order  and  Seidel  Aberrations 

The  aberration  function  can  be  expanded  in  terms  of  a  complete 

set  of  polynomials  called  the  Zernicke  polynomials.1  In  many  cases  it 

is  sufficient  to  expand  the  aberration  function  in  terms  of  primary 

aberrations-  Convenient  expressions  for  these  primary  aberrations 

were  made  by  Seidel  and  the  third  order  aberrations  now  bear  his  name. 

The  wavefront  aberration  representations  appear  in  Table  4.1.  The 

coefficients  Wa^c  are  constants  with  the  dimension  of  length.  0  is 

the  polar  coordinate  angle  of  points  in  the  reference  plane,  r  is  the 

radial  distance  in  the  reference  plane  normalized  to  a  given  pupil 

radius.  When  the  pupil  is  the  outer  edge  of  an  annular  aperture  r  = 
r! 

Similarly,  B  is  a  normalized  field  parameter.  When  the  source 
is  a  point  source  £  is  a  constant  and  can  be  absorbed  into  the  coeffi¬ 
cients  W  .  Table  4.2  lists  field  independent  wavefront  aberrations. 

at  c 

In  the  field  independent  situation  field  curvature  acts  like  defocus 


and  distortion  like  tilt. 


.  . ..  ,-t-.  l-  ^»  »  ^  tj  v;  r. 


84 


Table  4.1 

The  Primary  Wavefront  Aberrations 


Aberration  type 

Representation 

Defocus 

W 

020r 

First  order 

Tilt 

WlnBr  cos0 

aberrations 

Spherical  aberration 

W  r4 
"04  or 

Coma 

W131&r3cos0 

Third  order 

Astigmatism 

W222B2r2cos20 

(Seidel) 

Field  curvature 

aberrations 

Distortion 

W3n83r  cos0 

Table  4.2 

Field  Independent  Wavefront 

Aberrations 

Aberration  type 

Representation 

Defocus  /  Field  Curvature 

4W2  0rl 2/1)2 

Tilt  /  Distortion 

2Wj  jrjCos0/D 

Spherical  aberration 

lew^i-ZD- 

Coma 

8W31 r j  3cos0/D3 

Astigmatism 

4W22r j2c°s20/D2 

Aberrations  Produced  by  Plane  Parallel  Plates 

Consider  a  diverging  wavefront  incident  upon  a  tilted  plane 
parallel  plate.  The  tilt  angle  is  given  by  6  and  u  is  the  half  angle 


of  the  cone  of  illumination,  u  is  related  to  the  f-number  by  the 
2 

formula: 

F#  =  ^  (IV-1 ) 

The  thickness  of  the  plate  is  t.  It  is  clear  to  see  that  phase 
differences  result  from  passage  through  the  plate  by  breaking  up  the 
wavefront  into  rays.  Rays  corresponding  to  different  f-numbers  will 
strike  the  plate  at  different  angles.  By  Snell's  law  the  refraction 
angles  will  not  be  the  same  and  the  path  lengths  through  the  plate 
will  differ.  Because  light  moves  more  slowly  in  glass  than  in  air, 
transit  times  vary  and  phase  differences  ensue.  The  phase  differences 
appear  as  aberrations  in  the  wavefront.  A  table  of  the  primary  wave- 
front  aberrations  coefficients  as  functions  of  n,  u,  3  and  t  produced 


by  a  plane  parallel  plate  appears  below:^ 


Table  4.3 


Aberration  Coefficients  Due  to  Tilted 
Plane  Parallel  Plate 


Aberration  Coefficient 

Expression 

W 

22 

(Astigmatism) 

1 

2 

n2-l 

n3 

tu2e2 

W 

(Coma) 

1 

n2-l 

tu3e 

31 

2 

n3 

w 

(Spherical  Aberration) 

__  1 

n2-l 

tu4 

40 

8 

n3 

Notice  the  coefficient  for  coma  is  an  odd  function  of  the  tilt 
angle  0.  By  countertilting  two  plane  parallel  plates  of  equal  indices 
of  refraction  and  thickness  comatic  effects  can  be  neutralized.  At 
the  same  time  the  coefficients  for  astigmatism  and  spherical  aberra¬ 
tion  double. 

Table  4-4 

Aberration  Coefficients  for  Two  Countertilted 
Plane  Parallel  Plates 
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Typically  the  aberration  coefficients  are  very  small  numbers,  having  the 
magnitude  of  order  of  optical  wavelengths.  It  is  convenient  to  express 
the  aberration  coefficients  in  terns  of  wavelengths  of  the  incident 
light.  The  two  graphs  below  illustrate  the  aberration  coefficient 
curves  for  a  single  plate  and  two  counter  tilted  plates.  The  incident 
wavelength  is  taken  as  6328A,  the  thickness  of  one  plate  is  12.7mm  (1/2 
inch),  the  index  of  refraction  n  =  1.517  and  the  angle  u  =  .1  (5.74  ). 
The  aberration  coefficients  are  plotted  as  functions  of  the  plate  tilt 
angle  0., 


U 

A 

U 

E 


N 

C 

T 

H 


S 


ABERRATION  COEFFICIENTS  PROBUCEB  B'i  SINGLE  PLANF  PARALLEL  PLATE 


0  4  8  12 

TILT  ANGLE  (DECREES) 

SOLID  LINE  =  W22  DASHED  LINE  =  H46  DOTTED  LINE  =  W31 


Graph  4.1 


u 

A 

V 

E 

L 

E 

N 

G 

T 

H 

S 


ABO] 

S.0 

4.0 

3.0 

2.0 

1.6 

WIT  I  ON  C01 

EFFICIENTS 

1  PRODUCE! 

' 

i  BV  COUNT 

! 

ERTILTED 

i 

PLANE  PAR 

ALLEL  PLATES 

1 

1 

y 

i 

_ ! 

! 

/  ! 

1 

I 

i 

1 

| 

,y 

\ 

l 

I 

l _ 1 

y _ 

!  y 

i _ i _ 

s 

1 

[  !  -  -  | 

I 

i _ i _ 

i 

l _ ! _ 

i 

1 

L.^-  i  i 

t  ] 

i  i  i  i  ; 

1 - ■ - 1  -T  1 

l _ i _ : _ ! _ i _ : _ : _ 

0  4  8  12 

TILT  ANGLE  (DEGREES) 

SOLID  LINE  =  H22  DASHED  LINE  =  H40 


Graph  4.2 


Symmetries  in  the  Primary  Aberrations 

The  primary  aberrations  have  characteristic  symmetries  which  will 
ultimately  determine  their  contributions  to  the  diffraction  pattern. 

Both  defocus  and  spherical  aberration  are  independent  of  6.  They  are 
rotationally  symmetric  and  as  such  should  have  little  or  no  effect  on 
the  rotationally  symmetric  diffraction  patterns  discussed  in  Chapter  3. 
Tilt  and  coma  are  products  of  odd  functions  in  r  and  even  functions  in 
6.  Because  cos0  =  cos  (-0),  the  aberration  terms  for  tilt  and  coma  will 
have  equal  values  for  fixed  r  and  0.  This  means  there  is  one  plane  of 
symmetry  for  these  terms.  Astigmatism  is  the  product  of  an  even  func¬ 
tion  in  r  and  the  even  function  cos-' 0.  Because  costs')  =  cob^  (tt  -  0 )  = 
cos“(m  ♦  0)  =  cos^CPtt  -  0)  astigmatism  exhibits  two  planes  of  symmetry. 


Each  primary  aberration  will  be  discussed  separately. 


Defocus 


Defocus  can  best  be  discribed  as  a  shift  of  the  focal  point  on 
the  optical  axis.  In  effect  the  reference  spherical  wavefront  is 
changed  as  the  radius  of  curvature  changes.  If  the  focal  point  shift 
on  optical  axis  is  designated  as  Az  then  the  coefficient  of  defocus 
for  a  point  on  the  reference  spherical  wavefront  with  f-number  is:^ 


W20 


Az 

8(F#)2 


(IV-2) 


Case  6:  Circular  Obscuration 


The  Rayleigh-Sommerfeld  integral  of  interest  is: 


2tt 


U<r2,  t2.  z>  elkr22/2z  [  J  eikS  eik4W20jI 


0  d/2 


rlr2 

-ik  -  cos  (0.  -  0.) 

x  e  z  21 


2~A  lkz  lkr  2/2z 
-rr —  e  e  2 

lA 


d/2 


/krlr2\ 

x  Jo  l— ~ J  ridri 


(IV-3) 


The  solution  of  the  above  integral  follows  the  procedures  set  forth  in 


Chapter 


,,,  v  kA  zd2 

U(r  ,  z)  =  — - 

2  1AZ  ik(d2+SW20z) 


;ik2  elk  ~2z  eik  L  27  +  ~ J  T 


x  J, 


(5) 


-A 


ik(z 


l  +  3W20z 


+  £  +  K*°>  Jo 


(IV-4) 


The  intensity  distribution  is: 


I(r2,  2)  = 


(8W  z\2  0 


(*) 


(IV-5) 


If  no  defocus  is  present  W  =  0  and  the  above  expression  becomes  equiv 
alent  to  the  unaberrated  case  in  equation  III-26.  The  Bessel  function 
term  is  identical  to  that  for  the  unaberrated  case  —  as  expected  this 
rotationally  symmetric  aberration  does  not  affect  the  rotationally  sym¬ 
metric  diffraction  pattern.  Where  defocus  does  enter  is  in  the  overall 
intensity.  For  positive  defocus  the  factor  1  /  1 1  +  — ~~2~J  decreases 

as  the  distance  along  the  optical  axis  z  increases.  The  effect  of 
positive  defocus  is  to  decrease  the  overall  intensity  for  large  z. 
Something  interesting  happens  when  the  coefficient  of  defocus  is  nega¬ 
tive.  At  one  value  of  z  the  intensity  on  axis  goes  to  infinity.  This 
z  corresponds  the  focal  point  of  the  system.  To  see  this,  the  denom- 


inater  goes  to  zero  when: 
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The  obscuration  (or  aperture)  will  act  as  a  focusing  system.  With  no 
aberration  present  it  would  obey  the  lens  law. 


(1V-7) 


/■lere  s  is  the  distance  to  the  source,  i  is  the  distance  to  the  image 
and  f  is  the  focal  distance  of  the  system.  With  defocus  present  the 
lens  law  is  replaced  by:^ 


(IV-8) 


(IV-9) 


With  no  aberration  present  the  terms  s,  i  and  f  in  equation  IV-7  are 
all  infinite.  With  defocus  s  and  i  remain  infinite  and  f  changes. 
Then: 


(IV-10) 


Thus  the  focal  point  of  the  system  is  the  singularity  on  the  optical 
axis. 

As  the  magnitude  of  defocus  is  increased,  the  focal  point  moves 
in  from  infinity  and  approaches  the  aperture.  For  negative  defocus 
the  on  axis  intensity  will  increase  until  the  focal  point  is  reached. 


Then  it  will  decrease.  (If  the  coefficient  of  defocus  is  positive  the 
system  focal  point  occurs  at  negative  z  and  the  decreasing  intensity 
already  discussed  is  an  extension  of  the  decrease  from  the  focal 
point. ) 

The  singularity  intensity  is  reminiscent  of  the  geometrical  op¬ 
tics  limit  where  the  intensity  at  the  focal  point  is  infinite.  Appen¬ 
dix  5  demonstrates  the  infinite  intensity  is  a  result  of  an  infinite 
plane  wave  striking  a  circular  obscuration.  Physically  the  extent  of 

the  incident  wave  is  always  finite  so  the  intensity  at  z  =  -  tt —  will 

oV,2o 

never  go  to  infinity.  The  graphs  below  demonstrate  these  effects  for  on 
and  off  axis  intensity  plots. 


Graph  4*5 
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Graph  4.6 


ON  AXIS  INTENSITY  BEHIND  ANNULAR  APERTURE:  SPHERICAL  ABERRATION  EXACT  SOLUTION 


AXIAL  DISTANCE  CKETERS) 

OUTER  DIAMETER  =  .01  M  INNER  DIAMETER  =  .005  M  N40  ;  -.5  WAVELENGTHS 


Graph  4.7 


The  on  axis  intensity  is  characteristic  of  the  defocus  present. 
By  measuring  the  intensity  a  given  distance  z  along  the  optical  axis 
one  can  readily  determine  the  amount  of  defocus  in  the  incident  beam. 
From  equation  IV-5: 


(IV-1I) 


Case  7_:  Annular  Aperture 

In  the  case  of  the  annular  aperture  the  aberration  function  is 
normalized  to  the  aperture  outer  radius.  The  expression  for  W  is: 

r  2 

W  =  4W  — 

20  2 

D  ( IV- 12) 

Substituting  the  above  into  equation  IV-5  and  changing  the  upper 
integration  limit  for  r^  to  D/2  allows  one  to  write  out  the  distur¬ 
bance  function  from  equation  IV-4. 
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The  intensity  distribution  is: 


I(r2,  z)  = 


I  T  2 

Jo 


(■*  l 


*vf&)  -.(£) 


/kdr, 


X  J0  ^  2z 


|  cos  | 


8z 


(D2  -  d2)  +  kl 


■J'-Sl) 


(IV-14) 


The  cosine  terra  can  be  rewritten  in  a  somewhat  more  condensed  form: 


0 


cos  |  —  (D2  -  d2)  +  kW 

oz  20 


[1'l])'c°s(fe  <d2‘ 


d2) 


[■•?]) 


(IV- 15) 


When  no  defocus  is  present  the  intensity  distribution  reduces  to  that 
for  the  unaberrated  case. 

Defocus  will  affect  the  diffraction  pattern  of  an  annular  aper¬ 
ture  in  two  ways.  The  first  is  in  the  change  of  the  overall  intensi¬ 
ty.  The  intensity  distribution  is  multiplied  by  the  same  term  1  / 

(l  +  8W2 q z/D“ ) 'as  for  the  circular  obscuration  case.  The  coefficient 
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for  defocus  also  enters  the  cosine  term.  Thus  the  second  effect 
defocus  has  is  to  shift  the  positions  of  maxima  and  minima  on  axis. 

Maxima  will  occur  when  the  argument  of  the  cosine  term  equals  an 
odd  multiple  of  n. 


(2n  +  1)ti 


_ kP2(P2-d2) _ 

8m(2n+l)D2  -  8W20k(D2-d2) 


n  =  0,  1 ,  2 ,  .  .  . 


_ D2 ( 1-e2) _ 

4X(2n+l)  -  8W20(l-e2) 


(IV-16) 


Minima  occur  for: 


P2 ( 1-e2) 

z  “  - 

8nX  -  8W20(l-e2) 

The  last  maximum  on  axis  occurs  when  n  =  0  at  a  distance 
P2(l-e2) 

z  *  - - 

4X  -  8W20(1-e2) 


(IV- 17) 


(IV- 18) 


The  above  formula  suggests  the  last  maximum  on  axis  will  occur  at  infi¬ 
nite  z  when  the  denorainater  goes  to  zero,  or  when: 


W20 


X 


2 (1-e2) 


(IV- 19) 


9b 


But  there  is  a  second  effect  to  consider.  The  coefficient  1/(1  + 
8W2Qz/D2) ^  goes  to  zero  for  infinite  z.  This  term  is  the  upper  limit 
for  the  on-axis  intensity.  The  intensity  oscillates  between  0  and 
this  upper  bound.  The  result  that  the  last  maximum  on  axis  occurs  at 
infinite  z  is  a  statement  that  there  is  no  last  maximum  on  axis  behind 
the  aperture  —  the  oscillations  continue  to  infinity.  In  general 
this  will  be  the  case  whenever: 

8W20  (1  -  e2)  >  4X  (2n  +  1) 


X 


( IV— 2 1 


-  (D2  -  d2)  dz  =  2n 

8z2 


D2-d2  (1-e2)  D2 


The  following  graphs  demonstrate  the  on  axis  intensities  behind  a  one 
centimeter  outer  diameter  annular  aperture  with  e  *  .5  and  with  de¬ 
focus  coefficients  of  .5X*  -«5X,  .667\  and  2A.  (.667A  corresponds  to 

the  case  where  the  last  maximum  occurs  at  infinity.)  The  envelope 
1/(1+  8W2qz/d2)2  is  also  plotted. 


ON  AXIS  INTERS ITV  BEHIND  A  1  CM  OUTER  DIAMETER  ANNULAR  APERTURE 
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ON  AXIS  INTENSITY  BEHIND  A  1  CM  O'JTEF  DIAMETER  ANNULAR  APERTURE 


AXIAL  DISTANCE  (NET EPS) 

OUTER  DIAMETER  =  .01  M  INNER  DIAMETER  =  .BBS  M  N2B  =  2  WAVELENGTHS 


Graph  4.11 


The  annular  aperture  provides  two  methods  of  determining  the 
coefficient  of  defocus.  The  first  is  to  seek  out  an  on  axis  maximum 
and  measure  the  intensity.  This  intensity  must  lie  on  the  curve; 


(1V-22) 

The  amount  of  defocus  present  can  be  found  in  the  same  way  discussed 
for  the  circular  obscuration.  The  previous  graphs  show  the  limita¬ 
tions  of  measuring  the  intensity  at  a  maximum  to  find  the  coefficient 
for  defocus.  The  curve  defining  the  upper  bound  for  the  oscillations 
is  very  steep  around  the  focal  point.  This  curve  is  tangent  to  the 
on-axis  intensity  distribution  and  the  actual  intensity  maximum  does 
not  necessarily  coincide  with  this  tangent  point.  For  example,  from 
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Oraph  4*9  for  minus  one  half  wavelength  of  defocus  the  last  maximum  on 
axis  as  given  by  equation  IV-18  occurs  at: 


-PlL1— L2) —  =  16.93m 


4X  -  8W2o(l-e2) 


(IV-23) 


The  actual  (measured)  intensity  maximum  occurs  at  z  =  20ra  with  a  nor¬ 
malized  intensity  of  1 3* 6.  Putting  these  values  into  equation  IV-'  1 
yields  a  calculated  value  for  W.n: 


u  = 

20  8z 


-II  =  -4.6  x  10  7m  =  - . 72a 


(IV-24) 


—  a  44^  difference  from  the  actual  value.  When  using  this  technique 
to  find  the  coefficient  for  defocus  one  should  take  measurements  at 
values  of  z  far  away  from  the  effective  focal  point  of  the  system. 

The  second  method  is  to  measure  the  shift  of  extrema  on  axis.  With  no 
defocus  the  n^  maximum  occurs  at: 


=  P2n-e2) 

20  4X(2n+l) 


(IV-25) 


With  defocus  it  occurs  at: 


P2(l-e2) _ 

4X(2n+l)  -  8W20(l-e2) 


(IV-26) 


The  shift  on  axis  is  then: 


z  -  zo  = 


=  d2  (i-t2)  r 

"  4X (2n+l)  I 


4X(2n+l) _ 

4X(2n+l)  -  8W20(1-e2) 


-  1 


2W20 (1-e2) 

X(2n+1)  -  2W20 (1-e2) 


(IV-27) 


Az 

z’ 


(IV-28) 


W 


20 


(2n+l)A 
2 ( 1-e2) 


This  method  will  be  effective  if  the  maximum  of  order  n  for  the 
aberrated  case  can  be  determined.  A  good  criterion  for  this  is  if  the 
shift  on  axis  is  less  than  the  spacing  to  the  next  maximum.  Comparing 
equations  IV-21  and  IV-25,  one  must  satisfy  (W20  in  this  case  is 
negative) : 


2K20(l-e2) 


A(2n+1)  -  2W20(l-e2) 


8 

(1-e2) 


A 


W2Q(1  -  e2) 


>  d-e2)  .  , 

2U20  (2n+l)  X 


<  _ 

(1-e2) 


2n+l 

2n-l 


(IV-29) 


2n+-j_ 

For  large  n  - — —  =  1  and  the  coefficient  for  defocus  should  be  less 

2n-l 


than  half  the  minimum  value  in  equation  I V  —  1 9* 


Tilt 


When  tilt  is  present  in  the  incident  beam  the  center  of  the 
diffraction  pattern  moves  off  the  optical  axis  without  changing  its 
functional  form.  Tilt  will  not  be  investigated  in  the  experimental 


Spherical  Aberration 


Spherical  aberration  is  a  rotationally  symmetric  aberration.  As 
with  defocus  the  presence  of  spherical  aberration  should  not  a. feet 
the  overall  diffraction  pattern.  However  the  intensities  and  p:si- 
tions  of  on-axis  extrema  should  be  influenced. 


Case  S:  Circular  Obscuration 


The  Ray leigh-Sonmerfeld  integral  for  this  case  is: 

2  TT  oo  2 


.  w  x  A  ikz  ikr  2/2z 

U<r2’  z)  "TSe  e  2 


I L  •“ * 


ikl6W 


e  40  d4 


'0  d/2 


rir2 

-ik -  cos  ($-$,) 

x  e  z  2  1 


kA  ikz  ikr 
=  —  e  e  2 
iz 


r°°  r  2 

2/22  j  +  16\o-7Jr, 


d/2 


x  J, 


(IV-30) 


The  evaluation  of  the  integral  requires  several  steps.  Consider  the 
function: 


.  z  ikr ,  2  T  +  1 6K  n  — —  1 
x)  —  e  1  l2z  40  dJ 


g'O 


/ 

'  32W40rj2z  > 

1  ‘kIr 

r1  +  32“-r‘22] 

■  ( 

1  +  - 

'  d4  > 

)  V  2z 

L  d4  J 

( IV— 3 1 ) 


Then  rewrite  the  original  integral  as: 


(IV-32) 


The  first  integral  on  the  left  side  is  solved  as  in  Chapter  3  with  the 
solution: 


(IV-33) 


Consider  the  second  integral.  Let  g(rj)  be  the  same  as  in  IV —51  and 

/  kr  r  \ 

f(r  )  =  r  2Jn  (  — —  ) 

1  1  0\  z  j 


(IV-34) 


Then: 


C  ri 2  r  32Wunr  ,2zl 


/  If  v  V 


32W4Qz 

g'(rj)  f(rx)  drj  - -5 — 

d 


ik  -r — 
e  2z 


32W40r12z 


1  z 

_  f(ri)  ri3dri 


(IV-35) 


The  diffraction  effects  come  from  the  edge  where  r  =  d/2.  Making  the 
assumption  f(r  )  oscillates  slowly  compared  to  g'^  )  the  first  inte¬ 
gral  is  solved  by  the  method  of  Chapter  3- 


g’Cr^ 


,  .  z  d2  lk£^  [l+S^]  ,  /  kdr2 ) 

>  dri  ■  sr e  8z  L  a2  J 


(IV-36) 


Grouping  the  terms  thus  evaluated  yields  for  the  original  integral: 


d/2 

/  kdr  \  d2  [  8W4  0z  ~\  / 

-sJo(Tr)-lk^Ll+“J  (1- 


8W40Z 


'  f  d/2 

_  /krir2\  5^ 

Jo\^— ■ ry  ri  dri 


(IV-37) 


If  the  last  integral  is  evaluated  using  the  function: 


(IV-38) 


one  obtains  the  next  term  in  the  series  (1  -  (8W4gz/d 2)  +  .  .  .  ) 
which  is: 


(IV-39) 


The  expected  solution  to  the  original  integral  is: 


a.  j 

ik  0  \  2z  / 
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(IV-40) 


for 


8W  z 

“72 — 


<  1 


Notice  this  implies  z  < 


O'.' 

'-'•‘4  0 


z  =  -  ■qtt"  is  the  effective  focal 

O»40 


point  of  the  system.  The  disturbance  function  is  therefore: 


The  intensity  function  is: 


I(r2,  z) 


1  oo 


(IV-42) 


The  intensity  distribution  behind  a  circular  obscuration  has  an 
identical  form  spherical  aberration  and  defocus.  In  fact  if  both 
defocus  and  spherical  aberration  are  present  the  intensity  distribu¬ 
tion  is: 


Kr,,.  z)  = 


(IV-43) 


where  W'  =  W  +  W 

20  40 


(IV-44) 


All  previously  derived  results  for  defocus  in  this  case  hold  for 

*> 

spherical  aberration  for  the  region  of  validity  z  <  •  This 


Di 

UA4  0 


region  of  validity  includes  normal  laboratory  conditions  for  typical 
amounts  of  spherical  aberration.  If  an  upper  limit  of  W4Q  =  1 5 X  |  is 
assumed  for  a  one  centimeter  diameter  obscuration  then  equation  IV-42 
will  hold  for  z  <  4  meters. 


Case  9_:  Annular  Aperture 

The  aberration  function  for  spherical  aberration  for  the  annular 
aperture  case  is: 

V 

W  =  I6W40 

D1*  (IV-45) 

The  disturbance  function  is  found  from  expression  IV-Jb  with  the 
upper  limit  evaluation  of  r^  =  -y. 


U(r, 


2)  =  -Ae 
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The  intensity  distribution  is: 


x  (D2  -  d2)  +  kW4u 


- 


(IV-47) 


Intensity  on  axis.  The  intensity  on  axis  is  described  by  the 


following: 
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The  intensities  for  the  maxima  and  minima  lie  on  the  curves 


1(0,  z  )  =  A2 
max 


L  1+ 


8w,o2 

D2 


1+ 


8W4Qzd2 


1(0,  z  .  )  -  A2 
mm 
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8Wu02 


8W4Qzd2 
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(IV-49) 


As  with  the  case  of  the  annular  aperture  illuminated  by  a  plane  wave 
with  a  gaussian  intensity  distribution  the  on-axis  intensity  distribu¬ 
tion  oscillates  between  an  upper  and  lower  bound  with  the  lower  bound 
non-zero  (except  for  z  =  0  and  z  =  «■)•  For  positive  W  the  upper 
bound  is  a  decreasing  function  with  increasing  z.  The  lower  bound 
will  increase  to  a  maximum  value  then  go  to  zero  for  large  z.  This 
maximum  occurs  when: 


dl  . 

=  o  =  2AZ 
dz 
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D2 


1+ 


_1 _ 

8W4Qzd  2 
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1 1 1 


D3 

8W40d  (IV-50) 


Since  ^  >  1 


D3  >  D2 

8WHOd  8W4  0  (IV-51) 


Thus  the  maximum  lies  beyond  the  limit  where  equation  IV-43  is  valid. 
In  the  region  where  the  equation  holds  the  minimum  intensity  boundary 
curve  is  a  constantly  increasing  function. 

The  amount  of  spherical  aberration  produced  by  plane  parallel 
plates  is  very  small.  In  addition  the  experimental  phase  will  involve 
distances  z  on  the  order  of  one  meter.  In  all  cases  of  concern  the 
following  approximations  are  valid: 


8W4  0Z 
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8W  zd2 
40 


«  1 


(IV-52) 


With  these  restrictions  the  denominators  in  equation  IV-48  are 
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replaced  by  their  expansions  to  first  order.  The  intensity  on  axis 
then  is: 


1(0,  z) 


If  one  measures  the  intensity  at  a  maximum  on  axis  the  spherical 
aberration  coefficient  is  found  to  be: 


W 


40 


D2 

32z ( 1+e2) 
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(III-54) 


As  with  the  defocus  case,  spherical  aberration  also  causes  an 
axial  shift  of  maxima  and  minima.  Maxima  will  occur  when: 

(D2  _  d2)  +  kW^o  (1  _  e4)  =  (2n  +  1)lT 


D2(l-e2) 

z  =  - L - 

4A(2n+l)  -  8W40(l-e^) 


(IV-55) 


Minima  occur  when 


z 


D2 ( 1-e2) 

8mA  -  8W40(1  -e1*) 


(IV-56) 


The  spacing  on  axis  from  the  nth  order  maximum  to  the  n  +  1 order 
maximum  is: 


hz 


4  z2 
(1-e2)  D2 


X 


(IV-57) 


As  before  the  spacing  is  independent  of  aberration.  The  coefficient 
for  spherical  aberration  can  be  calculated  from  the  axial  shift  of 
maxima  on  axis.  Following  the  argument  in  the  section  on  defocus: 
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OH  AXIS  INTENSITY  BEHIND  ANNULAR  APERTURE:  SPHERICAL  ABERRATION  EXACT  SOLUTION 


0  .5  1  1.5  2  2.5  3  3.5  4 

AXIAL  DISTANCE  (KETFRS) 

OUTER  DIAHETEE  =  .81  M  INNER  DIAMETER  :  .825  M  W40  =  ,5  ViAl'ELENC-THS 


Graph  4. 1 2 
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IHTENSm  BEHIND  ANNULAR  APERTURE:  SPHERICAL  ABERRATION  APPROBATE  SOLUTION 


AXIAL  DISTANCE  (MEIERS) 

OUTER  DIAMETER  =  .Pi  M  INTO  DIAMETER  =  .005  N  W40  r  .5  WAVELENGTHS 


Graph  4.13 
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RADIAL  DISTANCE  <.ePi  M) 

OBSTACLE  DIAMETER  =  .01  M.  2  DIST.  =  2.00  M. 


W20  :  0.5  WAVELENGTHS 


Graph  4.14 


Graph  4.15 


Off  axis  intensity  distribution.  Off  axis,  the  intensity  behind 
the  annular  aperture  is: 


The  following  graphs 
with  WUQ  =  .25-\,  2.5^  and 


show  the  patterns  for  an  on  axis  maximum 
-2.5^.  There  is  very  little  change  in  the 


patterns. 
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RADIAL  DISTANCE  (.081  N) 

O.D.:  0.01  M.  I.D.  :  .003  M.  Z  DISI.  :  .502  K.  W40 


:  -2.5  WAVELENGTHS 


Graph  4.18 


Refocus  and  Spherical  Aberration  Combined 


8VJ  z 

U  Q 

Earlier  the  restriction  -  <<  1  was  imposed.  A  resonable 

D2  8W  z 

assumption  under  this  restriction  is  that  — 2—  <<  1  (the  amount  of 

D2 

defocus  is  of  a  similar  order  as  that  of  spherical  aberration).  When 
this  is  the  case  the  intensity  distribution  becomes: 


I(r,,  z)  =  A2  (1  -  16  —  (W2Q  +  W40))  JQ2 


—  (\ J  4-  U  —  ^  T  2 
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(D2  -  d2)  +  kW, 


(IV-61) 


To  separate  the  two  aberrations  one  must  capitalize  on  the  basic 

difference  between  them.  Spherical  aberration  has  a  ^  dependence 

2 

while  defocus  has  only  a  r  dependence.  Varying  the  obscuration 

ratio  (that  is  varying  the  inner  diameter  of  the  aperture,  keeping  the 

outer  diameter  constant)  will  change  the  intensity  of  maxima  on  axis. 

Measure  two  maxima  on  axis.  These  will  occur  at  z  and  z  .  The 

ldl  2 

first  measurement  will  be  I,  with  obscuration  ratio  £  =  —  and  the 

1  D 

^  7 

second  will  be  1  with  e  *  — Then: 

2  2  D 
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(IV-62) 


A  second  measurement  will  replace  I  ,  z  and  £  by  I  ,  z  and  - 

i  i  122  : 

Then  subtracting  the  two  and  solving  for  W  yields: 


W  —  (e  2  -  £  2) 
40  3  {  l  C2  ' 


-S(‘[viRlH]) 
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^  =  32z  z  (e  2-e  2)I  (4V22  '  V  +  (I2Z1 

12  1  2  0 


I  z  )) 
1  2  ' 
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Putting  this  result  back  into  IV-62  yields: 


J3  \  _  2  _ _ 

V  3  32ZiZ2(Ei2-E22)10 


(4 1  (z 
0  2 
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(1  z 
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-  1,0X1  +  e  2) 


1  2 
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48ZlZ2I0(El2~e22 


)  ^ 


a  +  Er)(4I0 '  I2)  + 


z2(l  +  e22)(li  -  4Iq) 
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The  separation  of  the  aberrations  is  also  possible  by  considering  axial 
shifts  of  maxima,  with  different  obscuration  ratios.  It  is  clear  by 
looking  at  equation  IV-61  that  the  term  in  the  cosine  k(W^Q  (l  —  e  )  + 
W40(l-c4))  acts  like  a  phase  factor.  Unless  one  varies  e  the  two 
cannot  be  separated. 

Maxima  on  axis  will  occur  when 

87"  <1  -  e2)  +  kW20(l-E2)  +  kW40(l  -  c4)  =  (2n  +  1 ) tt 


_ P2(  1-c2) _ 

4X(2n+l)  -  8W2Q(l-e2)  -  8W40(l-e4) 


(IV-65) 


x. 

If  the  nLn  maximum  is  measured  at  z  with  obscuration  ratio  ■  and  at 

1  1 

z  with  obscuration  ratio  then  the  respective  axial  shifts  from  the 
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unaberrated  case  are; 


D2  ( 1-e  . 2 ,  ‘ 

Az  .  =  - 7 — 

i  4 


X(2n+1)  -  2W£0  (1— e±2)  -  2W4Q(1 


_  1 

_e_4)  ~  X(2n+1) 


i  -  1.  2 


D2(l-e  ,2) 
X 


«2o  +  u,c<1+si2> 


X(2n+l)|\(2n+l)  -  2«2 0 < l-e±2) 
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After  much  algebra  (see  Appendix  6)  the  resolved  aberration  coeffi¬ 
cients  are: 

2X2 (2n+l)2 _ 


W 


20 


(e12-c22) (1 -Ej2)  (l-e22) 

x  ^D2  |az1  (l-e22)2(l-e12)  -  Az2  ( 1-e^)2 (l-e22)J - 
4X(2n+l) (e14-e24)  ^ 

v  (D-^Cl-Cj2)  +  AZj4X  (2n+l) )  (D2  (l-e22)  +  Az24X(2n+l)) 


(IV-67) 
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2X2 (2n+l)2 


(g12-c22) (l-e12) (l-e22) 


D2^Az. (l-e22)2  -  Azo(l-e12)2J  +  4X(2n+l)AziAz2(EJ2-e22) 
(D2(1-ci2)  +  Azi4X(2n+l))(D2(l-e22)  +  Az24X(2n+l)) 
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Solving  the  above  expressions  is  a  tedious  task.  There  is  another 


drawback  to  separating  the  aberration  coefficients  using  axial  shifts 
—  the  difficulty  of  locating  the  ntn  order  maximum.  With  the  re¬ 
strictions  at  the  beginning  of  this  section  imposed,  small  values  of  z 
are  implied.  Small  z  means  high  orders  of  maxima  and  small  spacing 
between  maxima.  To  verify  the  order  of  the  maxima  in  both  cases  will 
require  assiduous  measurements. 


Astigmatism 


The  aberration  function  for  astigmatism  is: 
rl2 

W  =  4W?  -  cos20. 

D2  (IV-69) 


The  Rayleigh-Sommer^eld  formula  for  the  disturbance  function  with 


2  -2*.  D/2 


It/  ,  A  ikz  ik 
U(r2’  Z)  =  U  6  6 
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r  r 
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g  2  g  Z  2  i' 
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The  above  must  be  evaluated  numerically.  As  explained  earlier  the  ex¬ 
pected  diffraction  pattern  will  have  two  planes  of  symmetry.  Astigma¬ 
tism  is  the  predominant  aberration  generated  by  tilting  parallel 
plates.  Thus  an  empirical  study  of  this  aberration  is  possible. 


Coma 


The  aberration  functic i  for  coma  is: 

ri3 

W  =  8W,,  -  cos0, 

D3  ( IV— 71) 


The  Rayleigh-Sommerfeld  formula  with  coma  is: 


.  A  ikz  ik  „ 
U(r2 »  2)  ~  e  e  2z 


2  .  2*  ,D/2 

r22  c  r  ..  n 


LI. 


rlr2 


ik  — r—  -ik  -  cos  (0,-0,) 

e  2z  e  z  z  * 


0  d/2 
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The  above  has  no  apparent  analytic  solution  and  must  be  evaluated 
numerically.  The  amount  of  coma  produced  by  the  plane  parallel  plates 
is  small  particularly  in  comparison  to  the  amount  of  astigmatism.  It 
is  doubtful  any  comatic  effects  will  be  observable  during  the  experi¬ 
mental  phase. 


Summary 

Chapters  3  and  4  included  many  formulas  for  intensity  distribu¬ 
tions  on  and  off  axis.  Below  is  a  review  of  the  formulas  and  equation 


numbers . 


Unaberrated  Cases 


Uniform  intensity  plane  wave  incident  upon  circular  obscurration: 


On  axis: 


1(0,  0,  z)  =  A2 


Off  axis: 


I(r,,  z)  =  A2Jq2 


z2+d2 


/kdr:\ 

\~) 


(III-12) 


(III-26) 


Uniform  intensity  plane  wave  incident  upon  annular  aperture: 


On  axis: 


),  0,  z)  =  A2 


2  -  2cos  (  -  (D2  -  d2) 


CIIT-39) 


Off  axis: 


I(r2,  2)  =  !u(r2,  z)l2 


=  A2  J* 


\  2z  y  0  \  2z  /  J0  ^  2z  y  J0  ^  2z  J 


l_ 

x  COS^  CD2  -  d2)^ 
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Gaussian  intensity  plane  wave  incident  upon  circular  obscuration: 


On  axis: 


1(0,  z)  = 


1+  d2/Az “ 


(TV-73) 


I 


Off  axis: 


I(r2,  z)  =  A2  e~  2  JQ2 


m 


(III-34) 


Gaussian  intensity  plane  wave  incident  upon  annular  obscuration: 


On  axis: 


1(0,  z)  =  A2  I  e  2  * 


c  d2 

-e"  2  7 


-  2e  4 


L  D^Jcos(- 


(D2  -  d2) 
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Off  axis: 


l(r„,  z)  = 


[_  £  /  kDr„\  c  d2  o  /kdr9 \ 

•  •  v (— )  - 

j  t“  3 (5) '.  (¥)•-(!;«• 
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Uniform  plane  wave  incident  upon  displaced  obstacle  and  aperture: 


On  axis: 


1(0,  Zj,  z2)  A2  2  -  2cos  | 


8  21  z2 


(IV-75) 


Off  axis: 


I(v  V  Z2}  = 


125 


■■■m  -(its-?])' 


(III-74) 


Aberrated  Cases 

Uniform  intensity  plane  wave  containing  defocus  incident  upon 
circular  obscuration: 


On  axis: 


1(0,  z)  = 


8W2  02Y 


(IV-76) 


Off  axis: 


I(r2.  z) 


8Wz\2  0 

20  1 


(*) 


(IV-5) 


Uniform  intensity  plane  wave  containing  defocus  incident  upon 
annular  aperture: 


On  axis: 


1(0,  z)  = 


8W23z\2 


2cos  Is!  (d2 


-  d2)  i  + 


8W20z 


(IV-77) 
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Off  axis: 


(IV-U) 


Uniform  intensity  plane  wave  containing  spherical  aberration  inci¬ 
dent  upon  circular  obscuration 


On  axis: 


1(0, 


16 


w,  z 
<10 

d2 


(IV-78) 


Off  axis: 


l(r2  ,  z)  =  A2 


[■  -  ”  ¥] 
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Uniform  intensity  plane  wave  containing  spherical  aberration  inci 
dent  upon  annular  apertue 


On  axis: 
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16W  z  , 
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Off  axis: 


I(r  .  z)  =  A2  1  -  16V 
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(IV-60) 


From  all  the  preceding,  one  can  construct  a  general  formula  con¬ 
taining  all  of  the  soluble  elements  discussed  in  the  two  chapters: 


l(r,  ,  z  )  =  A2 
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Notes 


^  Max  Born  and  Emil  Wolf,  Principles  of  Optics  6th  ed.  (Oxford: 
Pergamon,  1930),  p.  470. 

Warren  J.  Smith,  Modern  Optical  Engineering  (New  York:  McGraw 
Hill,  1966),  p.  132. 
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H.  H.  Hopkins,  Wave  Theory  of  Aberrations  (Oxford:  Oxford 
University,  1950),  p.  141. 

^  James  E.  Harvey,  "An  Introduction  to  Aberration  Theory,"  Class 
handout,  University  of  New  Mexico,  Albuquerque. 

^  Joseph  W.  Goodman,  "Introduction  to  Fourier  Optics,"  (San 
Francisco:  McGraw-Hill,  1968),  pp.  123-124. 


Chapter  5 


EXPERIMENTAL  PROCEDURES 

The  experiment  contained  three  phases.  The  first  phase  involved 
calibrating  the  aberration  generator  an_  determining  the  gaussian  con¬ 
stants  for  the  beams  used.  The  second  phase  included  preliminary  runs 
which  suggested  the  ultimate  tests  for  the  final  (testing)  phase. 


Laboratory  Equipment 

Below  is  a  photograph  of  the  laboratory  set-up.  The  next  page  is 
a  diagram  labelling  the  components. 


13 

The  laser  used  was  a  Spectra-physics  helium-neon  laser  with  a  . 5mW 
power  output.  The  spatial  filter  was  composed  of  either  a  20X  or  60X 
microscope  objective  (to  change  the  gaussian  distribution  constant  for 
the  incident  beam)  and  a  25  micron  pinhole.  The  plane  parallel  plates 
were  half-inch  thick  and  made  of  BK-17  glass  with  index  of  refraction 
1.517.  The  plates  were  mounted  on  posts  allowing  for  somewhat  precise 
(+_  half  degree)  tilt  angles.  The  collimating  lens  was  a  100mm  focal 
length  achromatic  doublet.  The  lens  could  be  translated  to  introduce 
defocus.  The  obstacles  were  ball  bearings  of  diameters  3/16,  5/16, 

1/2  and  3/4  inches  mounted  on  plane  parallel  plates.  Disks  drilled 
with  different  sized  holes  provided  the  apertures.  A  Gaertner  travel¬ 
ling  microscope  with  a  20X  objective  was  used  to  make  measurements  of 
the  diffraction  structure  on  and  off  the  optical  axis.  Photographs  of 
the  spot  and  ring  patterns  were  made  with  a  Bausch  and  Lomb  stereo 
viewer  with  Polaroid  camera  attachment.  The  camera  attachment  had 
adjustable  shutter  speed  and  used  Polaroid  3000  black  and  white  film. 

Calibration 

The  aberration  coefficients  for  given  parallel  plate  tilt  angles 
in  tables  4.3  and  4.4  were  verified  by  taking  and  analyzing  interfero- 
grams  produced  by  aberrated  wavefronts.  A  Smartt  point  diffraction 
interferometer  created  the  interferograms.  The  interf erograms  in¬ 
cluded  a  variety  of  tilt  angles  using  one  plate,  two  plates  counter- 
rotated  or  two  plates  rotated  in  the  same  direction.  Photographs  were 
made  of  the  interference  pattern  using  Polaroid  3000  black  and  white 
film.  The  photographs  were  digitized  and  run  through  an  analysis 
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program  called  FRINGE  at  the  Fringe  Reduction  Facility  at  the  Air 
Force  Weapons  Laboratory  at  Kirtland  AFB.  The  results  were  then 
compared  with  theory  (see  Chapter  6). 

The  gaussian  intensity  widths  of  the  two  microscope  objectives 
used  in  the  spatial  filter  were  investigated  using  two  methods.  The 
first  method  was  the  knife  edge  test.  A  collecting  lens  focused  the 
plane  wave  output  of  the  lab  set-up  (no  plane  parallel  plates  or 
target)  into  a  power  meter  (a  hv  6328a  meter).  A  razor  blade  on  a 
translation  stage  is  introduced  into  the  beam  before  the  collecting 
lens. 


plane  waves  collecting 

from  system  lens 

(250mm  focal  length) 

Figure  5.2 

Set-up  for  Determination  of  Gaussian  Constants 

As  the  razor  blade  is  moved  farther  into  the  beam  the  percent  of  the 
beam  still  reaching  the  power  meter  i3  measured  as  a  function  of  the 
distance  of  translation  of  the  blade.  A  best  fit  line  for  the  data 


points  is  calculated.  The  distance  between  the  position  where  84.2^ 

transmission  occurs  and  the  point  where  15.85b  transmission  occurs  is 

designated  ix.  This  difference  corresponds  to  the  radius  of  the  beam 

at  the  e~~  points  (the  radius  at  which  the  intensity  drops  to  e~~ 

times  the  maximum  value).  The  desired  form  of  the  gaussian  distribu- 

rj-_ 

tion  is  given  by  e-c  n-.  Equating  the  two  and  substituting  \x  for  r 


yields : 


-2 

e 


(V—  1 


As  expected  the  gaussian  constant  will  depend  on  the  size  of  the  aper¬ 
ture  D  taken  as  the  normalization  value. 

The  second  method  of  finding  the  gaussian  constant  is  a  method  by 
Stijns.1  A  wire  or  cylindrical  object  is  used  in  place  of  the  blade. 
Instead  of  wholescale  screening  of  the  incident  beam,  selective  shad¬ 
ing  is  possible.  The  wire  is  moved  across  the  beam  until  a  minimum  in 
power  registered  by  the  meter  occurs.  As  the  figures  below  suggest 
this  minimum  corresponds  to  the  point  at  which  the  wire  is  centered  in 


minimum  power 
Figure  5-3 


Wire  Method  of  Gaussian  Constant  Tetemination 
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the  beam.  The  ratio  of  the  power  at  the  minimum  to  the  unobstructed 
power  is  compared  to  an  error  function.  Use  of  a  graph  enables  one  to 

obtain  a  value  for  the  ratio  of  the  wire  diameter  to  the  radius  WQ  of 

_2 

the  e  points.  In  a  similar  method  to  equation  V-1 ,  the  gaussian 
constant  can  be  determined. 

Preliminary  Runs 

The  dry  run  phase  involved  taking  measurements  to  verify  theory 
and  to  perfect  experimental  techniques  for  the  final  trials. 

Final  Tests 

The  ultimate  experiment  involved: 

1)  Verification  of  theory  for  the  unaberrated  case. 

2)  Measurement  of  defocus  aberration  coefficients  by  axial 
shifting  of  on-axis  extrema. 

3)  Measurement  of  spherical  aberration  coefficients  by 
shifting  of  on-axis  extrema. 

4)  Study  of  astigmatism. 

Verification  of  theory.  The  testing  of  the  theory  presented  in 
Chapter  3  was  accomplished  by  varying  a  number  of  parameters.  These 
parameters  were  the  obstacle  diameter,  the  aperture  diameter  (for  the 
annular  aperture),  the  gaussian  constant,  the  distance  along  the 
optical  axis  and  the  separation  of  the  obstacle  and  aperture.  The 
travelling  microscope  allowed  for  on  and  off  axis  measurements.  In 
addition,  some  photographs  were  taken  to  compare  with  computer  gener¬ 
ated  intensity  contour  plots. 
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Measurement  of  defocus  coefficients .  The  positions  of  maxima  and 
minima  on  axis  behind  various  annular  apertures  were  measured,  marked 
and  their  orders  determined.  Defocus  was  introduced  by  moving  the 
collimating  lens.  The  shift  of  the  lens  position  becomes  the  iz  term 
in  equation  IV-2.  The  coefficient  of  defocus  thus  obtained  can  be 
compared  to  the  coefficient  found  by  equation  IV-28  from  the  shift  on 
axis. 

It  is  important  to  note  that  by  decollimating  the  incident  beam 
the  wavefront  is  no  longer  planar.  The  beam  will  diverge  or  converge 
depending  on  the  sign  of  defocus.  But  the  spherical  obscuration  is 
mounted  on  a  plane  parallel  plate.  Will  the  now  diverging  or  converg¬ 
ing  wave  contain  more  aberrations  upon  plate  passage?  Since  the  plate 
is  not  tilted  some  spherical  aberration  will  be  produced  but  only  a 
negligible  amount.  The  non-converging  rays  still  strixe  the  plate  at 
essentially  perpendicular  angles. 

Measurement  of  spherical  aberration  coefficients .  Spherical 
aberration  was  introduced  by  inserting  one  or  two  parallel  plates  in 
the  diverging  beam  from  the  spatial  filter.  Measurements  were  taken 
using  the  same  method  employed  for  defocus.  The  aberration  coeffi¬ 
cients  were  calculated  from  the  expression  for  spherical  aberration  in 
tables  4.3  and  4.4.  These  were  compared  with  the  results  obtained 
from  equations  IV-58  and  IV-59* 

Study  of  astigmatism.  Astigmatism  was  introduced  into  the  wave- 
front  by  countertilting  two  plane  parallel  plates.  Photographs  were 
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taken  of  the  diffraction  pattern  of  a  circular  obstacle  for  one  degree 
increments  in  tilt  angle.  This  became  the  standard  by  which  every 
other  photograph  was  evaluated.  Various  axial  distances  and  obstacle 
diameters  were  tried.  To  investigate  the  annular  aperture,  the  obsta¬ 
cle  was  purposely  offset  from  center  so  the  diffraction  patterns  from 
the  outer  and  inner  edges  were  separated.  The  changes  in  each  pattern 
by  the  addition  of  astigmatism  were  photographed.  The  effect  of  the 
aberration  was  then  observed  for  a  variety  of  annular  apertures  and 
these  were  compared  with  the  circular  obscuration  cases.  Finally  it 
was  demonstrated  that  astigmatism  also  shifts  the  position  of  extrema 
on  axis. 

Numerical  Studies 

The  unaberrated  cases  and  those  with  defocus  and  spherical  aber¬ 
ration  lent  themselves  to  study  using  an  IBM  PC  computer.  Appendix  8 
contains  program  listings  of  the  routines  used  to  generate  the  plots 
and  graphs  in  this  thesis.  Appendix  7  gi  .'°s  a  comparison  of  the  use 
of  single  precision  numbers  versus  double  precision  numbers  for  the 


analysis. 


Notes 


^  Erik  Stijns,  "Measuring  the  Spot  Size  of  Gaussian  Beam  with  an 


Oscillating  Wire."  IIEE  Jour-.al  of  Quantum  Electronics.  QE-1 6  (Decern 
ber  1980),  1298-1299. 


Chapter  6 


EXPERIMENTAL  RESULTS 

The  presentation  of  the  results  of  the  experiment  will  follow  the 
order  presented  in  Chapter  5* 

Determination  of  the  Beam  Parameters 

Two  different  microscope  objectives  were  used  in  the  spatial 
filter.  The  6GX  objective  produced  a  wide  intensity  pattern,  suggest¬ 
ing  a  low  gaussian  constant.  This  objective  was  employed  to  study 
annular  apertures.  If  no  aperture  was  used  then  the  13/16  inch  diame¬ 
ter  mount  for  the  collimating  lens  became  the  aperture  stop  for  the 
system.  To  minimize  the  edge-diffracted  effects  from  the  mount,  the 
objective  used  was  by  the  20X  objective.  This  produces  a  narrow 
gaussian  beam. 

The  methods  employed  to  find  the  gaussian  constants  were  the 
knife-edge  test  and  the  wire  test. 

Knife-edge  Test 

The  power  of  the  uninterrupted  beam  reaching  the  detector  was 
0.413mW.  The  table  below  gives  the  knife  edge  position,  the  power 
measurement  and  the  percent  transmission  past  the  knife  edge. 


HO 


($  transmission)  =  400  -  6.184  x  (knife  edge  position)  (VI-1) 

The  coefficient  of  determination  is  .9883.  The  15.8$  and  84.2$  trans¬ 
mission  points  occur  at  62.13  and  51.06mm  respectively  giving  a  x  of 
1 1 .07 mm.  Then  by  equation  V-1,  the  gaussian  constant  is: 

C,n  =  — —  =  16320  D2sT2  (VI-2) 

For  example,  a  one  centimeter  diameter  aperture  will  yield  a  gaussian 
constant  of  1.632. 

The  20X  objective  yielded  a  total  beam  power  of  .438mW.  The 
table  of  data  for  this  objective  appears  below: 

Table  6.2 

Knife  Edge  Test  Data  for  20X  Objective 


Knife  Edge 
Position  (mm) 

Power  (mW) 

$  Transmission 

61 

•  431 

98.4 

62 

.421 

96. 1 

63 

.405 

92.5 

64 

.379 

86.5 

65 

•  341 

77.9 

65 

.300 

68.5 

67 

.249 

56.8 

68 

.184 

42.0 

69 

.136 

31.1 

70 

.090 

20.5 

71 

.048 

1 1  .0 
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The  least  squares  linear  fit  through  the  data  points  from  65mm  to  70mm 
gives  a  linear  function: 

(90  transmission)  =  847.9  -  11.8  x  (knife  edge  position)  (VI-3) 

The  coefficient  of  determination  is  .9967.  The  x  term  is  5.87mm,  the 

difference  between  the  15.8$  point  at  70.43mm  and  the  34.2$  point  at 

64.56mm.  The  gaussian  constant  for  the  20X  objective  is: 

C  =  58044  D2m“  2  (VI-4) 

2  OX 

Wire  Test 

Four  cylindrical  objects  of  different  diameters  were  used  in  this 
test.  These  are  listed  in  the  table  below: 

Table  6.3 

Diameters  of  Objects  Used  in  Wire  Test 

_ Measured  Diameter _ 

Object  Inches  Millimeters 

Paper  clip  wire  .0345  ^  .0005  .876  _+  .013 

Nail  #1  .0950  _+  .0005  2.413  +_  .013 

Nail  #2  .1178  +.0005  2.992  ^.013 

Pen  insert  .1246  +  .0005  3.165  +  .013 


The  measured  unobstructed  power  for  the  60X  objective  was  .41 3mW. 


c. 


—  =  43510  D2m“2 
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This  differs  from  the  result  found  by  the  knife  edge  test  by  a  factor 
of  over  2.5* 

The  comparison  between  the  tw  -  tests  was  much  closer  for  the  20 X 
objective.  In  this  test  the  backgr ;un :  ;ower  was  .44SmW.  The  table 
below  shows  the  parameters  to  find  *  . 


Wire  Test  Parame 

ters  for  2 OX 

objective 

Object 

Minimum 

Power  ( mV ) 

?a  /P 

a/W 

Vf  ( mm ) 

Paperclip 

wire 

.339 

.868 

.160 

5.48 

Nail  #1 

.300 

.670 

.435 

5.55 

Nail  #2 

.271 

.  525 

.515 

5.51 

Pen  insert 

.  266 

.594 

.540 

5.86 

The  average  gaussian  radius  for  the  e~-  points  is  5.68mm  (+  ,19mm). 
Notice  the  radius  measured  using  the  knife  edge  method  (5.87mm)  is 
within  one  standard  deviation  of  the  above.  The  gaussian  constant 
determined  by  the  wire  test  is: 

C  =  61990  D~m“"  (VI-7) 


The  wire  test  method  gave  very  consistent  results.  Unfortunate¬ 
ly,  the  results  were  much  too  high.  For  a  lomm  aperture  the  gaussian 
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constant  for  a  6  OX  microscope  objective  is  11.14  using  the  wire  method 
value.  This  constant  is  sufficiently  large  to  eliminate  the  effects 
of  the  outer  aperture.  The  diffraction  pattern  should  be  just  that  of 
an  obstacle  alone.  This  is  not  what  is  observed.  There  are  diffrac¬ 
tion  effects  from  the  outer  edge.  One  possibility  for  the  difference 
between  the  two  tests  is  the  beam  might  not  be  perfectly  gaussian.  In 
this  case  the  results  from  both  tests  are  suspect.  For  the  purpose  of 
plotting  graphs,  the  values  of  the  gaussian  constants  from  the  knife- 
edge  test  will  be  used  as  the  standards. 


Aberration  Coefficient  Curves 

i 

i 

Defocus 

The  coefficient  of  defocus  as  a  function  of  the  shift  in  position 
J  of  the  collimating  lens  is  given  by  equation  IV-2. 


8(F/f‘)2 


(IV-2) 


The  f  number  is  determined  by  the  radius  of  aperture  used.  The  aber¬ 
ration  coefficient  curves  for  aperture  diameters  of  7mm  (solid  line), 
lOnun  (dashed  line)  and  l6mn  (dotted  line)  appear  below. 


T 

S 

1 


[« 


Astigmatism 

Early  in  the  experiment  tie  laboratory  ret-uj  was  found  to  con¬ 
tain  background  astigmatism.  The  majority  :.\.s  astigmatism  seemed 
to  come  from  the  plane  parallel  plates.  t :.e  :  late  •  removed  the 

diffraction  patterns  of  ail  obstacles  were  t cytometric  ring  patterns 
predicted  by  theory.  When  the  plates  were  inserted  t.ne  patterns  were 
aberrated  to  various  degrees  depending  on  the  diameter  of  the  obsta¬ 
cle.  The  residual  astigmatism  seemed  to  be  positive.  The  astigmatism 
produced  by  countertilting  two  plates  has  a  negative  sign.  Thus  one 
can  balance  out  the  background  positive  astigmatism  by  introducing 
negative  aberration.  By  measuring  these  required  plate  tilt  angles 
for  various  sized  obstacles  the  amount  of  background  astigmatism  can 
be  determined.  The  table  below  lists  the  countertilt  angles  using 
obstacles  of  diameters  5/1 6,  1/2  and  5/8  inches  and  the  associated 
aberration  coefficients.  The  error  term  in  the  aberration  coefficient 


Table  6.7 

Measure  of  Background  Astigmatism 


Obstacle 
Diameter  (nun) 


Tilt  angle  required 

to  balance  out  W2? 

background  aberration  (°)  (wavelengths) 


7.938  5  1-5  .090  +  .016 

12.7  3-5  ^  .5  .113  ^  -032 

2.5  +  »5  .09 0  +  .G3o 


15.875 
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-  .1 


is  the  average  deviation  based  upon  the  half  angle  uncertainty  in  the 
tilt  angle.  The  aberration  coefficient  for  the  background  astigmatism 
is  taken  as  .0975  ,  the  average  of  the  above  values.  This  is  added  to 
the  expression  from  Table  4.4  for  W  . 


W 


22 


.0975X  - 


(VI-8) 


Interferograms  using  the  point  diffraction  interferometer  were 
made  to  verify  the  above.  The  data  points  obtained  from  interferogran 
analysis  are  plotted  against  the  curve  of  equation  VI-3  for  a  13/16 
inch  diameter  aperture.  The  error  bars  reflect  the  uncertainty  in 
tilt  angles  and  the  associated  uncertainty  in  W  . 


-  • 
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AST  I CMfiT I  ABERRATION  COEFFICIENTS  FOR  13/16  INCH  flRERTL'RE 


Graph  6.2 


Seven  interferograms  were  selected  from  over  thirty  taken  to  be 


reduced.  These  are  shown  below.  The  first  is  of  one  plate  inserted 
with  no  tilt,  ^he  next  three  show  two  countertilted  plates  of  4.5°> 
10°  and  20°  each.  The  effect  of  increasing  astigmatism  is  seer,  to 
compress  the  interferogram  pattern  from  concentric  circles  to  concen¬ 
tric  ellipses.  The  last  three  show  the  effects  of  two  plates  lilted 
in  the  same  direction  at  angles  of  5°»  10°  and  15°.  Ccma  enters  into 
these  patterns  and  sweeps  the  fringes  to  one  side.  Photograph  6.6 
contains  a  dark  spot  in  the  center  of  the  pattern.  This  is  the  result 
of  the  pnotograph  being  taken  of  the  pattern  on  the  side  of  negative 
defocus.  All  the  other  photographs  show  patterns  on  the  positive 
defocus  side  of  the  null  (zero  focus)  position. 


Photograph  6.1 


One  plate  -  no  tilt 


Photograph  6.2 

Two  plates  -  4-5°  countertilt 


Photograph  6.3 

Two  plates  -  10°  countertilt 


c 
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In  each  photograph  the  center  of  the  pattern  was  measured,  (in 
the  case  of  Photograph  6.7  which  is  not  symmetric,  the  position  of 
center  had  to  be  estimated.)  This  became  the  center  of  the  pattern  to 
be  digitized  for  the  computer  analysis.  This  was  necessary  because  a 
shift  in  position  from  the  center  will  introduce  extraneous  aberra¬ 
tions.  To  see  this  take  a  simple  example.  Assume  the  only  aberration 
present  is  defocus.  Let  the  x-position  of  the  center  be  shifted  by 
Ax,  (”x  <<  1).  Then  the  aberration  function  becomes: 

W  =  W  r2  =  V  (x2  +  y2)  =}  W  ( (x  +  Ax)2  +  v2) 

20  20  7  20 

=  W  (x2  +  2xAx  +  0  4-  y2) 

20 

=  W  (x2  +  y2)  +  2AxW  x 
20  20 

=  W  r2  +  2AxW  r  cos0 
20  20 

=  V  r2  +  W  'r  cos0  (VI-9) 

20  11 

A  shift  in  the  x-direction  manifests  itself  by  the  introduction  of  an 
extraneous  tilt  term.  If  the  pattern  is  not  centered  when  put  into 
the  reduction  program  the  output  will  contain  many  of  these  manufac¬ 
tured  aberrations. 


Verification  of  Theory  for  the  Unaberrated  Case 


Three  features  of  the  diffraction  pattern  predicted  by  the  theory 
of  Chapter  3  were  investigated.  The  first  was  the  position  and  spacing 


of  the  rings  in  the  diffraction  pattern.  This  was  accomplished  in 
conjunction  with  the  observations  of  the  other  two  features.  The 
second  was  the  increase  of  the  overall  pattern  size  for  increased 
optical  axis  distance  Z.  The  third  was  the  decrease  in  the  overall 
pattern  size  for  increased  obstacle  diameter. 

Circular  Obstacle 

The  following  photographs  show  the  diffraction  patterns  behind  a 
3/16  inch  diameter  circular  obstacle  observed  at  20,  25,  50  and  75 
centimeters.  The  20X  microscope  objective  was  used  in  the  spatial 
filter.  The  stereo-viewer  was  set  at  a  magnification  of  6X.  The 
exposure  times  were  all  1/2  second.  The  photographs  demonstrate  not 
only  the  increase  in  pattern  size  with  increased  z,  but  the  regular 
spacing  of  the  ring  structure  and  the  brightness  of  the  central  spot 
in  comparison  to  the  rings. 

Quantitative  measurements  of  the  increased  pattern  size  came  from 
looking  behind  a  5/16  inch  diameter  ball  bearing  at  40cm  and  at  75cm. 
The  tables  below  compare  experiment  and  theory  for  these  cases.  The 
calibrations  on  the  travelling  microscope  allowed  for  measurement 
accuracy  within  +_  .005mm.  Additional  uncertainty  arose  from  the 
judgement  of  where  the  brightest  and  darkest  points  in  the  ring  struc¬ 


ture  occurred. 
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Table  6.8 


Comparison  of  Theory  and  Experiment: 
Positions  of  Ring  Minima  and  Maxima. 
5/16  Inch  Diameter  Obstacle  Viewed  at  40cra 


Distance  from 

Center  (mm) 

Ring  Structure 

Experiment 

Theory 

Percent 

Difference 

Minimum  #1 

.023 

.024 

4.2 

Maximum  #1 

.046 

.039 

17.9 

Minimum  #2 

.057 

.056 

1.8 

Maximum  #2 

.072 

.071 

1.4 

Minimum  #3 

.091 

.088 

3.4 

Maximum  #3 

.108 

.103 

4.9 

Minimum  #4 

.126 

.120 

5.0 

Maximum  ft 4 

.143 

.135 

5.9 

Minimum  #5 

.154 

.152 

1.3 

Maximum  #5 

.172 

.167 

6.0 

Minimum  #6 

.194 

.183 

6.0 

Maximum  #6 

.209 

•  199 

5.0 

Minimum  #7 


227 


215 


5.6 


157 


Table  6.9 


Comparison  of  Theory  and  Experiment: 
Positions  of  Ring  Minima  and  Maxima. 
5/16  Inch  Diameter  Obstacle  Viewed  at  75cm 


Distance  from 

Center  vam) 

Ring  Structure 

Experiment 

Theory 

Percent 

Difference 

Minimum  #1 

.045 

.046 

2.2 

Maximum  ft  1 

.070 

.073 

4.1 

Minimum  #2 

.103 

.105 

1.9 

Maximum  #2 

.134 

.134 

0.0 

Minimum  #5 

.162 

.164 

1.2 

Maximum  #5 

.169 

.194 

2.6 

Minimum  ffA 

.230 

.224 

2.7 

Maximum  rfA 

.256 

.254 

0.8 

Minimum  #5 

.287 

.284 

1.1 

Maximum  #5 

.317 

•  313 

1.3 

Minimum  #6 

.344 

.344 

0.0 

Maximum  #6 

•  374 

.373 

0.3 

Minimum  #7 

.405 

.404 

0.2 

The  above  tables  show  that  the  experiment  agrees  with  predicted 
values.  The  spacing  between  the  rings  shows  almost  a  100$  increase 
from  the  40cm  measurements  to  the  respective  75cm  measurements.  There 


were  larger  deviations  from  theory  for  the  pattern  viewed  at  40cm  than 
for  the  75cm  pattern.  This  is  reasonable  because  the  uncertainty  in 
measurements  results  in  a  larger  percentage  of  difference. 

The  photographs  on  the  next  page  show  the  diffraction  patterns 
behind  a  3/ 16  inch  diameter  obstacle  and  a  5/16  inch  diameter  obstacle 
viewed  at  a  distance  of  60cm.  A  20X  microscope  objective  was  used  and 
the  binocular  viewer  set  to  magnification  6X.  Photograph  6.12  had  an 
exposure  time  of  1/8  second  while  photograph  6.13  was  timed  for  1/2 
second.  The  diffraction  pattern  for  the  3/16  inch  obstacle  is  much 
larger  than  that  for  the  5/16  inch  obstacle. 


Photograph  6.12 

3/16"  diameter  obstacle  at  60cm 


Photograph  6. 1 3 

5/16"  diameter  obstacle  at  60cm 


The  positions  of  the  first  six  ring  intensity  maxima  were  taken 
for  the  above  cases  as  well  as  for  a  1/2  inch  diameter  obstacle 
comparisons  with  theory  appear  in  the  tables  below. 
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Table  6.10 

Comparison  of  Theory  and  Experiment: 
Positions  of  Ring  Maxima. 

3/16  Inch  Diameter  Obstacle  Viewed  at  60cm 


Ring  # 
(Maximum) 


Distance  from  Center  (mm) 
Experiment  Theory 


Percent 

Difference 


1 

.095 

.097 

2. 1 

2 

.180 

.178 

1  . 1 

3 

.262 

.258 

1.6 

4 

.345 

•  338 

2. 1 

5 

.415 

.418 

0.7 

6 

.511 

.493 

2.6 

Table  6.11 

Comparison  of  Theory  and  Experiment: 
Positions  of  Ring  Maxima. 

5/16  Inch  Diameter  Obstacle  Viewed  at  60cm 


Ring  # 
(Maximum) 


Distance  from  e'en  ter  (mm) 

Percent 

Experiment  Theory  Difference 


1 

.057 

.056 

1.7 

2 

.107 

.107 

0.0 

T 

J 

.153 

.155 

1.3 

4 

.203 

.203 

0.0 

5 

.252 

.251 

0.4 

6 

.  29Q 

.299 

0.0 

Table  6.12 
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Comparison  of  Theory  and  Experiment: 
Positions  of  King  Maxima. 

1/2  Inch  Diameter  Obstacle  Viewed  at  bOcm 


Ring  # 

Percent 

(Maximum) 

Experiment 

Theory 

Difference 

1 

.034 

.036 

5.6 

2 

.069 

.067 

3.0 

3 

.100 

.127 

3-1 

4 

.123 

.127 

3-1 

5 

.154 

.157 

1.9 

As  the  diam-  -  of  the  obstacle  increases  the  scale  of  the  pat¬ 
tern  diminishes.  As  previously  seer,  the  smallest  pattern  has  the 
largest  average  percent  difference  from  theory. 

Annular  Aperture 

To  test  the  theory  for  the  annular  aperture,  three  obstacle  sizes 
and  six  aperture  sizes  were  used.  The  table  below  lists  the  obscura¬ 
tion  ratios  •.  for  the  possible  combinations. 
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Table  6.13 
Obscuration  Ratios 


Obstacle  Diameters  (cm) 


Aperture 
Diameters  (cm) 

3/16"  (.47625) 

5/16"  (.79375) 

1/2"  (1.27) 

•  / 

.  66 

3/6"  (.9523) 

.5 

.833 

_ 

1.0 

.  47o 

.794 

_ 

1/2"  (1.27) 

.375 

.  625 

1  .CO 

3/o"  (1.5675) 

.3 

.5 

•8 

1 . 6 

.193 

.  496 

•  .  37*4 

The  following  photographs  ana  graphs  compare  the  experimental 
results  with  theoretical  predictions  for  a  variety  of  combinations. 

The  scale  factors  in  the  intensity  contour  plots  were  reduced  from 
those  used  in  Chapters  3  and  4  in  consideration  of  the  reduced  overall 
intensity  resulting  from  a  gaussian  intensity  distribution  of  the 
incident  wave.  They  were  reduced  again  by  half  in  the  contour  plots 
accompanying  the  last  two  photographs. 
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Photograph  6.14 


3/16"  diameter  obstacle  with  7mm  diameter  apertur 
Viewing  distance  =  62cm.  Obstacle-aperture  separation 
Magnification  5X.  Exposure  time  =  1/2  sec. 


Photograph  6.15 


3/16"  diameter  obstacle  with  10mm  diameter  aperture. 
Viewing  distance  =  59- 1cm.  Obstacle-aperture  separation  =  1.3c 
Magnification  6X.  Exposure  time  =  1  sec. 


Photograph  6.13 


1/2"  diameter  obstacle  with  5/8"  diameter  aperture. 
Viewing  distance  =  2.120m.  Obstacle-aperture  centered 
Magnification  6X.  Exposure  time  =  1  sec. 


16  M, 


RADIAL  DISTANCE  (.601  N) 

I.D.  =  .0127  H.  Z  DISI.  =  2.120  H.  GAUSS .  CONST 


Graph  6.12 


Intensity  profile  of  1/2"  diameter  obstacle 
with  5/6"  diameter  aperture 


Photograph  6.19 


1/2"  diameter  obstacle  with  5/6"  diameter  aperture. 
Viewing  distance  =  2.005m.  Obstacle-aperture  centered 
Magnification  6X.  Exposure  time  =  1  sec. 
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The  photographs  show  excellent  agreement  with  theory.  Photograph 

6.15  depicts  a  3/1 6  inch  diameter  obstacle  with  a  10mm  diameter  aper¬ 
ture  at  an  on-axis  minimum.  The  corresponding  -'-itensity  contour  plot 
does  not  have  a  dark  center.  This  is  due  to  the  resolution  of  the 
contour  plotting  routine.  The  central  minimum  in  this  -lot  is  too 
narrow  for  the  routine  to  observe.  It  does  she*  up  in  the  companion 
intensity  distribution  curve.  Of  interest  is  the  fact  that  the  on- 
axis  intensities  are  non-sero  in  the  minimum  on-axis  plots.  This  is  a 
consequence  of  the  gaussian  nature  of  the  incident  beam.  This  effect 
is  most  easily  seen  in  photograph  6.19  where  the  center  of  the  dif¬ 
fraction  pattern  is  not  dark.  Photographs  6.16  and  6.1”  show  very 
similar  patterns.  This  is  expected  because  of  similar  obscuration 
ratios.  Photograph  6.15  corresponds  to  a  ratio  of  .43  and  photograph 

6.16  to  a  ratio  of  .5. 

A  quantitative  comparison  of  radial  positions  of  diffraction 
pattern  minima  and  maxima  for  theory  and  experiment  follows. 


Table  6.14 

Comparison  of  Theory  and  Experiment: 
Positions  of  Ring  Maxima  and  Minima. 

3/16  Inch  Diameter  Obstacle  with 
1cm  Diameter  Aperture,  Viewing  Distance  =  50.2cm 
Obstacle-Aperture  Separation  =  1cm 


Qualitative 

Description 


Normalized  _ Cente 

Intensity 

'Theory)  Experiment 


Distance  from 
Center  (mm) 


Theory 


Percent 

Difference 


Center  (minimim) 

.061 

.000 

.000 

Maximum  ring  1 

.475 

.030 

.031 

Minimum 

.000 

.059 

.055 

Maximum  ring  2 

.277 

.073 

.076 

Dark  band 

~  .0 

.092-. 154 

. : 13-.130 

Maximum  ring  3 

.132 

.171 

.165 

Minimum 

.000 

.190 

.188 

Maximum  ring  4 
(faint) 

.071 

.208 

.209 

Dark  band 

~  .0 

.21 6-. 282 

.242-. 271 

Maximum  ring  5 

.086 

.296 

.298 

Minimum 

.000 

.317 

.323 

Maximum  ring  6 

( faint) 

.023 

.332 

•  341 

Table  6.15 


Comparison  of  Theory  and  Experiment: 

Positions  of  Ring  Maxima  and  Minima. 

3/16  Inch  Diameter  Obstacle  with  1.6cm  Diameter  Aperture, 
Viewing  Distance  =  59*9cm,  Obstacle-Aperture  Separation  =  .8cm 


Normalized 

Distance  from 

Center  (mm) 

Qualitative 

Description 

Intensity 

(Theory) 

Experiment 

Theory 

Percent 

Difference 

Center  (mininim) 

•  3H 

.000 

.000 

-- 

Maximum  ring  1 

.697 

.027 

.023 

17.4 

Minimum 

.000 

.051 

.057 

10.5 

Maximum  ring  2 

.196 

.102 

.099 

3.0 

Minimum 

.000 

.159 

.147 

8.2 

Maximum  ring  3 

.103 

.172 

.174 

1 . 1 

Minimum 

.000 

.227 

.226 

0.4 

Maximum  ring  4 

.056 

.245 

.248 

1.2 

Minimum 

.000 

.304 

.303 

0.3 

Maximum  ring  5 

.026 

.322 

.323 

0.3 

Minimum 


000 


3TS 


380 


0.5 
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Table  6.16 

Comparison  of  Theory  and  Experiment: 

Positions  of  Ring  Maxima  and  Minima. 

5/16  Inch  Diameter  Obstacle  with  5/8  Inch  Diameter  Aperture, 
Viewing  Distance  =  62.8cm 


Distance 

from 

No rmalized 

Center  ( 

mm) 

Qualitative 

Intensity 

Percent 

Description 

(Theory) 

Experiment 

Theory 

Difference 

Center  (rainimim) 

1 . 280 

.000 

.000 

_ 

Minimum 

.  000 

.035 

.034 

2.9 

Maximum  ring  1 

.079 

.073 

.071 

2.8 

Minimum 

.  0.  0 

.  090 

.090 

0.0 

Maximum  ring  2 

'  -  4 

•  • 

.111 

.109 

1.8 

Minimum 

.... 

.1?' 

.132 

0.3 

Maximum  ring  ? 

•  0'5? 

.174 

.172 

1.2 

Minimum 

.000 

.1  9~ 

•  1  ^2 

0.  0 

Maximum  ring  4 

.045 

•  c.  ^70 

.  209 

0.5 

Maximum  ring  5 

.020 

•  274 

.272 

0.7 

Minimum 

.000 

.291 

.  290 

0.3 

Maximum  i  ■  ~  6 

.032 

.310 

.?■.- 

0.3 

Maximu... 

.015 

.371 

7  *“• 

•  >  1  C 

2.3 

Minimum 

.000 

•  s  -  -■ 

0.3 

Maximum  ring  8 


024 


0 


The  above  tables  demonstrate  close  agreement  between  experiment  and 
theory.  An  interesting  feature  is  the  relative  nature  of  minima  and 
maxima.  Table  6.15  shows  the  central  minimum  has  an  intensity  (theo¬ 
retical)  higher  than  the  maximum  intensity  in  most  of  the  surrounding 
rings!  Clearly  the  appearance  of  a  maximum  or  minimum  depends  or.  the 
immediate  intensity  background. 

The  final  study  for  the  unaberrated  case  involved  the  positions 
and  spacing  of  maxima  and  minima  on  axis  behind  annular  apertures. 
Three  combinations  of  apertures  and  obstacles  were  employed.  The 
first  was  a  3/16  inch  diameter  obstacle  with  a  7mm  diameter  aperture. 

A  maximum  or  minimum  on  axis  was  found  and  its  postion  measured  to  the 
nearest  centimeter.  The  corresponding  position  on  the  scale  of  the 
travelling  microscope  was  noted.  The  microscope  was  translated  along 
the  optical  axis  and  the  spacings  between  the  extrema  measured.  (The 
microscope  could  be  moved  a  total  distance  of  5.5cm.)  The  scale  al¬ 
lowed  for  a  maximum  uncertainty  in  measurement  of  approximately  .01cm. 
The  first  case  includes  measurements  starting  at  various  positions 
along  the  optical  axis.  The  last  two  involve  just  one  starting  posi¬ 
tion.  A  quirk  of  the  scale  on  the  travelling  microscope  is  that  scale 
measurements  are  inversely  related  to  optical  axis  distances.  A 


decrease  of  position  by  2cm  on  the  microscope  scale  means  an  increase 
of  2cm  along  the  optical  axis. 


ISO 


Table  6.17 

Comparison  of  Theory  and  Experiment: 

Positions  and  Spacing  of  On  Axis  Extrema. 

16  Inch  Diameter  Obstacle  with  7mm  Aperture  Separated  by  1cm 


Measured  Starting  Position  (  +  «5cm):  33cm 


Experiment _  _ Theory 


Extremum 

Type 

Scale 

Reading 

(cm) 

Distance  from 
Starting  point 
(cm) 

Optical  axis 
Distance 
(cm) 

Distance  from 
Starting  Point 
(cm) 

% 

Difference 

Minimum 

4.85 

0.0 

33.31 

0.0 

_ 

Maximum 

3-94 

.91 

34.36 

1.05 

13.3 

Minimum 

2.68 

2.17 

35.48 

2.17 

0.0 

Maximum 

1.60 

3.25 

36.67 

3-36 

3-3 

Minimum 

.33 

4.52 

37.96 

4.65 

2.8 

Measured  Starting  Position  (i  .5cm):  40cm 


Experiment _  _ Theory 


Extremum 

Type 

Scale 

Reading 

(cm) 

Distance  from 
Starting  point 
(cm) 

Optical  axis 
Distance 
(cm) 

Distance  from 
Starting  Point 
(cm) 

% 

Difference 

Maximum 

5-22 

0.0 

39.33 

0.0 

__ 

Minimum 

3.76 

1.46 

40.82 

1.49 

2.0 

Maximum 

2.07 

3-17 

42.42 

3.09 

2.6 

Minimum 

.37 

4.87 

44.15 

4.82 

1.0 

Minimum 


-  -  /•.  r. 


pUT\"  '  "  ■" 


k 


U 
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Table  6.17  (Cont.) 

Comparison  of  Theory  and  Experiment: 

Positions  and  Spacing  of  On  Axis  Extrema. 

3/1 6  Inch  Diameter  Obstacle  with  7mm  Aperture  Separated  by  1cm 


Measured  Starting  Position  (  +  .5cm):  50cm 


Experiment 


Theory 


Extremum 

Type 

Scale 

Reading 

(cm) 

Distance  from 
Starting  point 
(cm) 

Optical  axis 
Distance 
(cm) 

Distance  frcm 
Starting  Point 
(cm) 

Difference 

Maximum 

5.85 

0.0 

50.34 

0.0 

tmmm 

Minimum 

3.20 

2.65 

52.82 

2.48 

6.6 

Maximum 

.64 

5-21 

55.56 

5.22 

0.2 

*4 


Positions  and  Spacing  of  On  Axis  Extrema. 

3/16  Inch  Diameter  Obstacle  with  5/8  Inch  Diameter  Aperture 
Separated  by  .25cm.  Starting  Position  =  Minimum  at  76cm  (+  .5 


Experiment _  _ Theory 


Extremum 

Type 

Scale 

Reading 

(cm) 

Distance  from 
Starting  point 
(cm) 

Optical  axis 
Distance 
(cm) 

Distance  from 
Starting  Point 
(cm) 

at 

to 

Difference 

Minimum 

5-34 

0.0 

75.53 

0.0 

_ 

Maximum 

4.71 

.63 

76.16 

.63 

0.0 

Minimum 

4.09 

1.25 

76.81 

1.28 

2.3 

Maximum 

3-58 

1.96 

77.46 

1.93 

1.6 

Minimum 

2.74 

2.60 

78.13 

2.60 

0.0 

Maximum 

2.03 

3-31 

78.81 

3.28 

0.9 

Minimum 

1.35 

3-99 

79.50 

3.97 

0.5 

Maximum 

.57 

4.77 

80.20 

4.67 

2.1 
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Table  6.19 

Comparison  of  Theory  and  Experiment: 

Positions  and  Spacing  of  On  Axis  Extrema. 

5/16  Inch  Diameter  Obstacle  Centered  with  5/8  Inch  Diameter  Aperture. 
Starting  Position  =  Minimum  at  60cm  (+  .5cm) 


Extremum 

.Type 

Experiment 

Theory 

% 

Difference 

Scale 

Reading 

(cm) 

Distance  from 
Starting  point 
(cm) 

Optical  axis 
Distance 
(cm) 

Distance  from 
Starting  Point 
(cm) 

Minimum 

5-43 

0.0 

60.22 

0.0 

— 

Maximum 

4.94 

.49 

60.71 

.49 

0.0 

Minimum 

4.43 

1.00 

61.21 

•  99 

1.0 

Maximum 

3.94 

1.49 

61.71 

1.49 

O.C 

Minimum 

3.41 

2.02 

62.23 

2.01 

0.5 

Maximum 

2.87 

2.56 

62.75 

2.53 

1.2 

Minimum 

2.35 

3.08 

63.28 

3.06 

0.7 

Maximum 

1.84 

3-59 

63.82 

3.60 

0.3 

Minimum 

1.23 

4.20 

64.37 

4.15 

1.2 

Maximum 

.73 

4.70 

64.93 

4.71 

0.2 

Minimum 

.14 

5.29 

65.50 

5.28 

0.2 

The  above  tables  show  excellent  agreement  between  theory  and 
experiment.  Table  6.17  demonstrates  the  increased  spacing  between 
extrema  with  increased  distance  along  the  optical  axis.  This  table 
had  the  highest  average  percent  difference  as  well  as  the  largest 
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single  deviation  from  theory  (almost  30%  difference  between  theory  and 
experiment  for  the  spacing  of  extrema  at  75cm).  This  table  illus¬ 
trates  the  major  difficulty  in  measuring  the  positions  of  on  axis 
extrema.  When  the  spacing  between  adjacent  maxima  and  minima  is  large 
(more  than  2cm),  it  becomes  difficult  to  judge  the  exact  position  of 
the  extremum.  The  pattern  changes  so  slowly  with  z  translation  large 
amounts  of  uncertainty  in  measurements  are  introduced.  This  source  of 
potential  error  can  be  avoided  by  limiting  measurements  to  the  region 
on  the  optical  axis  where  extremum  spacing  is  small.  The  extent  of 
this  region  increases  with  larger  aperture  size  and  smaller  obscura¬ 
tion  ratio. 


Verification  of  Theory  for  the  Aberrated  Case 


Defocus 

Defocus  was  introduced  by  axially  translating  the  collimating 
lens.  For  each  shift  of  the  lens  the  change  in  the  on-axis  position 
of  a  minimum  of  known  order  was  measured.  This  was  possible  by  shift¬ 
ing  the  lens  in  . 1  mm  movements  so  as  to  make  the  resulting  on-axis 
change  be  less  than  the  spacing  between  adjacent  minima.  The  aberra¬ 
tion  coefficient  W  for  each  measurement  was  calculated  from  equation 
IV-28.  These  points  were  plotted  against  the  aberration  coefficient 
curve  given  by  equation  IV-2.  There  was  an  uncertainty  of  .05cm  in 
the  measurements  of  the  on-axis  extrema.  The  error  bars  based  on  this 


uncertainty  were  included  in  the  plots,  but  the  error  range  for  any 


given  point  is  practically  negligible. 

The  following  graphs  illustrate  the  comparison  of  theory  and 
experiment  for  three  cases:  The  first  was  a  3/16  inch  diameter  obsta 
cle  with  a  7mm  diameter  aperture.  The  measurements  were  based  on  che 
7^  order  minimum  at  75.11cm  behind  the  aperture.  The  second  was  a 
3/16  inch  diameter  obstacle  with  a  10mm  diameter  aperture.  The  mea¬ 
surements  were  based  on  the  20  ^  order  minimum  at  76.65cm.  The  final 
case  involved  a  3/16  inch  diameter  obstacle  and  a  16mm  diameter  aper- 
ture.  The  60Ln  order  minimum  at  76.69cm  was  used. 


DEFOCUS:  3/16  INCH  OESIftCLE  WITH  7  KH  ft PEP.IU RE 
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The  graphs  show  remarkable  agreement  between  theory  and  experi¬ 
ment.  The  shift  of  extrema  on  axis  is  a  very  effective  method  of 
measuring  defocus. 


Spherical  Aberration 

The  attempts  to  verify  the  theory  for  spherical  aberration  were 
mostly  unsuccessful.  Because  the  coefficient  of  spherical  aberration 
is  proportional  to  the  fourth  power  of  the  diameter  of  the  outer  edge 
of  the  annular  aperture  the  largest  sized  apertures  were  used.  A  16mm 
diameter  aperture  had  equal  to  -.0766  wavelengths  when  two  half 
inch  thick  plates  were  inserted  in  the  system.  From  equation  IV-58, 
the  shift  on  axis  for  the  ntJl  order  maximum  with  spherical  aberration 
present  i3: 


Az 


(2n+l) X 


(VI-10) 


A  1/2  inch  diameter  obstacle  with  a  16mm  diameter  aperture  has 
its  9^  order  maximum  at  1.97m.  The  addition  of  two  plates  will 
result  in  a  shift  of  position  of  nearly  1cm.  Even  a  shift  as  large  as 
this  is  hard  to  distinguish  when  one  considers  the  distance  to  the 
nearest  minimum  is  8cm.  To  pinpoint  the  position  of  the  maximum  is 
very  difficult  without  intensity  measurements.  The  uncertainty  is  at 
least  on  the  order  of  a  centimeter.  In  short,  the  best  case  gives  a 
predicted  shift  less  than  the  experimental  error. 

Many  of  the  results  were  conflicting,  e.g.,  the  insertion  of  one 
plate  seemed  to  move  the  1 1 order  maximum  at  1.7m  a  distance  of 


I  0.4cm.  This  gave  a  value  for  W  of  . 043 \  —  a  1  J>%  error.  When 

uu 

remeasured,  the  shifts  for  one  and  two  plates  appeared  to  be  1.18cm 

and  3.1cm,  giving  values  for  W  of  . 1 27 \  and  .  33S\  —  4  times  larger 

-  0 

than  actual. 

Additional  uncertainty  arose  from  the  need  to  reposition  the 
collimating  lens  when  a  plate  was  inserted  or  removed.  For  each  plate 
the  lens  had  to  be  moved: ^ 


n 


1.317-1 

1.517 


1 . 27cm 


=  .433cn  (VI-11) 

This  displacement  was  verified  by  observing  the  size  of  the  beam  at 
various  distances  along  the  optical  axis,  and  by  observing  the  inter¬ 
ference  pa;,:-,  m  produced  by  a  shear  plate  (when  the  beam  is  colli¬ 
mated,  one  obtains  a  null  pattern).  As  demonstrated  in  the  previous 
section  on  defocus,  a  change  in  the  position  of  the  collimating  lens 
introduces  defocus  which  also  shifts  the  on-axis  position  of  minima. 
For  a  16mm  diameter  aperture,  an  error  of  .1mm  in  the  position  of  the 
collimating  lens  from  the  proper  position  will  yield  a  defocus  term  on 
the  order  of  . U  which  will  shift  the  positions  of  extrema  at  1.7m  by 
nearly  2cm. 

One  can  introduce  a  substantial  amount  of  spherical  aberration  by 


reversing  the  achromatic  collimating  lens.  Unfortunately,  this  involves 
realigning  the  system.  This  technique  was  considered  impractical. 
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Study  of  Astigmatism 


Circular  Aperture 

To  study  the  effect  of  astigmatism  on  the  diffraction  pattern  of 
a  circular  obstacle,  a  standard  had  to  be  established.  The  5/16  inch 
diameter  obstacle  provided  this  standard.  There  were  two  benefits  to 
using  this  obstacle.  Because  this  diameter  fell  in  the  center  of  the 
range  of  obstacle  diameters  used  (3/16  inch,  5/16  inch  and  1/2  inch), 
one  could  investigate  the  change  in  pattern  by  increasing  or  decreas¬ 
ing  the  diameter.  In  addition,  the  diffraction  pattern  of  the  5/16 
inch  diameter  obstacle  changed  at  a  convenient  rate  with  the  increase 
of  parallel  plate  tilt  angle.  An  increase  of  one  degree  in  tilt  angle 
resulted  in  a  noticeable,  but  not  dramatic,  change.  A  series  of  photo¬ 
graphs  of  degree  increments  in  tilt  angle  offered  a  wide  spectrum  of 
astigmatic  effects  and  yet  allowed  one  to  observe  subtle  differences. 

The  following  photographs  are  of  a  5/16  inch  diameter  obstacle 
observed  at  60cm.  The  first  photograph  used  magnification  of  6X  and 
an  exposure  time  of  1/2  seconds.  All  the  other  photographs  were 
magnified  seven  times  and  had  one  second  exposure  times.  The  astigma¬ 
tism  was  produced  by  two  half-inch  plane  parallel  plates.  In  the 
range  from  tilt  angles  of  10°  to  18°  photographs  were  made  of  plates 
countertilted  and  tilted  in  the  same  direction.  This  was  to  observe 
whether  coma  had  any  effect  on  the  diffraction  pattern.  Each  caption 
lists  the  amount  of  astigmatism  present  (given  by  equation  VI-3)  as 
well  as  the  upper  and  lower  bounds  of  the  aberration  (in  parentheses) 
based  upon  an  uncertainty  of  +_  .5°  in  the  tilt  angle. 


Photograph  6.26 

13°  countertilt.  W  =  -.509\  (-.463X,  -.5571) 


Photograph  6.27 


13° 


tilt  same  direction.  ^  =  -.509'  (-.463 


Photograph  6.32 


Photograph  6.34 


17°  tilt  same  direction.  W 


-.940'  (-.830' 


Photograph  6.36 


Photograph  6.33 


Photograph  6.40 


Photograph  6.44 


Photograph  6.4 6 


Photograph  6.47 


There  are  a  number  of  important  features  in  these  photographs. 

The  first  thing  one  notices  is  the  lack  of  distinction  between  pat¬ 
terns  produced  by  countertilting  or  tilting  in  the  same  direction. 

Coma  has  no  noticable  effect  on  the  images.  This  may  be  due  to  the 
difference  in  magnitude.  At  18°  tilt  the  amount  of  astigmatism  is 
-1.065^.  The  amount  of  coma  at  that  same  angle  is  a  mere  -.155\. 
Astigmatism  breaks  up  the  ring  structure  into  separate  nodules.  The 
disintegration  of  the  rings  begins  in  the  center  of  the  diffraction 
pattern.  The  greater  the  amount  of  astigmatism  the  greater  the  number 
of  fractured  rings. 

A  ring  will  break  up  in  a  predictable  way.  If  the  number  of  the 
ring  (counting  from  the  center)  is  n,  then  the  ring  will  split  into  4n 
nodules  and  form  a  square  of  n  +  1  nodules  per  side. 

To  describe  the  process  of  ring  break  up,  it  is  convenient  to 
define  some  new  terms.  The  first  effect  of  astigmatism  on  a  given 
ring  is  to  distort  it  from  its  circular  shape.  The  ring  is  pushed  out 
at  points  on  the  top  and  bottom  and  left  and  right.  These  points  are 
spaced  by  right  angles.  As  such,  they  suggest  the  directions  north, 
south,  east  and  west,  and  will  be  referred  to  as  compass  points.  The 
compass  point  at  the  top  of  the  photograph  (the  arbitrary  "up"  direc¬ 
tion)  will  be  designated  north.  A  ring  e  '  jiting  departure  from 
circular  shape  with  the  appearance  of  compass  points  will  be  said  to 
be  distorted. 

After  the  ring  becomes  distorted,  an  increase  in  astigmatism  will 
cause  lobing.  The  width  of  the  ring  is  no  longer  uniform,  but  concen¬ 
trated  in  regions  which  will  become  nodules  —  remnants  of  the  ring 
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when  finally  and  fully  split.  The  actual  point  of  ring  break-up  will 
be  defined  as  when  at  least  half  of  the  ring  is  established  into 
separate  nodules.  The  region  inside  the  most  newly  broken  ring  will 
be  called  the  aberrated  area. 

The  analysis  of  some  sample  photographs  will  fortify  these  con¬ 
cepts.  Photograph  6.20  (0°  tilt)  shows  the  first  ring  slightly  dis¬ 
torted.  This  is  evidence  of  presence  of  background  (positive)  astig¬ 
matism.  Photograph  6.21  is  a  very  symmetric  set  of  concentric  rings 
indicating  near  absence  of  astigmatism.  Photographs  6.22  through 
6.27  chronicle  the  increased  lobing  in  the  first  ring  to  eventual 
ring  break-up  at  13°  tilt.  The  difference  between  lobing  and  ring 
break-up  is  seen  in  photographs  6.38  and  6.39.  At  19°  tilt  the 
fourth  ring  shows  definite  lobing.  At  20°  tilt  the  fourth  ring  has 
broken  up  into  distinct  nodules.  Notice  that  for  large  amounts  of 
astigmatism  (over  1)0,  distortion  extends  far  beyond  the  region  of 
broken  rings. 

The  breakup  of  the  ring  structure  has  a  linear  dependence  on  the 
amount  of  astigmatism  at  the  outer  edge  of  the  obstacle  (w22)*  A 
review  of  the  photographs  shows  the  first  ring  disintegrates  between 
12°  and  13°  tilt.  The  table  below  lists  the  angles  for  which  dif¬ 
ferent  rings  break  up.  If  one  plots  the  ring  number  versus  the  aber¬ 
ration  coefficient  W  at  which  the  ring  breaks  up,  one  sees  the 
linear  relationship.  (The  error  bars  represent  the  +_  .5°  tilt  angle 
uncertainty.)  This  plot,  though,  is  misleading.  The  abscissa  is 
uniformly  spaced.  Spatially  the  rings  are  not  evenly  distributed. 
Plotting  the  required  astigmatism  at  ring  breakup  as  a  function  of 
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the  radial  distance  from  the  center  of  the  spot  provides  a  much 
better  perspective  on  the  domain  of  influence  of  the  aberration. 


Table  6.20 

Required  Countertilt  Angles  for  Diffraction 
Pattern  Ring  Break-up 


Ring  Number  Tilt  Angle 


(wavelengths) 


1 

2 

3 

4 

5 

6 

7 

8 


12  -  13° 
15° 

18° 

20° 

22° 

24° 

25° 

27° 

29° 


-.419  -  -.509 
-.710 
-1 .065 
-1.338 
-1.64 
-1 .97 
-2.H6 
-2.519 


9 


-2.921 


The  plotted  line  is  the  best  fit  to  the  data  by  minimizing  the  root 
mean  square  deviation.  One  can  multiply  the  radial  distance  by  the 
constant  kd/2z  and  obtain  the  argument  for  the  J,  Bessel  function  for 

the  unaberrated  case.  A  graph  of  V  as  a  function  of  the  Bessel 

22 

function  argument  appears  below.  Graphs  6.19  and  6.20  demonstrate  the 
radius  of  the  region  of  ring  break-up  expands  linearly  with  increased 


M22  UERSL1S  BESSEL  FUNCTION  ftRGtlRENT 


Graph  6.20 


astigmatism.  In  particular  when  the  aberration  coefficient  V _ ,  corre¬ 
sponds  to  a  parameter  x=  ~~  which  is  a  maximum  of  JQ(x),  the  ring  at 
distance  r  will  break  up. 

The  position  of  the  outermost  broken  ring  is  a  very  noticable  and 


convenient  measure  of  astigmatism  present  in  the  incident  beam.  There 
are,  however,  subtle  astigmatic  effects  which  bear  mention.  In  the 
unaberrated  diffraction  pattern  the  region  of  highest  intensity  is  the 


'■'»  j  »'J^ '.!  ■  j  ^  w  ^  j  .  »  ■  r  .  w*  \"  w?  . ■.  -.% 


210 


central  spot.  As  astigmatism  is  introduced  this  region  moves  out  from 
the  center.  For  example,  Photographs  6.24  through  6.29  chronicle  the 
shift  of  maximum  intensity  from  the  center  to  the  first  (broken)  ring. 
At  W  =  -.71  X  (Photograph  6.30)  the  central  spot  has  entirely  disap¬ 
peared.  Photographs  6.34  and  6.35  show  at  W22  =  -.94^  the  second 

(broken)  ring  is  becoming  the  most  intense  area.  The  central  spot  is 

reappearing  and  the  first  ring  is  dimming.  At  W22  =  —  1  •  065 (Photo¬ 
graph  6.36),  the  third  ring  has  split,  the  second  ring  is  the  most 
intense,  the  first  ring  has  disappeared  and  the  central  maximum  is 
bright.  This  pattern  continues  through  the  remaining  photographs.  As 
the  amount  of  astigmatism  increases  the  area  of  greatest  intensity 
spreads  out  from  the  center  like  the  crest  of  a  wave.  Ahead  of  this 
crest  distortion  occurs  and  immediately  in  front  of  it  ring  break-up 
takes  place.  This  crest  is  composed  of  either  one  or  two  rings. 
Immediately  behind  this  maximum  intensity  "wave,"  there  is  a  dark 
region.  The  broken  rings  left  behind  by  the  crest  form  a  checkerboard 
pattern.  The  center  and  even  rings,  and  the  odd  rings  wax  and  wane 

alternately.  All  these  areas  of  the  pattern  are  characteristic  of  the 

amount  of  astigmatism  present.  This  will  be  important  when  studying 
the  diffraction  patterns  of  annular  apertures  where  only  the  region 
around  the  central  spot  remains  intact. 

The  central  spot  first  becomes  dark  when  =  -.71  A.  The  cycle 
of  having  a  dark  center  occurs  with  added  multiples  of  -.808X  (^  .  035'-). 

An  interesting  feature  of  Photographs  6.40  -  6.49  is  the  concen¬ 
tration  of  intensity  at  the  compass  points  of  distorted  rings.  This 


concentration  extends  beyond  the  region  where  ring  break-up  takes 
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place.  The  ultimate  astigmatic  diffraction  pattern  seems  to  be  a 
bright  diamond  outline  with  spikes  at  the  compass  points. 

The  photographed  patterns  seem  to  be  slightly  canted.  Using  the 
compass  analogy,  they  don't  point  to  true  north.  A  certain  amount  of 
the  angular  deviation  may  stem  from  the  optics  of  the  Bausch  and  Lomb 
stereo  viewer.  However,  photographs  taken  using  just  a  mounted  camera 
back  also  revealed  a  slight  angle  of  tilt.  The  origin  of  this  tilt 
may  lie  in  the  background  astigmatism.  The  orientations  of  the  planes 
of  symmetry  of  the  astigmatism  produced  by  the  plane  parallel  plates 
lie  in  the  vertical  and  horizontal  directions.  The  background  astig¬ 
matism  is  not  necessarily  oriented  in  these  directions.  The  diffrac¬ 
tion  patterns  are  the  result  of  the  combination  of  these  two  aberra¬ 
tion  sources.  A  difference  in  orientations  may  create  a  tilted  pat¬ 
tern.  If  this  is  true  and  if  the  tilt  angle  is  characteristic  of  the 
individual  orientations,  this  may  provide  the  means  of  determining  an 
unknown  angle  of  orientation.  One  could  introduce  astigmatism  of 
known  orientation  and  observe  the  combination. 

The  diffraction  pattern  behind  a  circular  obstacle  for  the  unab¬ 
errated  case  does  not  change  shape  when  viewed  at  different  optical 
axis  distances.  The  pattern  with  astigmatism  present  is  also  constant 
in  shape.  The  following  photographs  show  the  diffraction  patterns  for 
three  different  distances,  .35m,  .60m  and  3m.  The  first  photograph 
was  taken  with  magnification  7X;  the  remaining  two  with  magnification 
6X.  The  patterns  are  identical. 
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Photograph  6. 51 

•60m  behind  5/16"  diameter  obstacle. 

15°  countertilt.  =  -.701 \.  1/2  sec.  exposure 


Photograph  6.52 


3m  behind  5/16"  diameter  obstacle. 

15°  countertilt.  W  =  -.701  X.  1/2  sec.  exposure. 


The  fact  the  shape  of  the  diffraction  pattern  for  the  circular 
obstacle  is  independent  of  optical  axis  distance  z  means  Graph  6.20 
(W  as  a  function  of  the  Bessel  function  argument  kdr/2z)  is  valid 
for  all  z.  One  can  determine  the  amount  of  astigmatism  present  by 
observing  the  pattern  for  a  given  obstacle  at  any  distance. 

Research  reveals  that  the  shape  of  tne  diffraction  pattern  is 
also  independent  of  obstacle  diameter.  What  determines  the  eventual 
diffraction  image  is  the  amount  of  astigmatism  at  the  edge  of  the 
obstacle  (W  ,).  The  following  photographs  of  patterns  produced  by 
the  3/16  inch  and  1/2  inch  diameter  obstacles  demonstrate  this.  The 
first  two  photographs  show  patterns  very  similar  to  Photograph  6.30 


(W  „  *  -.701 X).  The  remaining  two  photographs  correspond  closely  to 
Photograph  6.40  (W  =  -1.485^). 


If*/; 

4  T‘  ’»  A  It 


rJw*  •  «■*  •»*  '<>*<.* 

A/  a  ■■•  VM 

If- *  /*'  v  „  \  ^  % 

i^^SjgggS 


w 


'  *v  _ 

'*s> 


%**$* 


fcinre 


[♦!*] 


Photograph  6.54 


5m  behind  1/2"  diameter  obstacle.  9°  countertilt. 

-.647  1  (-.566X,  -.732X),  Magnification  4X.  1/2  sec.  exposure. 


Photograph  6.55 

•50m  behind  3/16"  diameter  obstacle.  35°  countertilt. 

1 . 485-'  (-1.44H,  -1.5311),  Magnification  bX.  1/8  sec.  exposure 


Photograph  6.56 


3m  behind  1/2"  diameter  obstacle.  13°  countertilt. 

=  -1.455X  (-1.333a,  -1.577a),  Magnification  6X.  1/2  sec.  exposure 


There  is  close  agreement  in  the  amount  of  astigmatism  for  the  three 
obstacle  diameters.  The  first  series  of  photographs  give  W  within 

t-  £ 

,05a  and  the  second  series  within  .03\.  Comparisons  of  photographs 
not  reproduced  here  gave  results  consistently  within  .051. 

There  is  an  obvious  application  of  the  above.  A  catalog  of 
diffraction  patterns  for  some  arbitrary  obstacle  diameter  with  the 
corresponding  amounts  of  astigmatism  can  be  produced  very  simply. 
(For  example,  Photographs  6.20  -  6.49  could  be  used.)  To  find  the 
amount  of  astigmatism  in  a  test  beam,  one  need  only  compare  the  dif¬ 
fraction  pattern  behind  any  sized  circular  obstacle  to  those  in  the 
catalog.  The  coefficient  W  is  readily  determined. 
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Annular  Aperture 

The  diffraction  pattern  behind  an  annular  aperture  is  a  combina¬ 
tion  of  diffraction  effects  from  the  inner  and  outer  edges.  Each 
edge  contributes  its  own  Arago  spot  and  ring  pattern.  It  is  instru":- 


Photograph  6.58 


Photograph  6.60 


1/2"  diameter  obstacle  with  5/8"  diameter  aperture 
8°  countertilt. 


Photograph  6.61 


1/2"  diameter  obstacle  with  5/8"  diameter  aperture 
9°  countertilt. 


Table  6.21 


Astigmatism  Coefficients  of  Off-set  Diffraction  Patterns 


1/2" 

obstacle 

5/16"  obstacle 

comparison 

Countertilt 

Angle 

*22 

(wavelength) 

Countertilt 

Angle 

*2  2 

(wavelength) 

6° 

-.233 

10° 

-.261 

7° 

-.352 

11° 

-.337 

8° 

-.492 

13° 

-.509 

9° 

-.647 

15° 

-.710 

10° 

-.821 

17° 

-.940 

5/6" 

aperture 

5/16"  obstacle 

comparison 

Countertilt 

Angle 

w22 

(wavelength) 

Countertilt 

Angle 

W,, 

(wavelength) 

6° 

-.419 

12° 

-.419 

7° 

-.606 

14° 

-.606 

9° 

-1 .065 

19° 

-1.198 

diffraction  patterns  from  the  aperture  for  countertilt  angles  of 
and  10°  could  not  be  determined.)  The  individual  patterns  are 


Photograph  6.64 

3/16  diameter  obstacle  with  5/8"  diameter  aperture, 
e  =  .3.  Viewing  distance:  2  -  2.2m.  Magnification  4X. 
Exposure  time  1/4  second.  7°  countertilt.  Maximum  on  axis 


Photograph  6.65 

3/16"  diameter  obstacle  with  5/8"  diameter  aperture. 
e  -  .3.  Viewing  distance:  2  -  2.2m.  Magnification  4X. 
Exposure  time  1/4  second.  10°  countertilt.  Minimum  on  axis. 
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Photograph  6.68 

3/16"  diameter  obstacle  with  5/8"  diameter  aperture, 
e  =  .3.  Viewing  distance:  2  -  2.2m.  Magnification  4X. 
Exposure  time  1/4  second.  12°  countertilt.  Maximum  on  axis 
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Photograph  6.69 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture. 

;  =  .625.  Viewing  distance:  2  -  2.05m.  Magnification  6X 
Exposure  time  1  second.  10°  countertilt.  Minimum  on  axis 


Photograph  6.70 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture. 

=  .625.  Viewing  distance:  2  -  2.05m.  Magnification  6X 
Exposure  time  1  second.  11°  countertilt.  Maximum  on  axis 


Photograph  6.71 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture. 


=  .625.  Viewing  distance:  2  -  2.05m.  Magnification  bX 
Exposure  time  1  second.  12°  countertilt.  Maximum  on  axis 


Photograph  6.72 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture, 
e  =  .625.  Viewing  distance:  2  -  2.05m.  Magnification  6X 
Exposure  time  1  second.  13°  countertilt.  Minimum  on  axis 


Photograph  6.73 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture, 
e  =  .625.  Viewing  distance:  2  -  2.05m.  Magnification  6X 
Exposure  time  1  second.  13°  countertilt.  Maximum  on  axis 
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Photograph  6.74 

5/16"  diameter  obstacle  with  1/2"  diameter  aperture. 

=  .625.  Viewing  distance:  2  -  2.05m.  Magnification  6X 
Exposure  time  1  second.  14°  countertilt.  Minimum  on  axis 
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Photograph  b. / 5 

5/16”  diameter  obstacle  with  1/i"  diameter  aperture. 

=  .625*  Viewing  distance:  2  -  ~.G5m.  Magnification  6X 
Exposure  time  1  second.  14°  countertilt.  Maximum  on  axis 
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Photograph  6.78 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 

=  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X. 


Exposure  time  1  second. 


countertilt.  Minimum  on  axis. 
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Photograph  6.79 

3/16"  diameter  obstacle  with  7mm  diameter  aperture, 
f  »  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  22°  countertilt.  Maximum  on  axis 


Photograph  6.80 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 
e  =  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  24°  countertilt.  Minimum  on  axis 


Photograph  6.81 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 

£  =  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  24°  countertilt.  Maximum  on  axis 


Photograph  6.32 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 

£  =  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  26°  countertilt.  Minimum  on  axis 


Photograph  6.83 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 

=  .63.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  26°  countertilt.  Maximum  on  axis 


Photograph  6.84 

3/16"  diameter  obstacle  with  7mm  diameter  aperture. 
r.  =  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  28°  countertilt.  Maximum  on  axis 


Photograph  6.85 

3 / 1  6 "  diameter  obstacle  with  7mr.  diameter  aperture. 

>  =  .68.  Viewing  distance:  .6  -  .65m.  Magnification  4X 
Exposure  time  1  second.  36°  countertilt.  Maximum  on  axis 


Photograph  6.b7 

1/2"  diameter  obstacle  with  5/b"  diam-: 
•  =  .3.  Viewing  distance:  2.J5  -  2.15m. 
Exposure  time  1  second.  10°  countertilt. 


ter  aperture. 
Magnification  6X. 
Minimum  on  axis. 
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Photograph  6.88 

1/2  diameter  obstacle  with  5/8"  diameter  aperture. 

-  .8.  Viewing  distance:  2.05  -  2.15m.  Magnification  6X. 
xposure  time  1  second.  10  countertilt.  Maximum  on  axis. 
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Photograph  6.39 

1/2"  diameter  obstacle  with  5/8”  diameter  aperture. 

=  .8.  Viewing  distance:  2.05  -  2.15m.  Magnification  6X. 
xposure  time  1  second.  12  countertilt.  Maximum  on  axis. 
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The  diffraction  patterns  behind  annular  apertures  with  high  ob¬ 
scuration  ratios  are  quite  similar  to  those  produced  by  a  circular 
obstacle.  The  resemblance  decreases  as  the  obscuration  ratio  does. 

The  outer  areas  of  the  aberrated  region  cease  to  show  lobing  and  ring 
break-up  but  give  way  to  intricate  symmetries.  When  the  obscuration 
ratio  is  as  low  as  .3  none  of  the  structures  seen  in  the  circular 
obscuration  case  are  manifested. 

An  effect  seen  in  the  unaberrated  case  seems  to  be  exaggerated 
when  astigmatism  is  present.  This  is  the  dramatic  change  in  the 
overall  pattern  from  the  position  at  a  maximum  on  axis  to  a  minimum  on 
axis.  Photographs  6.65  and  6.66,  6.67  and  6.68,  6.74  and  6.75,  and 
6.80  and  6.81  are  good  examples  of  this  effect.  Because  the  patterns 
can  vary  to  such  a  large  extent  a  catalog  of  diffraction  patterns,  as 
suggested  for  the  circular  obscuration,  is  impractical.  The  changes 
go  through  cycles  corresponding  to  the  cycles  of  on-axis  extrema.  The 
pattern  at  any  viewing  distance  z  will  be  associated  with  a  given 
point  in  the  cycle  and  so  the  shape  of  the  pattern  will  be  repeated 
many  times  as  one  moves  along  the  optical  axis. 

This  change  from  maximum  to  minimum  must  be  taken  into  account 
when  comparing  the  annular  aperture  photographs  to  the  circular  obsta¬ 
cle  photographs.  For  example  Photograph  6.87  shows  the  pattern  behind 
an  annular  aperture  of  obscuration  ratio  .8  at  an  on-axis  minimum. 

(it  is  clear  this  is  a  minimum  on  axis  because  the  number  and  extent 
of  tne  bright  rings  in  the  center  is  maximum.  Compare  this  to  Photo¬ 
graph  6.19.)  The  prominent  central  feature  is  the  absence  of  the 
bright  central  spot  and  the  four  very  bright  nodules  composing  the 
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split  first  ring.  This  is  reminiscent  of  Photograph  6. 50  corresponding 
to  W,,2  =  -.71'..  However,  this  is  incorrect.  The  last  broken  ring  in 
Photograph  6.87  is  the  third  and  the  fourth  ring  shows  severe  lobing. 
This  is  suggestive  of  Photograph  6.53  corresponding  to  W  =  -1.198'.. 

A  comparison  of  Photograph  6.88  to  the  circular  obstacle  photographs 
confirms  this.  In  this  photograph  the  center  spot  is  bright  and  the 
similarity  to  Photograph  6.58  is  immediately  apparent.  A  good  rule  of 
thumb  when  comparing  photographs  of  annular  aperture  diffraction  pat¬ 
terns  to  circular  obscuration  diffraction  patterns  is  to  count  broken 
rings  if  the  pattern  is  at  an  on-axis  minimum  and  match  interior 
patterns  when  at  an  on-axis  maximum.  Table  6.22  lists  the  results  of 
a  comparison  of  some  photographs  from  6.70  -  6.90  to  photographs  6.20 
-  6.49.  The  figures  in  parentheses  represent  the  high  and  low  limits 
to  W  based  upon  an  uncertainty  in  tilt  angle  measurement  of  +_  .5°. 

The  preceeding  table  shows  that  the  astigmatism  at  the  outer  edge  of 
the  annular  aperture  drives  the  diffraction  pattern. 

Apertures  of  similar  obscuration  ratios  exhibit  similar  diffrac¬ 
tion  patterns  with  similar  amounts  of  astigmatism.  Photograph  6.74 
corresponds  to  e  =  .625  and  W,,„  =  -1.704).  at  the  outer  edge.  This 
pattern  compares  very  closely  to  Photograph  6.81  with  t  =  .68  and  W„., 

*  -1.51X  at  the  outer  edge.  Photograph  6.72  with  =  -1.97'  is 

similar  to  Photograph  6.85  with  W  =  -1.79'- 

£-  «i 

It  was  suggested  in  Chapter  5  that  annular  apertures  with  high 
obscuration  ratios  are  better  suited  for  aberration  studies.  This  is 
due  to  the  larger  number  of  rings  in  the  beat  envelopes.  More  of  the 
ring  structure  stays  intact  between  modulations  so  any  deviations  in 


Table  6.22 
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the  ring  pattern  is  easier  to  spot.  This  raises  the  question,  can  the 
amount  of  astigmatism  be  determined  from  patterns  of  an  annular  aper¬ 
ture  with  a  low  obscuration  ratio?  The  answer  is  a  guarded  yes.  A 
rough  measure  of  the  aberration  present  is  possible.  The  region  of 
aberration  extends  a  distance  out  from  the  center  of  the  diffraction 
pattern.  This  radial  distance  can  be  measured  and  compared  to  the 
radial  distance  for  ring  breakup  found  in  the  previous  section. 

To  prepare  the  following  table,  the  average  radii  of  the  regions 
of  aberration  for  photographs  of  the  3/1 6  inch  diameter  obstacle  with 
a  5/8  inch  diameter  aperture  were  measured.  (The  aberrated  region  was 
taken  as  the  area  out  to  the  last  bright  band  which  showed  looping  or 
separation  in  the  band.  For  example,  in  Photograph  6.66  this  was  the 
third  band  from  the  center.)  These  were  compared  to  the  known  ring 


Table  6.23 

Radii  of  Aberrated  Areas  in  Diffraction  Patterns 
Behind  an  Annular  Aperture  with  e  =  .3 


Tilt  angle 

W22  at 

outer  edge  (*) 

Radius  of  aberrated 
area  (mm) 

kdr 

2z 

7 

-.606  (-.509, -.710) 

.17 

6.69 

8 

-.821  (-.710, -.940) 

.195 

7.68 

10 

-1.353  (-1.193,-1.485) 

.33 

13.00 

12 

-1.970  (-1.801,-2.146) 

.45 

17.73 

14 

-2.717  (-2.519,-2.921  ) 

•  56 

22.10 

(innn 
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spacing  for  the  annular  aperture  with  no  aberrations  present.  The 
actual  radial  distance  of  the  extent  of  the  region  of  aberration  was 
calculated.  The  Bessel  function  argument  kdr/2z  was  calculated  also. 

The  results  from  the  table  are  plotted  against  the  line  found  in 
Graph  6.20.  The  points  suggest  a  line  of  steeper  slope  than  that  from 
Graph  6.20.  Recall  the  line  plotted  in  Graph  6.2C  corresponded  to  the 
radius  of  the  aberrated  area  when  defined  as  the  region  of  broken 
rings.  The  distance  which  correlates  to  this  radius  in  the  annular 
aperture  case  seems  to  be  different  from  that  measured  in  Table  6.  23. 

A  new  line  with  a  different  slope  should  enable  one  to  determine  1#'2 , 
from  measurements  of  the  diffraction  patterns  of  an  annular  aperture 
with  a  low  obscuration  ratio. 


H22  VERSUS  BESSEL  FUNCTION  ARGUMENT  (hdr/2z:  r  :  ABERRATE*  AREA  RAIL'S) 


Graph  6.21 


2. 


An  interesting  feature  is  seen  in  the  photographs  for  the  annular 
aperture.  The  center  of  the  diffraction  pattern  often  contains  a  +  or 
x  structure.  These  are  seen  quite  readily  in  the  photographs  below. 


Photograph  6.90 


3/16"  diameter  obstacle  with  5/8"  diameter  aperture. 
Viewing  distance  2m.  8°  countertilt. 


Photograph  6.91 


3/1t>"  diameter  obstacle  with  5/8"  diameter  aperture. 
Viewing  distance  2.05m.  8°  countertilt. 


The  structures  are  also  seen  in  Photographs  6.63,  6.64,  6.67,  6.63, 
6.78  and  6.79.  In  these  last  two  instances  the  +  and  x  are  dark 
structures  at  the  center.  These  seem  to  be  related  to  the  positions 
of  maxima  and  minima  on  axis.  A  bright  ♦  appears  to  be  composed  of 
the  bright  central  spot  and  the  first  broken  ring.  A  bright  x  looks 
to  be  made  of  the  bright  central  spot  and  the  northeast,  southeast, 
southwest  and  northwest  nodules  in  the  second  broken  ring.  Minimum 
intensities  in  these  areas  give  rise  to  the  dark  structures. 

Astigmatism,  like  spherical  aberration  and  defocus,  shifts  the 
position  of  on-axis  extrema.  This  was  verified  using  the  same  proce¬ 
dure  as  for  tne  previously  mentioned  aberrations.  The  first  run 
involved  the  shift  of  the  7tn  order  minimum  at  65cm  behind  an  annular 


aperture  composed  of  a  3/16  inch  diameter  obstacle  with  a  7mm  diameter 


aperture  separated  by  1cm.  The  seccr.d  run  involves  tne  shift  of  the 
23rd  order  minimum  at  66.5cm  behind  a  3/1 6  inch  diameter  obstacle  with 
a  10mm  diameter  aperture  separated  by  .3cm.  The  graphs  below  plot  the 
data  points  for  each  run.  The  line  is  the  line  of  best  fit  by  mini¬ 
mizing  the  root  mean  square  deviations. 


ASTIGMATISM:  3/16  INCH  DIAMETEP.  OBSTACLE  WITH  7  MH  DIAMETER  APERTURE 


Graph  6.22 


ASTIGMATISM:  3/16  INCH  DIAMETER  OBSTACLE  WITH  18  MM  DIAMETER  APERTURE 


Graph  6.23 


The  coefficients  W  from  Graph  6.23  seems  to  go  as  the  square  of  the 


axial  shift.  However,  the  linear  plots  have  an  attractive  feature. 

The  ratio  of  the  radii  for  the  two  apertures  used  is  .7.  The  ratio  of 
slopes  of  the  best  fit  lines  in  the  graphs  is  .723*  The  data  does  not 
fit  a  formula  like  that  found  for  defocus  and  spherical  aberration: 


W  =  £  A* 

22  C  z' 


(VI-12) 


where  n  is  the  order  of  the  minimum  and  c  is  a  constant  presumably 
dependent  upon  r. 

The  complex  astigmatic  patterns  produced  wit:,  high  angles  of 
countertiit  make  measurements  of  on-axis  shifts  very  difficult.  A 
more  conclusive  test  of  the  effect  of  astigmatism  on  the  positions  of 


axial  extrema  should  limit  itself  to  small  amounts  of  the  aberration. 


For  now  it  is  sufficient  to  note  negative  astigmatism  causes  a  shift 
on  axis  in  the  negative  z  direction. 

There  are  two  reasons  not  to  use  the  shift  of  extrema  to  r.easur 
astigmatism.  The  first  reason  is  that  the  change  in  the  overall 
diffraction  pattern  is  an  accurate  and  sufficient  gauge  of  aberratio 
present.  The  second  is  that  both  defocus  and  spherical  aberration, 
which  have  little  effect  on  the  pattern,  cause  axial  shifting,  ’’en¬ 
suring  the  shift  alone  could  lead  to  an  erroneous  measure  for  W  -  if 
W20  and  are  not  accounted  for. 


i 

t 
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Chapter  7 


CONCLUSIONS 

The  most  important  result  of  the  theoretical  and  experimental 
treatment  in  this  thesis  is  almost  trivial.  The  ultimate  diffraction 
pattern  behind  a  circular  obstacle  or  an  annular  aperture  depends 
entirely  on  the  physical  state  of  the  incident  wave  at  the  edges. 

This  is  seen  repeatedly  in  the  formulas  from  chapters  3  and  4-  The 
disturbance  function  at  points  in  image  space  contains  phase  terms, 
Gaussian  terms,  aberration  terms  and  zero-order  Bessel  functions  which 
all  have  arguments  dependent  upon  the  radii  of  the  edges.  By  under¬ 
standing  the  physics  at  these  boundaries  one  gains  insight  into  the 
diffraction  process  and  its  results. 


Summary 


The  diffraction  patterns  predicted  by  the  theory  of  chapters  3 
and  4  agreed  very  closely  to  the  experimental  findings,  typically  with 
less  than  5%  deviation.  The  intensity  distribution  in  the  near  field 
region  close  to  the  optical  axis  can  be  described  with  excellent 
accuracy  by  the  formulas  reviewed  at  the  end  of  Chapter  4.  These 
formulas  depend  only  on  the  zero  order  Bessel  function  so  the  computa¬ 
tion  of  the  intensity  at  a  given  point  is  substantially  reduced  from 
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that  requiring  the  use  of  unwieldy  Lommel  functions.  The  periodicity 
of  the  Bessel  function  provides  an  explanation  of  the  modulated  ring 
structure  observed  in  annular  aperture  diffraction  patterns  in  terms 
of  beat  envelopes. 

The  theory  of  Chapter  4  was  very  successful  in  predicting  the 
shift  of  on-axis  extrema  behind  an  annular  aperture  with  the  introduc¬ 
tion  of  defocus  in  the  incident  beam.  Measuring  this  shift  is  a  very 
effective  method  of  determining  defocus.  This  method  should  work  for 
spherical  aberration  also,  but  the  maximum  amount  of  spherical  aberra¬ 
tion  created  from  the  apparatus  produced  expected  shifts  less  than  the 
experimental  error.  It  was  found  that  astigmatism  will  shift  the 
positions  of  the  axial  maxima  and  minima  as  well.  In  all  the  above 
cases,  if  the  aberration  coefficient  is  negative  the  shift  will  be  in 
the  minus  z  direction. 

Astigmatism  produces  very  distinct  changes  in  the  diffraction 
patterns  of  a  circular  obstacle.  The  effect  of  the  aberration  is 
seen  in  the  first  ring  of  diffraction  pattern  at  values  of  Wn m  as  low 
as  -.lx.  This  is  when  the  first  ring  starts  to  concentrate  its 
intensity  at  four  compass  points  —  the  vertices  of  a  diamond.  As 
the  aberration  is  increased,  the  first  ring  breaks  up  into  four 
separate  spots  when  W  is  approximately  -.4,\.  In  general  the  nl< 
ring  from  the  center  will  break  up  into  4n  points  forming  a  diamond 
with  n+1  points  on  a  side.  This  region  of  broken  rings  constitutes 
an  area  where  aberration  effects  are  manifest.  From  the  radius  of 
this  region,  the  amount  of  astigmatism  can  be  calculated.  The  radius 
and  W_^  are  related  linearly.  The  same  amount  of  astigmatism  at  the 
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edge  of  any  sized  circular  obstacle  will  produce  the  same  aberrated 
image . 

An  annular  aperture  with  a  high  obscuration  ratio  (  >  .6)  sub¬ 
jected  to  an  astigmatic  beam  will  give  a  diffraction  pattern  similar 
to  that  for  the  circular  obstacle.  In  this  case  it  is  the  amount  of 
astigmatism  at  the  outer  edge  that  determines  the  pattern.  Apertures 
of  nearly  equal  obscuration  ratios  produce  similar  diffractions  with 
similar  amounts  of  astigmatism  at  the  outer  edges.  W?0  can  be  calcu¬ 
lated  by  comparing  the  annular  aperture  diffraction  pattern  to  those 
for  the  circular  obstacle.  When  the  obscuration  ratio  is  low,  the 
diffraction  pattern  is  very  complex.  In  this  case  one  can  measure 
the  extent  of  the  aberrated  region  and,  in  an  analogous  manner  to 
that  for  the  circular  obstacle,  find  the  amount  of  aberration  present. 

The  astigmatism  produced  by  tilting  plane  parallel  plates  drowned 
out  any  effects  which  might  have  come  from  coma.  A  numerical  study  of 
coma  is  necessary.  Coma  should  produce  a  diffraction  image  with  one 
plane  of  symmetry. 

The  circular  obstacle  seems  best  suited  to  detect  the  presence  of 
astigmatism  and,  likely,  coma.  An  annular  aperture  is  required  to 
detect  the  rotationally  symmetric  aberrations,  defocus  and  spherical 
aberration. 

Many  of  the  theoretical  and  experimental  results  found  in  this 
thesis  suggest  embarkation  points  for  future  study. 
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Suggestions  for  Further  Research 

The  experimental  apparatus  did  not  include  any  method  of  mea¬ 
suring  the  intensity  of  Arago's  spot  and  the  surrounding  ring  struc¬ 
ture.  Many  theoretical  predictions  from  chapters  3  and  4  involve  the 
intensity  of  the  spot.  For  example,  measuring  the  intensities  of  two 
on-axis  maxima  will  allow  one  to  separate  defocus  and  spherical  aber¬ 
ration.  A  detector  to  measure  this  region  near  the  optical  axis  must 
be  very  small.  Its  dimensions  will  be  guided  by  the  dimensions  of  the 
aperture.  The  smaller  the  aperture,  the  larger  the  detector  can  be. 

The  detector  can  be  made  by  masking  all  but  a  small  area  (diameter  on 
the  order  of  10  to  50  microns)  on  the  detector  surface.  The  detector 
should  allow  translation  along  the  z-axis.  The  detector  should  be 
allowed  to  move  in  a  radial  direction  to  find  the  intensity  peak. 

This  will  allow  study  of  the  ring  structure  as  well. 

The  experimental  apparatus  contained  some  background  astigmatism. 
The  measure  of  the  tilt  angles  had  a  half-degree  uncertainty.  A 
definitive  study  of  the  astigmatic  case  where  the  tilt  angles  can  be 
measured  to  high  accuracy  should  be  made.  A  catalog  of  diffraction 
patterns  can  be  compiled  and  used  in  aberration  detection. 

Another  failing  in  the  experimental  equipment  was  the  size  of  the 
aperture  stop  at  the  collimating  lens  ( 1 3 / 1 6  inch).  If  the  parallel 
plates  were  tilted  in  the  same  direction  the  transmitted  beam  began  to 
hit  the  collimating  lens  off  center.  If  the  tilt  angle  was  much 
greater  than  20°  the  beam  began  to  miss  tie  lens.  A  larger  collima¬ 
ting  lens  and  holder  should  be  used  in  future  set-ups. 


The  Spot  of  Arago  has  a  long  and  interesting  history,  but  so  far 
its  practical  applications  have  been  few.  The  use  of  the  spot  to 
detect  and  measure  aberrations  will  renew  its  importance. 
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Appendix  1 

THE  LIFE  AND  WORKS  OF  ARAGO 

The  great  men  and  women  of  science  are  known  typically  only  by 
their  theorems,  equations  or  experiments.  Textbooks  rarely  have  the 
space  to  provide  biographical  information  beyond  an  occasional  foot¬ 
note  containing  dates  of  birch  and  death  or  an  amusing  anecdote.  A 
sense  of  history,  the  personal  side  of  science,  is  missed  in  the 
process.  Many  scientists  led  colorful  and  controversial  lives.  It 
would  be  a  particular  shame  to  overlook  the  life  of  Dominique  Francois 
Jean  Arago. 

Arago  was  born  on  February  26,  1786  at  Estagel,  a  small  village 
in  the  eastern  Pyrenees  near  the  Mediterranean  coast.  The  village  was 
a  halting  station  for  French  troops  enroute  to  the  army  of  the  Pyre¬ 
nees.  Young  Arago  was  fascinated  by  the  soldiers  and  once  asked  an 
officer  what  he  needed  to  do  to  be  able  to  wear  such  a  fine  uniform. 
The  reply  was  to  become  a  student  at  the  Polytechnic  School/  To 
prepare  for  the  entrance  examination  Arago  taught  himself  mathematics 
from  the  works  of  Euler,  Legendre,  Lagrange  and  Laplace.  He  so  im¬ 
pressed  his  examiner,  Louis  Monge,  that  he  was  entered  into  the  school 

2 

in  1805  with  highest  commendation. 

In  1804  he  was  appointed  as  secretary  to  the  Observatory  of 
Paris.  Here  he  met  and  worked  with  Laplace  and  Biot.  Arago  and  Biot 


were  commissioned  in  1806  to  complete  the  measurement  of  the  meri¬ 
dional  arc.  (The  National  Convention  earlier  had  adopted  the  measure 
of  the  one  ten-millionth  part  of  the  quadrant  of  the  meridian  passing 
through  Paris  as  the  length  of  the  meter.  The  complete  quadranc  could 
not  be  measured  but  it  was  felt  by  measuring  the  arc  from  Dunkirk 
through  Barcelona  to  the  Balearic  Islands  and  knowing  the  relative 
latitudes  the  total  quadrant  length  could  be  deduced.)^  The  two 
Frenchmen,  with  the  aid  of  two  Spanish  scientists  Chaix  and  Ro.riguez, 
first  measured  the  distance  from  the  coast  of  Spain  to  the  island  of 
Ibize.  Here  Biot  and  Arago  spent  months  on  high  mountaintops  sending 
signals  by  fire  to  their  Spanish  counterparts  on  peaks  over  a  hundred 
miles  distant.  They  braved  severe  windstorms,  bandits,  fatigue  and 
scarce  supplies.  Biot  returned  to  Paris  in  1807  with  the  preliminary 
measurements.  Arago  continued  measurements,  moving  to  the  summit  of 
Mount  Galatzo.^ 

In  1808  Napoleon  put  his  brother  Joseph  on  the  throne  of  Spain 
which  precipitated  a  war.^  The  natives  of  the  Spanish  owned  islands 
where  Arago  was  positioned  became  convinced  the  strange  fire  signals 
on  their  mountain  peaks  at  night  were  messages  meant  for  an  invading 
army.  A  mob  of  islanders  was  dispatched  to  capture  the  hapless 
Frenchman.  An  assistant  got  warning  to  Arago  who  was  able  to  disguise 
himself  as  a  peasant.  He  escaped  by  joining  the  ranks  of  the  crowd  in 
pursuit  of  him.  He  entrusted  his  equipment  and  papers  to  his  still- 
faithful  servants.  He  made  his  way  to  Palma,  the  island's  port. 
Although  he  had  a  ship  waiting  to  take  him  to  safety  he  stayed  until 
his  precious  recordings  and  equipment  were  returned.  By  this  time  the 
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mob  had  discovered  his  charade  and  to  avoid  their  wrath  Arago  was 
obliged  to  turn  himself  in  to  the  local  authorities.  He  spent  the 
next  two  months  in  prison  and  used  the  time  to  complete  calculations 
on  his  measurements.^ 

In  July  of  1808,  Rodriguez  interceded  on  Arago's  behalf  and  ob¬ 
tained  permission  for  his  passage  to  Algiers.  In  Algiers  he  befriend¬ 
ed  the  Dey,  the  governor  of  the  city.  Soon  the  young  Frenchman  set 
sail  for  Marseilles.  The  ship's  cargo  included  two  lions  meant  as 
presents  to  the  emperor  Napoleon.  Within  sight  of  the  French  coast,  a 
Spanish  privateer  captured  the  vessel.  A  casualty  of  the  seizure  was 
one  of  the  lions.  Arago  was  returned  to  Spain  and  imprisoned  for 
three  months.  The  Dey  learned  of  the  capture  and  the  lion's  demise 
and  was  outraged.  He  threatened  war  if  the  crew  and  ship  were  not  set 
free.  Spain  had  its  hands  full  with  France  and  was  not  anxious  to 

enter  into  any  other  entanglements.  The  ship  departed  once  more  for 
7 

Marseilles. 

Before  the  vessel  made  port,  a  violent  storm  forced  it  back  out 
to  sea  and  onto  the  craggy  coast  of  Sardinia.  Sardinia  was  at  war 
with  Algiers  and  refused  to  let  the  ship  land  to  make  repairs.  The 
fractured  craft  found  its  way  to  Bougia,  a  port  in  northern  Africa. 

The  winds  were  unfavorable  for  passage  to  Algiers  so  Arago  decided  to 
travel  by  land.  To  make  the  perilous  journey  he  donned  the  disguise 
of  a  Bedouin.  Before  Arago  arrived  in  Algiers  the  Dey  was  overthrown. 
The  new  governor  was  not  friendly  to  the  French  and  the  young  scien¬ 
tist  avoided  incarceration  only  by  the  intercession  of  the  Danish 
consul.  Still  he  suffered  through  a  series  of  maddening  delays. 


Finally  in  July  of  1809,  three  years  after  he  and  Biot  departed,  Arago 

Q 

returned  to  his  native  land. 

As  a  reward  for  his  diligence  and  adventures  in  the  name  of 
science  Arago  was  nominated  for  entry  into  the  Academie  des  Sciences 
when  a  vacancy  arose.  His  competition  was  Simeon  Denis  Poisson  who 
was  being  sponsored  by  Pierre  Simon  Laplace.  Joseph  Louis  Lagrange 
led  the  support  for  Arago  and  eventually  convinced  Laplace  to  change 
his  endorsement.  Arago  won  and  at  age  23  became  the  youngest  member 
ever  to  be  elected  to  the  prestigious  body.^ 

The  scientific  career  of  Arago  was  quite  varied,  involving 
theoretical  works,  experimentation,  teaching  and  administration.  He 
taught  at  the  Ecole  Polytechnic.  As  an  astronomer  at  the  Royal  Obser¬ 
vatory  he  delivered  public  lectures  from  1812  to  1845.  In  1816  Arago 
and  Joseph  Louis  Gay-Lussac  founded  and  co-edited  the  publication 
Annalen  de  Chimie  et  de  Physique.  He  also  helped  establish  Compte 
Rendu3.  When  Jean  Baptiste  Joseph  Fourier  died  in  1830  Arago  took 
over  as  permanent  secretary  of  the  Academie  des  Sciences.1^  Part  of 
hi3  duties  involved  writing  eloges  on  members  of  the  Academie.  Much 
of  what  is  known  of  the  lives  of  Augustin  Fresnel,  Thomas  Young, 
Alessandro  Volta,  Andre  Marie  Ampere,  Fourier,  Nicolas  Leonard  Sadi 
Carnot  and  others  stem  from  these  writings.^ 

Arago  was  very  active  politically.  From  1830  to  1844  he  served 
in  the  Chamber  of  Deputies.  During  his  term  he  supported  the  develop¬ 
ment  of  steam  engines,  railroads  and  lighthouses.  He  suggested  a 
patent  law,  boring  the  Artesian  wells  at  Crenelle  and  the  construction 
of  the  museum  at  Cluny.  To  promote  science  he  sponsored  Louis  Jacques 
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Mande  Dageurre  and  the  new  field  of  photography  and  arranged  for  the 
publication  of  the  works  of  Pierre  de  Fermat  and  Laplace.  In  1848 
France  became  a  republic.  Arago  was  appointed  the  Minister  of  War  and 
Marine  (a  unique  combination  of  two  offices).  As  Minister  he  put  an 
end  to  the  practice  of  flogging,  improved  naval  rations  and  abolished 
slavery  in  French  colonies.  A  staunch  republican,  Arago  refused  to 
sign  an  obligatory  oath  of  allegiance  to  the  new  emperor  Napoleon  III 

in  1852.  The  emperor  made  an  exception  in  Arago's  case  and  the  scien- 

1 2 

tist  maintained  his  offices  until  his  death  on  October  2,  1853* 

The  collected  writings  and  scientific  papers  of  Arago  occupy  ten 
volumes.  Below  is  a  short  account  of  some  of  his  scientific  achieve¬ 
ments.  In  1810  he  discovered  polarization  of  light  when  passed 

through  sheets  of  mica.  He  and  Fresnel  showed  polarized  rays  can 

1  "5 

interfere  only  when  polarized  in  transverse  planes.  He  demonstrated 
that  refraction  of  light  is  independent  of  the  motion  of  the  refract¬ 
ing  body.  This  prompted  Fresnel  to  suggest  the  Earth  drags  ether  with 
it.^  In  1833  Arago  proposed  the  comparison  of  velocities  of  light  in 
different  media.  The  corpuscular  theory  predicted  light  should  travel 
more  rapidly  in  a  medium  such  as  water  than  in  air.  The  wave  theory 
suggested  the  opposite.  Light  would  be  sent  through  two  tubes,  one 
filled  with  air,  the  other  with  water  and  onto  a  rapidly  rotating 
mirror.  The  difference  in  transit  times  through  the  tubes  would 
result  in  different  incident  angles  upon  the  rotating  mirror.  Jean 
Bernard  Leon  Foucault  conducted  the  experiment  and  found  light  moved 
more  slowly  in  water.  This  was  powerful  evidence  in  favor  of  the  wave 
theory.^  He  studied  the  bands  of  Saturn  and  Jupiter  as  well  as  the 
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polar  ice  caps  of  Mars.  He  made  photometric  measurements  of  numerous 
astronomical  bodies.  Arago  showed  a  correlation  between  the  changes 
in  the  earth's  magnetic  field  and  the  appearance  of  the  aurora  bore¬ 
alis.  He  was  the  first  to  demonstrate  that  steel  could  be  permanently 
magnetized  by  using  an  electric  current  and  iron  filings  temporarily 
so.  The  Royal  Society  of  London  awarded  Aragc  the  Copley  Medal  in 
1825  for  his  exhibition  of  magnetism  in  non-ferrous  materials.  (Other 
recipients  of  the  award  include  Humphrey  Davy.  Hans  Christian  Oersted, 
Benjamin  Franklin,  Joseph  Priestley  and  Poisson.)^ 

In  the  course  of  his  scientific  work  Arago  had  many  collabora¬ 
tors.  Mention  has  been  made  about  his  work  with  Young,  Poisson, 
Fresnel,  Dageurre,  Foucault  and  Biot.  His  relation  with  Biot  soured 
when,  while  Arago  was  on  a  trip  away  from  Paris,  Biot  attempted  to 

publish  results  of  a  joint  effort  under  only  his  name  with  but  a 

17 

cursory  reference  to  Arago's  help.  ‘  From  1824  to  1830  Arago  and 

Pierre  Louis  Dulong  performed  a  series  of  dangerous  experiments  on  the 

elastic  force  of  steam.  The  tests  involved  heating  boilers  to  the 
1  ft 

point  of  bursting.  He  was  the  mentor  for  many  scientists  in  the 
early  nineteenth  century.  In  addition  to  Foucault  and  Daguerre,  Arago 
sponsored  Urban  Jean  Joseph  Leverrier  (who  discovered  Neptune  by 
studying  pertubation  on  the  planet  Uranus)  and  Arraand  Fizeau.^ 

One  intriguing  collaboration  never  transpired.  After  Napoleon 
Bonaparte  suffered  his  defeat  at  Waterloo,  he  decided  to  leave  Europe 
and  travel  to  America  to  devote  the  remainder  of  his  life  to  science. 
He  wanted  a  companion  well  versed  in  the  sciences  to  travel  with  him. 
His  choice  was  Arago.  Arago,  though,  bordered  on  chauvinism  in  the 
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traditional  sense  of  the  word.  He  refused  Napoleon’s  offer  because  he 
wanted  France  to  remain  in  the  forefront  of  science.  The  presence  of 
Napoleon  in  America  would  sanction  the  growing  field  of  science  there 
which  could  only  diminish  the  French  scientific  authority. 

Arago's  deeply  felt  sense  of  patriotism  is  his  most  telling 
characteristic.  It  was  the  impetus  for  his  entry  into  the  Ecole 
Polytechnic,  the  launching  point  for  his  scientific  career.  It  ren¬ 
dered  inevitable  his  political  service  to  his  country.  His  extraordi¬ 
nary  life  and  works  are  extensions  of  that  spirit. 


.  A-i u: 
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Appendix  2 

THE  SPOT  OF  ARAGO  IN  THE  SHADOW  DURING  A  SOLAR  ECLIPSE 

During  a  solar  eclipse  the  moon  acts  as  an  obstacle  to  the  essen¬ 
tially  plane  waves  coning  from  the  sun.  The  configuration  is  similar 
to  the  circular  obstacle  case  described  in  Chapter  3«  Does  this  mean 
the  Spot  of  Arago  should  be  present  in  the  lunar  shadow  during  a  solar 
eclipse?  Certainly  theory  predicts  its  existence,  but  this  phenomenon 
has  never  been  recorded.  Some  possible  explanations  for  this  include 
the  spectral  distribution  of  light  from  the  sun,  the  difference  in 
spot  intensity  between  a  spherical  and  circular  obstacle,  atmospheric 
effects,  and  spot  sice. 

The  light  from  the  sun  is  not  monochromatic.  It  produces  a  dim 
multi-colored  spot  (see  Potter's  description  of  the  spot  in  Chapter  2) 
which  for  small  obstacles  must  be  viewed  through  a  magnifying  system. 
Because  the  spot  produced  by  sunlight  is  so  insubstantial  no  one  had 
cause  to  look  for  it  until  Arago.  If  this  carries  over  into  the  case 
of  a  very  large  obstacle  illuminated  by  sunlight  and  observed  at  very 
large  distance  then  the  situation  may  be  similar  to  that  in  1818. 
Because  the  spot  is  unnoticible  to  the  naked  eye  no  one  has  looked  for 
it. 

The  intensity  of  the  spot  is  also  dependent  on  whether  the  obsta¬ 
cle  was  circular  or  spherical  (see  the  section  on  Raman  and  Krishnan's 
work  in  Chapter  2).  For  example,  at  50cm  behind  a  1/4  inch  diameter 


obstacle  the  spot  intensity  for  a  spherical  obstacle  is  half  that  for 
a  circular  obstacle.  At  10cm  the  spherical  obstacle  intensity  is  nine 


times  less  that  for  a  circular  obstacle.^  To  compare  cne  above  re¬ 
sults  to  the  case  of  the  moon  introduce  the  quantity  —  where  d  is 
the  diameter  of  the  obstacle  and  2  is  the  observation  distance.  The 
moon's  average  diameter  is  1733.3km  and  its  average  distance  from  the 
earth  is  3.S44  x  10-km.^ 


1738.3km 
3.844  x  105km 


=  4.522  x  10 


(A2-i; 


Compare  the  above  to  the  case  of  the  1/4  inch  obstacle  observed  at 
50  cm. 


,625cm 

50cm 


.0125 


(A2-2) 


As  x  decreases  the  difference  in  intensities  becomes  less.  Since  x 
for  the  moon  is  almost  three  times  less  than  the  result  in  A2-2  the 
intensity  of  the  spot  produced  by  the  moon  should  be  between  one  and 
one  half  the  intensity  produced  by  an  equal  sized  circular  obstacle. 
(This  assumes  the  results  by  Raman  and  Krishnan  apply  to  any  sized 
obstacle.)  The  girth  of  the  moon  should  have  little  effect  on  the 
observed  spot  intensity. 

Before  reaching  the  surface  of  the  earth  the  rays  diffracted  by 
the  lunar  edges  must  pass  through  the  earth's  atmosphere.  Variations 
in  atmospheric  density  will  introduce  phase  differences.  The  cumula¬ 
tive  effect  of  differences  may  result  in  destructive  interference 
which  diminishes  spot  intensity.  A  model  to  test  atmospheric  effects 


is  depicted  below.  The  flame  introduces  air  currents  which  cause 
variations  in  air  density  with  resultant  phase  differences. 


Figure  A2. 1 

Model  to  Test  Atmospheric  Effects  in  Spot 


The  size  of  the  spot  will  be  affected  by  tie  diameter  of  the 
obstacle  ar.d  the  observation  distance.  Taking  the  radius  of  the 
as  the  radial  distance  to  the  first  intensity  minimum  in  the  ring 
structure,  then  for  the  lunar  case, 
kr  d 


2 

2z 


2.4048 


spo 


_  4 . 8096z 

r2  kd  (A2- 

If  one  takes  ■  =  5000A  as  an  average  optical  wavelength  then 

=  8.46  x  10 

The  ratio  of  the  spot  size  to  the  shadow  size  is  on  the  order  of  5  x 

10"  1 :. 


This  makes  the  spot  impossibly  small  to  see  unless  the 
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observation  is  made  where  the  center  of  the  shadow  hits  the  eartn's  I 

surface  with  some  sort  of  magnifying  device. 

The  spot  should  exist  at  the  center  of  the  moon's  shadow  but  the 
effect  will  be  significant  only  for  points  in  space  a  much  greater  ! 

distance  behind  the  moon  than  the  earth. 


i 


i 


i 


Appendix  3 


ON  AXIS  INTENSITY  IMMEDIATELY  BEHIND  AN  ANNULAR  APERTURE 


The  maximum  intensity  on  axis  behind  an  annular  aperture  is  found 
from  equation  III-37.  The  argument  of  the  cosine  factor  for  maxima  is 
of  the  form  (2n  +  1)^  where  n  is  an  integer.  Then  the  intensity  is: 


I 

max 


(0,  0, 


z  ) 
max 


=  A2 


|_^1+  d2/4z2 


+ 


^1+  d2/4z2 


(A3-1) 


with  the  restriction  that  Z  satisfies: 


Zmax  4A(2n+l)  n  -  0,  I,  2  .  .  .  (A3-2) 

Expression  A3— 1  forms  an  upper  limit  to  the  oscillations  of  the  inten¬ 
sity  on  axis. 

Likewise  the  minimum  intensities  on  axis  satisfy  the  relation: 


I  .  (0,  0,  z  .  )  =  A2 

min  mm  i  i—  —y  r 

LV1+ 17  V 


I 


1+ 


_d£ 

4z2 


(A3-3) 


with  the  restriction  for  Z: 
D2-d2 


z  . 


min  8mA 


m  =  1 ,  2 ,  3,  .  . 


In  the  region  where  -  <<  1  and  -  <<  1  then  the  upper  and 

Lo2  Lz- 


lower  bounds  to  the  oscillatory  envelope  are  4IQ  and  0. 

Consider  the  region  immediately  behind  the  aperture.  Here  — — 

r,2  k*2 

and  — —  are  not  negligible.  For  very  small  z  these  factors  iecome 
Uz" 

very  large  —  becoming  the  dominant  terms  in  the  square  roois  in  the 
denominators  of  A3-1  and  A3-3.  The  upper  and  lower  limits  cf  the 
intensity  oscillation  envelope  both  go  to  zero  as  Z  gets  infinitesi¬ 
mally  small. 

If  the  obscuration  ratio  of  the  aperture  is  very  small  'D  >>  d) 
then  there  is  a  region  along  the  optical  axis  where  the  effect  from 

the  outer  edge  is  negligible  compared  to  that  from  the  inner  edge. 

d2 

Here  —  >>  1  but  —  is  not.  Then  A3-1  and  A3-3  become 


I  „  A2  - - - 

max, min  -  o 

1+  — 

4z2  (A3- 

which  is  the  result  for  the  circular  obscuration  alone. 

Physically  this  result  is  easily  explained.  As  depicted  below 
the  cosine  obliquity  factor  for  outer  edge  diffracted  rays  is  much 


smaller  than  that  for  inner  edge  diffracted  rays.  The  amplitudes  of 


Figure  A3J 

Ray  Diffracted  to  Hear  On-axis  Point 


The  following  two  graphs  show  the  on  axis  intensities  for  the 
region  immediately  behind  annular  apertures  with  obscuration  ratios  of 


1/100  and  1/1000. 


ON  AXIS  INTENSITY 


Graph  A3. 1 


ON  AXIS  INTENSITY 


AXIAL  DISTANCE  (METERS) 

OBSTACLE  DIAMETER  =  .01  M.  APERTURE  DIAMETER  =  10  N. 


Graph  A3.2 


The  graphs  demonstrate  that  on  axis  immediately  behind  an  annular 
aperture  (of  sufficiently  small  e)  the  intensity  pattern  is  character¬ 
istic  of  a  circular  obstacle  alone.  As  z  increases,  the  intensity 
oscillation  envelope  unfolds  and  the  pattern  becomes  characteristic  of 
the  annular  (or  circular)  aperture. 

The  off-axis  intensity  plots  for  an  annular  aperture  of  obscura¬ 
tion  ratio  1/1000  are  given  before  for  optical  axis  distances  of  .5m 
and  20m.  This  dramatically  shows  the  circular  obscuration  behavior 
for  small  z  and  the  circular  aperture  behavior  for  large  z. 


A _  ,,  RADIAL  DISTANCE  (.661  M> 

OBST.  DIAH.r  .010  H,  APES.  DIAM.  r  10.00  M.  2  BIST.  =  0.50  M. 


Graph  AJ.J 

Intensity  Distribution  for  Small  z 


Appendix  4 


MH4ICING  OF  UNABERRATED  BEHAVIOR  3Y  WAVES  CONTAINING 
DEFOC'JS  OR  SPHERICAL  ABERRATION 


One  of  the  important  results  from  Chapter  4  was  that  the  spacing 
between  maxima  and  minima  on  axis  behinu  an  annular  apertue  was  inde¬ 
pendent  of  any  aberration  present.  It  was  also  shown  that  defocus  and 
spherical  aberration  shift  the  positions  of  extrema.  Consider  the 
case  where  the  aberration  shifts  a  maximum  (call  it  the  m^  order 
maximum)  onto  the  position  of  the  ntn  order  maximum  for  the  unaber¬ 
rated  case.  Can  the  amount  of  aberration  present  be  determined? 

Assume  the  aberration  present  is  defocus  with  aberration  coeffi- 

cient  W  .  Then  if  the  position  of  the  mt:  order  maximum  for  aber- 
20 

rated  case  (z')  is  the  same  as  the  position  of  the  unaberrated  nth 
order  maximum  (z),  one  can  equate  equations  III-40  and  IV— 1 6 : 


P2(l-e2)  =  P2(l-e2) _ 

4A(2n+l)  4A (2m+l)  -  8W2Q(l-e2) 


(A4-1) 


4A(2m  +  1)  -  8W20(1  -  e2)  =  4A(2n  +  1) 

w  =  I*™!.  A 

(1-e2)  (A4-2) 

Both  m  and  n  are  integers  so  let  p  =  m  -  n  where  p  must  be  an  integer. 
Then  if  defocus  is  of  the  form: 
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the  extrema  for  the  aberrated  case  are  exactly  superimposed  on  the 
extrema  for  the  unaberrated  case.  (Equation  A4-3  is  also  the  expres¬ 
sion  obtained  by  considering  the  m^  order  aberrated  minimum  and  the 
n*'*1  order  unaberrated  minimum,  m  and  n  are  strictly  arbitrary. ) 

Thus  if  one  takes  measurements  of  axial  positions  of  extrema  and 
compares  them  to  the  unaberrated  positions  and  finds  them  identical, 
then  the  proper  conclusion  drawn  is  that  any  defocus  present  has  the 
form  given  by  equation  A4-3*  (Clearly  the  unaberrated  case  corre¬ 
sponds  to  p  =  0.)  In  this  case  the  amount  of  defocus  can  be  deter¬ 
mined  best  by  intensity  considerations. 

The  corresponding  condition  for  the  superposition  of  extrema  for 
spherical  aberration  is: 


W 


40 


P 

(1-e4) 


A 


(A4-4) 


ON-AXIS  INTENSITY  AT  Z 


SEHINO  A  CIRCULAR 


OBSCURATION  AND  ANNULAR  APERTURE 


On  axis  the  Rayleigh-Sommerfeld  integral  gi  :en  by  equation  IV-J 
is  reduced  to: 


13(0,  2) 


oo 

Ak  ikz  f  ik  f  tj—  +  4  — —  r  2 
=  —  e  I  e  2z  ,2  r  dr 

iz  J  L  d  -1  ii 


(A5-i; 


-r— —  ,  for  W  <0,  the  above  is: 
8W20  20 


-)  -  - 


8AkW2u  -ikd2/8W. 


r  4w2o 

l  a2 


+  — —  r 


x  Tjdrj 


8AkW20  -ikd2/8W 


20  /  ridri 


JLl  e-ikd  /8W2(j  2 

,2  1 


(A5-2! 


For  the  circular  obscuration  case  the  disturbance  function  (and  there 
fore  the  intensity)  takes  on  an  infinite  value  at  this  particular  z. 


If  the  upper  limit  is  not  infinite  then  (A5-2)  becomes: 


U(0, 


D2 


8W 


20 


4AkW20  -ikd2/8K  ( rl 


-  -  4  klJ  <1  -  t2)  e-ikd2/8W20 
1  20 


(A5-: 


The  intensity  is: 


^  -  8S“>" 


20 


U<0’  -  8§-> 


20 


=  A2k2W2Q2(l  -  e2) 


(A5- 


Thus  the  intensity  has  a  finite  value. 

An  alternate,  but  related,  question  is  can  an  on  axis  maximum 

n  2 

occur  at  z  =  - -r~ —  (thus  implying  infinite  intensity)?  For  this  to 

'"2  0 

happen  one  must  satisfy: 


z  = 


P2(l-e2) 


4A(2n+l)  -  8W2u(l-e2) 


8W 


20 


8W2QD2(1  -  e2)  =  -4AD2 (2n  +  1)  +  8W2QD2(1  -  c2)  (A5- 

or 

4AP2 (2n  +  1)  =  0  (A5- 


which  is  impossible.  Thus  the  intensity  should  remain  finite  in  the 
annular  aperture  case,  in  agreement  with  A5-4* 


a 


2B2 


Appendix  6 


SEPARATION  OF  DEFOCUS  AND  Si  HER I CAL  ABERRATION  BY  AXIAL 
SHIFTING  OF  MAXIMA  WITH  DIFFERENT  03SCURATI0N  RATIOS 


From  equation  IV-66  the  axial  shift  from  the  unaberrated  case  for 
the  nth  maximum  behind  an  annular  aperture  when  both  defocus  and 
spherical  aberration  are  present  is: 


Az 


=  d2(i-e2)2 


W  +  W  (1+e2) 
_ 2  0 _  4  0 


X  (2n+l)  X(2n+1)  -  2K  (1-e2)  -  2W  (l-e4)l 

L  20  40  J 

(IV-66) 


Az  X(2n  +  1)  - 


2W2  0(1  " 


C2) 


"  2W40(1 


D  (I-e2)2 
2  (2n+l) 
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+  W  (1  +  e2) 
40 


W 

20 


+  W  (1  +  e2)  =  2Az\2  (2n-bl)2  _  Az4X(2n-H) 
40  D2 (1-e2) 2  D2 (1-e2) 


W 

20 


Az4X (2n+l) 
D2 (1-e2) 2 


W 


20 


(1  +  t2) 


(W  +W  (1+e2))  P^(l-£2)  +  Az4X(2n+l) 

K  20  4 0 ' 1  £  )}  '  : - 

D^ ( 1-e2) 


2Az\2 (2n+l ) 2 
D2 (1 -c2) 


W  +  W  (1  +  e2) 

20  40 


_ 2Az'1  2  (2n+l)2 _ 

(D2  ;i-e.2)  Az4X(2n+l))  (1-c2) 


(A( 


If  the  axial  shifts  from  two  different  obscuration  ratios  are 
measured,  there  will  be  two  equations  with  which  to  find  the  two 

unknowns  (V  and  W  ).  Subtracting  the  two  resulting  equations  a 

2  0  4  0 

solving  for  W  yields: 

40 


w 

40 


2AZj  X2 (2n+l)2 

(D2 (1-Ej  2)  +  Az  4X(2n+l))(l-e12) 


2Az_,X2(2n+l)2 

(D2(l-e22)  +  Az24X(2n+l))  (l-c22) 
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1 ’  ' 


Az,(l-e„-)z  -  Az_(l-e_z)^|  +  4A(2n+l)Az  Az  (e 


(D2 (l-ei2)  +  Azi4A(2n+l))(D2(l-e22)  +  Az  4A(2n+l)) 


The  above  unwieldy  expression  has  no  readily  apparent  simplification 
Substituting  back  into  equation  Ao-1  yields: 

2A2 (2n+l) 2  AZj (l-c22)  (e12-c22) 
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Appendix  7 


COMPARISON  OF  CALCULATIONS  USING  SINGLE  PRECISION 
AND  DOUBLE  PRECISION  NUMBERS 


Most  of  the  graphs  in  this  thesis  were  calculated  and  plotted 
using  an  IBM  Personal  Computer.  The  IBM  can  store  numbers  with  either 
single  or  double  precision.  Using  single  precision,  only  six  digits 
are  accurate.  With  double  precision,  17  digits  are  accurate.^  The 
drawback  to  using  double  precision  numbers  is  the  longer  program 
execution  time  required.  For  example,  a  typical  run  to  produce  an  off 
axis  intensity  contour  plot  took  approximately  ten  minutes  with  single 
precision.  Double  precision  required  almost  twenty  minutes  for  the 
same  run.  Because  of  the  large  time  differences,  it  is  important  to 
ask  whether  the  loss  of  accuracy  using  single  precision  numbers  is 
significant. 

To  test  this,  four  runs  were  made.  The  program  calculating  the 
off-axis  intensity  distribution  takes  the  most  time  and  has  the  high¬ 
est  number  of  calculations  of  the  programs  used.  The  maximum  and 
minimum  values  for  the  intensity  distributions  produced  using  double 
and  single  precision  were  compared.  The  first  two  runs  used  the  case 
of  a  one  centimeter  circular  obscuration  viewed  at  a  distance  of  one 
meter.  The  last  two  runs  used  the  case  of  an  annular  aperture  of  outer 
diameter  one  centimeter  and  inner  diameter  .3352  centimeters  viewed  at 
the  position  of  an  on-axis  minimum  at  .5017  meters. 


*7  r  7  ^  r  7  "  *' 


'  J-w  f  ■  <r 


V  »."  V*’ 
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For  the  circular  obscuration  case,  22  points  corresponding  to  the 
radial  distances  for  maxima  and  minima  up  to  .7  millimeters  from 
center  were  compared.  The  radial  distances  were  calculated  to  the 
nearest  one  and  one  half  microns  and  the  corresponding  normalized 
intensities  to  the  nearest  .0001.  Of  the  22  radial  distances  only  two 
differed.  Likewise  only  two  intensities  differed  in  value.  These 
differences  are  listed  below. 

Table  A7.1 

Comparison  of  Single  Precision  Versus  Double 
Precision  for  Circular  Obscuration  Run 


s 

Extremum 

Type 

Single  Precision 
Calculated 
Position  (mm) 

Double  Precision 
Calculated 
Position  (mm) 

Difference 

(mm) 

Minimum  #4 

.237 

.236 

.001 

i 

Minimum  #5 

C\J 

O 

KN 

• 

o 

o 

-.002 

Extremum 
Type  . 

Single  Precision 
Calculated 
Intensity 

Double  Precision 
Calculated 
Intensity 

Difference 

1- 

4 

Maximum  #5 

.4076 

.0477 

.0001 

► 

Maximum  #5 

.0000 

.0001 

.0001 

4 

( 

I 

The  plots  for 

each  case  reveal  no 

apparent  difference. 

•  J 


r4 


h 


For  the  annular  aperture  case,  the  first  46  extrema  out  to  a 
radial  distance  of  1.1  millimeters  were  compared.  Of  the  46  positions 
10  differed  each  by  only  .001  millimeters.  Only  two  intensities 
varied.  The  differences  are  listed  below. 

Table  A7.2 

Comparison  of  Single  Precision  Versus  Double 
Precision  for  Annular  Aperture  Run 


Extremum 

Type 


Single  Precision 
Calculated 
Position  (mm) 


Double  Precision 
Calculated 
Position  (mm) 


Extremum 

Type 


Single  Precision 
Calculated 
Intensity 


Double  Precision 
Calculated 
Intensity 


Difference 

(mm) 


Difference 


Maximum  #3 
Maximum  £9 


.3156 

.0996 


.515 

.099 


.0004 

.0002 
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ed. 


Notes 


^  BASIC  IBM  Personal  Computer  Hardware  Reference  Library,  1st 
pp.  311. 
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Appendix  8 


COMPUTER  PROGRAM  LISTINGS 


The  programs  listed  below  were  used  to  study  on  and  off  axis 
intensity  patterns  behind  circular  obscurations  and  annular  apertures. 
The  programs  are  written  in  advanced  3ASIC  using  DOS  1.1  for  the  IBM 
Personal  Computer.  The  advanced  features  of  BASIC  allowed  for  graph¬ 
ics  with  a  screen  resolution  of  200  x  640  pixels. 

The  first  program  listed  is  KAXMI'J.  This  program  lists  the 
positions  for  the  maxima  and  minima  on  axis  behind  an  annular  aper¬ 
ture.  The  computer  prompts  (inputs)  are  self-explanatory.  Allowance 
is  made  for  the  presence  of  defocus,  spherical  aberration  and  a  dis¬ 
placed  obstacle  and  aperture.  The  distance  along  the  optical  axis  is 
measured  from  the  aperture.  If  Az  is  the  separation  between  the 
aperture  and  obstacle,  the  nth  order  maximum  satisfies: 


z 

max 


JL 

2 


_ D2  -  d2 _ 

4X(2n+l)  -  8W  (1-e2)  -  8W  (1-e4) 
2  0  *+0 


-  Az  + 


The  nth  order  minimum  has  the  same  form  with  (2n  +  1  )  replaced  by  2n. 

The  second  program  is  AXIS.03S.  This  program  plots  the  on  axis 
intensity  function  behind  a  circular  obscurration,  using  equations  IV' 
43  and  III-34. 

The  third  program  is  AXIS. AM.  This  program  is  the  analog  of 
AXIS.0B3  for  the  annular  aperture.  It  employs  equation  IV-72  for  r^  : 

0. 

The  fourth  program  is  ARAG0.03S.  This  provides  an  intensity 
contour  plot  and  a  radial  intensity  function  plot  for  the  diffraction 
pattern  at  a  given  distance  behind  a  circular  obstacle.  This  program 
uses  equations  IV-43  and  III-34. 

The  last  program  is  ARAGO.AM  which  yields  the  plots  mentioned 
above  for  an  annular  aperture.  Equation  IV-72  is  used.  This  program 
as  well  as  ARAGO.OBS,  includes  a  subroutine  to  find  Bessel  functions 
of  the  first  kind.  The  subroutine  comes  from  BASIC  Programs  for 
Scientists  and  Engineers  by  Killer. 1 


5  d  [ m  fiAx - 

10  ’  MAXIMA  AMD  MINI!'! A  OH  AXIS . UNIFORM  BEAM,  NO  ABERRATIONS 

20  re-.'  f 

40  PR  TAT 

30  PRINT  "  THIS  PROGRAM  PROVIDES  A  LISTING  Or  MAXIMA  AND  Ml  Mi  HA  ON  A  X 1 3 
60  PRINT  "BEHIND  A  TARGET  CONSISTING  OF  A  DISPLACED  OBSTACLE  AMD  APERTURE. 
70  PRINT  "ME  ASSUME  A  WAVELENGTH  OF  6320  ANGSTROMS. 

80  PRINT 

90  PRINT  "INPUTS  ARE  SEPARATION  OF  OBSTACLE  AND  APERTURE,  DIAMETERS  OF 
100  PRINT  "OBSTACLE  AND  APERTURE,  AND  ABERRATION  COEFFICIENTS. 

1 1 0  PR  I NT 

120  INPUT  "INPUT  THE  SEPARATION  OF  OBSTACLE  AND  APERTURE  T!l  METERS" 5  DELE 
130  PRINT 


140 

INPUT 

"INPUT 

THE  DIAMETER  OF  THE  APERTURE  IN  METERS" 5  DO 

1 50 

PRINT 

1 60 

I NPUT 

"  I  NPUT 

THE  DIAMETER  OF  THE  OBSTACLE  IN  METERS";!. 

165 

EPS  = 

D/DD 

1 70 

PRINT 

172 
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AMOUNT 

OF  DEFOCUS  IN  TERMS  OF  WAVELENGTHS” 5 

174 

PR  I  NT 

176 

I  NPUT 

"  INPUT 

AMOUNT 

OF  SPHERICAL  ABERRATION  IN  TERMS  Of 

1  78 

PR  1  NT 

ISO 

I  NPUT 

"DO  YOU  WISH 

A  PRINT  OUT  ON  SCREEN  <Y/N) "5 PRT* 

1 90 

IF  PIT 

T$  =  "Y 

"  THEN 

PRT  *  1  ELSE  PRT  =  0 

191 

1. 72 
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...CJLATC 

MAXIMA 
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=»  DFC*L 

AM 

V0 1 0  ■■■■■■■  3A*LAM 
700  CO EE  1  ■  DD*DD  ~  D*D 
710  C  0  E  F  2  D  D  *  0  D 

:  1 5  CCE F 7  -  2  *  M 0 2 0 ■  *  (  1  - E P S * E P S  >  +  2 * W 0 4 0 *  ( 1  -  •  E P S * E P 3  * E P S * E P  S ) 

220  FOR  N  «  0  TO  40 

222  C0EF5  *  LAM-*  < 2-*N *•  1  >  -  C0EF3 

224  C0EF6  ■  LAM*2*N  -  C0EF3 

250  MAX  (N  i-  l )  .  (C0EF1/  <4*C0EF5>  -  DELZ  +  SOP  (  (DELZ . C0EF1  /  (4*C0EF5)  )  A2  +  C0EF2* 

DELZ/C0EF5) ) /2 
755  IF  M  >  0  GOTO  270 

260  M I N  < N • 1  •  1  >  =  ( CHEF  1  /  ( 4*-C0EF6 )  -  DELZ  +  SC!R  (  (DELZ-C0EF1  /  C4*C0EF6)  )  ""'2  •+•  C0EF2* 
DELZ/C0EF6; ) /2 
770  IF  PRT  —  0  GOTO  300* 

200  r-’P  I  NT  N MAX  ( 1 1*  1 )  ,  M I N  1 N+  I.  > 

700  NEXT  M 


302  ’  PRINT  OUT  RESULTS 

,  y  )  •’ 

!.  0  '...PRINT  "MIX  M-!A  AND  MINIMA  ON  AXIS: 

32v  I..PPINT 

730  LPRIMT  "SFPA-:: 'vr  '  •  if  f  PF  OBSTACLE  AND  APERTURE  =  "  ;  DELZ :  "  M.  " 


365  LPRINT  "BEFOOLS  ~  "•'"""5"  WAVELENGTHS 
ELENGTHS" 

366  LPRINT 

370  LPRINT  "ORDER  ( N )  Z (MAXIMA) 

380  LPRINT  " -  - 

390  FOR  N  -  1  TO  -U 
400  LPRINT  USING  "  ## 

410  NEXT  N 
420  PRINT 

430  INPUT  "DO  YOU  WISH  TO  MAKE  ANOTHER  RUN 
440  IP  0*  =  "Y"  THEN  GOTO  10  ELBE  END 


20  KEY  OFF 
70  PRINT 

40  PRINT  ”  HUS  PROGRAM  CALCULATES  THE  INTENSITY  ON  AXIS  BEHIND  AM  CIRCULAR 
5  0  P  F:  I N  T  ' '  0 B  S T A C  L E 
60  PRINT 

70  PRINT  "HE  ASSUME  A  WAVELENGTH  OF  6328  ANGSTROMS. 

30  PRINT 
100  ’ 

101  ’  INPUT  PARAMETERS 

1  (  I  •* 

110  INPUT  "WHAT  IS  THE  MAXIMUM  DISTANCE  BEHIND  THE  SCREEN  TO  BE  CALC  'LATEO" 5  ZMAX 


140  PRINT 

145  INPUT  " NH AT  IS  THE  MINIMUM  DISTANCE  BEHIND  THE  SCREEN 


"□  BE  CALCULATED" 5 ZMIN 


146  PRINT 

147  Z1  -  2 MAX  -  ZMIN 

150  ’ 

151  ’  CALCULATE  THE  RANGE  ALONG  THE  OPTICAL  AXIS  “0  BE  PLOTTED 

l  '5  2  ' 

155  Zt  *  ZMAX  -  ZMIN 

160  INPUT  "WHAT  IS  THE  DIAMETER;  CF  THE  OBSTACLE  "METERS*  ”  5  D 
170  PRINT 

190  PRINT:  INPUT  "WHAT  IS  THE  CCEFF  1 3 1  FT  IT  FOR  DEFOCUS  IN  TERMS  CF  WAVELENGTH "  5  Wi 
195  W020  *  6. 32SE--07  *  W1 

200  PRINT:  INPUT  "WHAT  IS  THE  COEFFICIENT  FOR  SPHERICAL  ABERRATION  IN  TERMS  OF  W 

AVELENGTH" : W2 

210  WO 40  ■  6 . 3 2 BE -07  *  W2 

220  PRINT:  INPUT  "WHAT  IS  THE  GAUSSIAN  CONSTANT" ;GC 

230  ’ 

231  ’  CALCULATE  THE  GAUSSIAN  FACTOR 

*--»  -T  ? 

240  GF1  »  EXP  <— GC/4) 

(390  •’ 

891  ’  TRANSFER  TO  COLOR  DISPLAY 
392  ’ 

900  PRINT:  INPUT  "MAKE  SURE  VIDEO  SCREEN  IS  ON  AND  HIT  RETURN  WHEN  READY'  ARB* 
''•HO  DEF  SEG  ~  0 

?  2  0  P  0 1 <  E  0  H  4 1 0  ,  ( P  E  E  K  <  &.  H  4 1 0 )  A  N  D  0  H  C  F  *  0  R  h  H 1 0 

930  SCREEN  1,0, 0,0  :  SCREEN  0:  WIDTH  80 

940  LOCATE  ,,,1,6/7  DEF  SE3 

950  GOTO  1000 

960  ’ 

L  ’  mANT'EF  BACK  70  i  lO'C.IH  .ViS 

" 

970  OEF  SEG  •=  0 

HSO  POKE  CM 4  1 0  ,  (PEEK  ( LH4 1 0  >  OR  TH30> 

■-9.>  CCPIIEM  o:  WIDTH  40:  WIDTH  SO:  LOCATE  .  ,1,12,  13:  DEF  SEG:  GOTO  :  770 

1  •  )i  >o  ’ 

'  •  ‘  1  L  ■  F’!...OT  "HE  AXES 


1010  SCREEN  2 
’  020  'T0E  I  •-  l  TO  20 
1030  LOCATE  I,  12 


1  <.>40 
1 050 
1060 
1 070 
1030 
1090 
1  1 00 
1  1  .1 0 
1  120 
1121 
1122 
1  125 
1  130 
1140 
1150 
1160 
1  1  70 
1  1  GO 

4  1  *?'•  * 
1 200 
1 2 1 0 


1240 


PRINT  CHRT <  1 7°) ; 
NEXT  1 

LOCATE  21,12 
PRINT  CNR  1  ( 192: 
FOR  I  =  1  TO  66 
LOCATE  21,121 
PRINT  CHR$ ( 196) ; 


LABEL  Y  AXIS 


1  300 
]  3 1 0 
1  33  0 
1  330 
;  2 

1350 


LOCATE  16,9 
PRINT  1 ; CHR$ (197) 

LOCATE  11,9 
PRINT  2;;  CHRft  (197) 

LOCATE  6,9 
PRINT  7-  C!  !R$  ( 197) 

LOCATE  1 , 9 
PRINT  >  i  Cl  iE*  (197) 

FOR  I  -  1  TO  9 
LOCATE  6+1,3 

ON  I  GOTO  1 230 , 1 250 , 1 270 , 1 290 , 1 250 , 1 3 1 0 „ i 230 , 1 270 , 1 330 
PRINT  "I" 

GOTO  1340 
PRINT  "N " 

GOTO  1340 
PRINT  "T" 

GOTO  i 340 
PRINT  "E" 

GOTO  1340 
PRINT  "S" 

GOTO  1340 
PRINT  "Y" 

NEX  r  I 

CALCULATION  AMD  PLOTTING  OF  THE  INTENSITY  FUNCTION  USING  EQUATION 


1355  Z  »  Z 1/540 
1 360  FOR  I  -  1  TO  540 

1405  ABE  -  <  1-16*W040*  (ZMIN  +  2  +  1 )  /  (D*D>  )  /  ( 1+8-SW020*  (ZMIN+Z*!  ) 

1410  X  -  93  +  I 

1420  V  =  164  -  Cl  NT  ' AB2*AE»2*GF  1  *GF1  ■“■40) 

1425  IF  3:  =  1  GOTO  144  0 

L 4 30  LINE  (  X  ,  XO,  YO) 


( D-*D)  ) 


Y  0  " 

NEXT  i 

LABEL  ,<  AXIS 


!.  590 


pon  i nt (log  e: i :  leg  :  iO  ’) 

SCALE.  -  Zl/1. 0'T'OVj 

IF  SCALE  :  4  i"l  ip; 1  EACr  -  2  EL 

(  1MIT  n  :  <.,,(•  •  /  c-  r  . ■  r-  «.  •  p'T ' 

FOR  I  -■  1  iO  IN!  S'"  iLE  -  PACT 
LOCA""E  71,12  +  I+UNir 
PR  IN  I"  OHR‘5  <  177) 


1600  LOCATE  22,7  I*UNIT 
i  j£j  i  o  p  pynjT  Z  !i  l ! !  1  ■*  i  o  p  n  y  /  F  ^  C  T 
1620  NEXT  I 
1630  LOCATE  27,70 

1640  PRINT  "AXIAL  DISTANCE  (METERS! " 
1.650  LOCATE  7  1,5 

1660  PRINT  "DIAMETER  OF  2 OCCUPATION  =  " 
1670  LOCATE  25,5 

1  l  MT  "'  j  io  -  "  s  ’■■•Z:  "LArlRDA  E-AU3S 

1730  L.CCATE  1,1 
1740  ’ 

1.741  ’  PAUCE  TO  ALLOW  T  PINT:  NO  CF  PLOT 
1742  ’ 

1745  "Oil  I  '  TO  3000 

1750  NEXT  I 

1  7.60  GOTO  9”  > 

1  770  71...  S 

17S0  INPUT  "DO  YOU  WISH  TO  MAKE  ANOTHER 
1790  IF  D*  =  "N”  THEN  END  ELSE  GOTO  1.30 


THE  INTENSITY  ON  AXIS  BEHIND  AN  A!  'NULAR 


V  T  ~  V  Ji'  I 


10  CL 3 
20  KEY  OFF 
30  PRINT 

40  PRINT  "THIS  PROGRAM  CALCULATES 
TO  P R IN"  " AF'ERTU F: E 


*;j  O  f 

='R  i  m  r 

70  PRINT 

HO  1 

~’R 1  NT 

1 00 

1 0  :i. 

7  IN! 

102 

1  1 0 

I NPUT 

•  40 

PRINT 

145 

INPUT 

1.46 

PRINT 

147 

Z1  = 

150 

* 

151 

•'  CAI 

152 

.*  55 

Z 1  ~= 

1 60 

INPUT 

1  70 

PRINT 

i  90 

PR  I  NT 

75 

WO  20 

200 

PRINT 

,v,;r! 

LENGTH 

2 1 0 

WO  40 

220 

PR  I  NT 

2  30 
231 

7  CAL! 

-r  •**> 

7!  40 

(3! 

.760 

0! 

■490 

■191 

’  TRO1 

692 

1  •  t )  *  * 

• 1 0 

PR  I  NT 
DEF  S! 

■?20 

:>7>0 

PC):  E  : 
SCREE 

940 

L.QCA  t: 

'’SO 

'  2) '  ' 

GGTG 

'  S  I. 
r  "■ /.  "" 

TF'.'ii 

;,,j  j 

pi  E 

;On 

,  T r-rr 

ZMAX  -  Z M I N 


SEb 


:  U '  DTM  4 :  HjnrH  !?<"•:  I  ..OF  ATE  ,.1,12,13:  DEF  SE6:  GOTO  1770 
OT  THE  AXES 


r 


1020  FOR  I  *  1  TO  20 
1030  LOCATE  1,12 
1 0 4 0  P PINT  C H R  *1(173)  • 

1050  NEXT  I 

1060  LOCATE  21.12 

1070  PRINT  CNF  £  192) 

1030  FOR  I  =  1  TO  66 
1 090  LOCATE  2 1 ,  1 2  + 1 
1 1 00  PR  I  NT  CHR  p  •:  1  96 )  ; 

1110  NEXT  I 
1 1 20  ’ 

1121  •’  LABEL  V  AXIS 

1122  ’ 

1125  LOCATE  16,9 
1 1 30  PR  I  NT  1 CHR$  (197) 

1140  LOCATE  11,7 
1 1 50  PR  It  iT  2 .!  CURT  <  1 97 ) 

1 1 60  LOCATE  6,9 
1 1 70  PR  I NT  3 5 CHRf ( 1 97 ) 
t 130  LOCATE  1,9 
1 1 90  PR I NT  4 ; OHP$ ( 1 97 ) 

1 200  GT>  =  "  NORMAL  I  Z  E  "  :  HT>  «  "  I N7GM5 I  TV  " 

3.210  FDR  1  =  1  TO  9 
1220  LOCATE  6+1,3 
1 230  PR I NT  M 1 04  < 6$ ,1,1); 

'710  LOCATE  6+1,6 

1250  PRINT  MID$’(H«,  I,  1)  5 

1.260  NEXT  I 

1270  LOCATE  16,3 

1 230  PRINT  " D " 5 

1350  ’ 

1331  ’  CALCULATION  AND  PLOTTING  OF  THE  INTENSITY  FUNCTION  USING  EQUATION  A 10-4 
1352  ’  KK  =  WAVE  NUMBER 

: 355  Z  -  Z 1/540 
1360  FOR  I  =  1  TO  540 

1370  C0EF4  -  4*Z*Z*GC*GC/  ( DD *DD*DD*DD*KK*KK ) 

1390  C  -  COS  ( 1241147  !  *  (  <DD*DD  -  D*D)  /  < ZMIN+Z*  I  >  +  B*W020*  <  1-D*D/  (DD*DD)  >  + 

32*W040* ( 1-D*D*D*D/ (DD*DD*DD*DD) ) ) ) 

1O0  A3 1  «  <1-1 6*W040*  <  ZM I N+Z* I )  /  (DD*DD)  )  /  <  (1+8*W020*  < ZMIN+Z*. I )  /  (DD*DD>  )  *SQR 

1  i-DO+DD/  <  4+  (ZMIN+Z* I )  *  (  ZMIN+Z* I )  )  )  ) 

1  105  AB2  -  <  1  —  1 6 *14040*  <  ZMIN+Z*  I )  *D*D/  <DD*DD*DD*DD)  )  /  (  <1+8*W020*  <  ZMIN+Z*I)  *D* 

rV  ( D D D D * D D * D D ) ) *SQR ( 1 +D*D/ (4* (ZMIN+Z* I ) * ( ZMIN+Z* I ) ) ) ) 

’  ■’’•6  00EF1  =  ABL  -*AB1  *GF1*GF1 

!  107  C0EF2  »  AB2*AB2*(3F2*GF2 

!.  C0EF7  --  9  »GF  1  *BF2  * AB 1  * AB2*C 

Y  rr  93  ••  T 

' 420  Y  -  164  -  Cl NT ( (COEF 1  +  C0EF2  -  C0EF3) / < 1+C0EF4) *40) 

'  *  73  IF  I  =  1  GOTO  1440 

'•  11  "0  LINE  <X, (X0,Y0> 


VO  =  y 


LABEL 


1550 

IF  SEAL 

E  .  1  Ti  iEH  FACT  =  2  : 

!  560 

UNIT  = 

66/1  SCALE  -*  fact: 

!.  570 

~CF:  f  - 

i  '  jT  ^  L  j  ” 1  j  *"  1  L  ^  1  | 

1 530 

LOCATE 

2  U  1 2  I  ••UNIT 

;.  590 

F:  T  N  r  .. 

*  ■»-* 

1  oO*,‘ 

~2 ,  iO  ■  1  ” 

1  6  1. 0 

?r  J  :  |  r  7 

•i  r ;  | ;  ;i  PC: L  /F  AC  r 

;  i,  ’«  ') 

l  G  Tr 

1 1 , 30 

;  64('> 

PRINT  " 

a  v :  al  r; :  ::t  -.nce  •  metei 

i.  650 

LOCATE 

■>  rrr  — ? 

-  -J  H 

1660 

PRINT  U 

3 1  MO  "CBS.  DIAL!.  ~  # 

j"r.  .l,  „  DD.GC; 

!  380  'PRINT  "  WO. 20  -  "  3  W1 5  "LAMBDA  WO  40  "?W2;  "LAMBDA  GAUSSIAN  CONG  TAN 

r  ==  ">GC; 

1  730  LOCATE  1  ,  J. 

L  740  ’ 

1741  1  PAUSE  TO  ALLOW  PRINTING  OF  PLOT 

*742 

1745  FOR  I  -=  1  TO  3000 
1750  NEXT  I 
i.  760  GOTO  970 

i  ’7  70  P’  P 

1.  730  INPUT  "DO  YOU  WISH  TO  STORE  THE  PLOT  ON  DISK  (Y/N)  "  3  A* 

1 790  IF  Al>  "M"  GOTO  1900 
1  300  ' 

1001  i  ROUTINE  TO  STORE  PLOT 


t  SOS  7  '•  I N  r 

I  .!7UT  "WHAT  IS  THE  NAME  YOU  WISH  TO  STORE  THE  PLOT  UNDER"  "FIT 
•  3  2  o  F  7 !  N  A  i '!  %  ' '  B ;  "+  F 1  * 

’  b:zf  see  «  s-HBaoo 

40  PSAVE  F I  LNAH$ ,  0 ,  ?<H4000 

‘900  >  NPIJ  r  "DO  YOU  WISH  TO  MAKE  ANOTHER  RUN  < Y/N)  "  3  D* 
c- !•  >  IF  m  "N"  THEN  END  ELSE  GOTO  130 


1 


PROGRAM :  ARAGO . CDS 


10  I  'EY  OFF 
15  DIM  IMT2N<  ’  *0 
20  CL 3 
GO  PFTU 

40  PR  I i  !T  "ARAGO  PLOTTING  ROUTINE: 

50  PRINT 

60  P~,T!!T  "THIS  PROGRAM  WILL  GIVE  A  CONTOUR  PLOT  OF  THE  SPOT  OF  ARAGO  AT  A 
70  PRINT  " CHOGEN  DISTANCE  BEHIND  AN  OPAQUE  CIRCULAR  TARGET. 

30  PRINT 

•?0  PRIM!’  "THIS  PROGRAM  ASSUMES  A  WAVELENGTH  OF  631...  ANGSTROMS.  INPUTS  ARE 

100  PRINT  "THE  OBSTACLE  DIAMETER.,  THE  TOTAL  RADIAL  DISTANCE  OF F  THE 

*10  C'DT‘F  "OPTICAL  AXIS  FOR  COMPUTATION,  THE  ABERRATION  COEFFICIENTS  FOR  DEFOCUS 

ICO  PRINT  "AMD  SPHERICAL  ABERRATION,  THE  GAUSSIAN  CONSTANT  AND  THE  DISTANCES  Z1 
125  PRINT  "AMD  Z2  (MEASURED  FROM  THE  APERTURE  AMD  OBSTACLE  RESPECTIVELY) . 

! 30  PRINT 
140  PRINT 

150  INPUT  "INPUT  THE  TOTAL  RADIAL  DISTANCE  TO  BE  CALCULATED  (IN  METERS)" 5 R1 
160  PRINT 

174  INPUT  "INPUT  THE  OBSTACLE  DIAMETER" ;D 
190  PRINT 

1 95  INPUT  "INPUT  THE  DISTANCE  Z2  FROM  OBSTACLE  TO  VIEWING  PLANE  (IN  METERS)" 5 Z2 
200  PRINT 

210  INPUT  "INPUT  THE  COEFFICIENT  OF  DEFOCUS  IN  TERMS  OF  WAVELENGTH " j  W 1 
220  WO 20  «  6.328E-07  *  W1 
230  PRINT 

240  INPUT  "INPUT  THE  COEFFICIENT  OF  SPHERICAL  ABERRATION  IN  TERMS  OF  WAVELENGTH" 
■  W2 

.'2  WO  :K»  «  6.32SE . 07  *  W2 

r"  r  T  MT 

270  I -.'PUT  "INPUT  THE  GAUSSIAN  CONSTANT" ;  GC 
280  PRINT 

29'.'  INPUT  "DO  YOU  WANT  TO  SAVE  THE  CONTOUR  PLOT  ON  DISK  <  Y/N)  "  *  CPLOTS 
'00  IF  CPLOT*  «  "N"  GOTO  330 

310  PRINT  :  INPUT  "WHAT  IS  THE  NAME  OF  THE  FILE  TO  STORE  THE  PLOT  UNDER"; FI  $ 

320  FU_1T  «  "B:  FIT 
330  PRINT 

740  INPUT  "DO  YOU  WANT  TO  SAVE  THE  INTENSITY  PLOT  ON  DISK  (Y/N) "5 IPLOT* 

350  : F  I PLOT %  «  "N"  GOTO  410 

360  PRINT:  INPUT  "WHAT  13  THE  NAME  OF  THE  FILE  TO  STORE  THE  PLOT  UNDER"  ?  F2iS 
370  F  I!  2-5  ~  "p;  ”  j- 
730  ON  ERROR  GOTO  390 
TO  RESUME  !  TXT 

••’■'O  PRINT:  INPUT  "DO  YOU  WISH  TO  BYPASS  THE  CONTOUR  PLOT  (Y/N)  " ;  PBY$ 

101  * 

302  ’  CALCULATE  ABERRATION  AND  GAUSSIAN  FACTORS 


•  1 0  P I  -  3 .  1. 4 1. 59 
1  0  A D 2  -  l  /  1 1  i- 8 * W 0 2 0 *  Z2/  (D*D )  ) 
'!20  A Ei 4  a  1-1 6-*W040*Z2/  ( D*D ) 

■2i  or::  «  cxr  :-gc/4> 

‘25  C0EF2  =  AB4*GF2*GF2/ (AB2*AE<2) 
429  10  *  C0EF2 

■‘30  PRINT  :  PRINT  "COMPUTING  ..." 
I  o  o  D?  -  O 


531  ’  CALCULATE  ARGUMENT  FOR  BESSEL 
FUNCTION 


■UNCTION  AMD  EVALUATE  RADIAL  INTNEIT 


535  X  -  (D  -  r  L  *  PI  ■  ' 
540  GO CUD  00 'DO 
550  INTEN(I)  -  C0EF2 * J"7-* J7 
560  NEXT  T 
570  CL 3 


091  7  TRANSFER  TO  COLOR  DISPLAY 
392  ’ 

900  PRINT:  [N'-'UT  "  MAKE  CURE  VIDEO  EC  REE. -I  IS  CN  AMP  MIT  RETURN  WHEN  REV"/"" 
°10  DEF  PEG  »  0 

920  POKE  5.H410,  \  PEEK  ( &H4  1 0 )  AND  OHCF'  OR  -vH  1 0 
930  SCREEN  1,0, 0,0  :  SCREEN  0:  WIDTH  SO 
940  LOCATE  ,  ,1,6,7  :  DEF  '3 EG 
930  GOTO  9°5 

960  ' 

961  ’  TRANSFER  SACK  TO  MONOCHROME 

962  ’ 

O70  DEF  SEG  •=  0 

"’90  POKE  5!  1410,  (PEEK  (S<H410>  OR  &H30> 

990  SCREEN  0*.  WIDTH  40:  WIDTH  80:  LOCATE  ,  ,  1,  12,  13:  DEF  SEG:  GOTO  1950 
995  SCREEN  2:  WIDTH  80 
tOOQ  IF  PBY*  -  "Y"  THEN  GOTO  1200 
1001  •’ 

1002  ’  CONTOUR  PLOTTING  ROUTINE,  FIND  INTENSITY  AT  RADIAL  DISTANCE  FROM 

SCREEN  CENTER  ( 1 00 , 320 ) . 

1003  ’ 

1004  FOR  [  =  l  TO  00 

1005  FOR  J  =  1  TO  24 
1007  J1  ■■=  ,J  + 1.  ,7 

1  o  t  o  R2  -  C I  NT  ( SQR  (  <  1 70- J  1  *3 )  *  <  1 70- J 1  *8 )  ( 320- 1  *8 )  *  ( 320- 1  *8 )  )  ) 

1032  IF  R2  <  1  THEN  R2  «  1 
1015  COSUB  2030 

1.020  ON  DUN  GOTO  1030,  1050,  1070,  1090,  1110 
1030  A*  -CHRf(O) 

1040  GOTO  1150 
1050  AT  -  CHR5 ( 176) 

1060  GOTO  1150 
1070  AT  =  CH.RT < 177) 

1030  GOTO  1150 
1090  A3  •-  CURS  U78) 

1100  GOTO  1150 
1110  A3  -  CHF:3  <  219) 

1150  LOCATE  J, I 
1160  PRINT  AT? 

1170  NEXT  ,7 
.1.175  NEXT  I 
1 1 77  GOSUD  2000 


t.  J  00  ' 

1131  ’  PAUSE  FOR  PRINTING  AMD  STORING  OF  PLOT. 

1 132  ’ 

1183  FOR  K  -  1  TO  5000:  NEXT  K 
1185  IF  CPLOTT  ==  "N"  GOTO  1200 
1190  DEF  SEG -  &H3800 
1 1 95  BSAVE  F I L 1 T ,  0 , &H4000 


DRAW  AXES 


1200  ' 
1201  ’ 
1 202  ’ 


1 205 

CLS 

1 2 1  o 

FOR  I 

;  i-n 

1 220 

LOCATE 

I  .  12 

1  230 

PRINT  r’ 

'  T'+  ■  1 

1240 

NEXT  I 

1250 

LOCATE 

21,12 

i  ^  c. 

PRINT 

■  i  ~  t-  ,■  * 

•  r'7i') 

FOR  I  == 

1  TO 

1.230 

LOCATE 

i  i 

1 290 

PR  I NT  ' 

HR*  •.  1 

1 300 

y 

'  30  L 

lab: 

v  ) 

1305  NEXT  I 

1310  LOCATE  16,7 

1 320  PR  I  NT  ..  25 1  CURT.  «'  IP'!') 

1330  LOCATE  11,3 
1 3 4 0  P  R INT  .51 C MR  * ; 1 9 " > 

1350  LOCATE  6,7 
1  36  0  PR  I  NT  .75;  Cl  1RT>  <  1  ^7  > 

! 770  LOCATE  1 , 9 

1330  PRINT  1  ;  CMR *  (  i. 97 > 

1 7?0  FOR  I  ■=  1  TO  9 

:  400  LABI*  -=  "NORMALIZE” ;  LAB2*  =  "INTENSITY" 

1410  FOR  I  -  1  TO  9 

14  2')  LOCATE  6+1,3 

!  4  PR  I  NT  MID*  ( LAB  1  *  ,  1 ,  1 )  1 

1.4  40  LOCATE  6+1,6 

■  iso  PR  I NT  M I D* ( LAB2* ,1,1)5 

■  4cO  NEXT  I 

i 4 70  LOCATE  16, 3 
1430  PRINT  »D"; 
i  530  ’ 


1531 

LABEL  X 

AXIS 

1  cr  -■ 

!• 

1  540 

LOCATE  24, 

25 

1  350 

FRINT  " RAD 

IAL  DISTANCE  (.001  M ) " 5 

1560 

1 56 1 

’  PLOT  RA 

D I AL  I NTENS I TY  FUNCT I ON 

156.2 

" 

i  565 

ODD  •=  1 250 

*  p;  1 

1 570 

ROW  -  t  1  ;t  ■ 

LOO 'ODD) /LOG (10) ) 

1  530 

SCALE  -  DO 

o/io  -paw 

i  r5c;'0 

IF  SCALE 

4  THEN  FACT  »  2  ELSE  FACT 

i.  600 

UNIT  -  .50' 

•SCALE  *  FACT ) 

1  6  1 0 

FOR  I  =  1 

TO  INT (SCALE  *  FACT) 

1620 

LOCATE  21, 

12  -  I- UNIT 

1  630 

PRINT  Cl-'R* 

1 1  97  > 

1.  Cl '  \  U 

LOCATE  22, 

!0  ~  I *UN I T 

1  650 

PRINT  1*10 

POW/FACT 

1  660 

NEXT  I 

1 665 

YO  =  164  - 

Cl NT ( 10*1 60) 

l  666 

XO  =  93 

V. 

V 


1 670 
1680 
1635 
1695 
1 700 
1705 
1 7 1 0 
1720 
1730 

If  • 

1735 

1740 

1745 

1750 

1760 

1761 

1762 
1765 
1770 
1S00 
1810 
1 320 
1830 
1840 


i860 
1 365 
1869 
i  370 
1330 
1890 
1 900 
1910 
1 920 
1 930 
1940 
1  > 
1 960 
1970 
2000 
200 1 
2002 
2005 
201 0 


Z  D1ST .  =  ##„###  M. 


FOR  I  =  1  TO  440 
Y  =  1 64  -  CINT ( INTEN ( I )  *  160) 

X  3s  93  +•  t 

LINE  -  X,  Y)  --  CXO, Y0> 

XO  =  X 
YO  =  Y 
NEXT  I 
LOCATE  25,2 

PRINT  USING  "  OBSTACLE  DIAMETER  =  „ #####  M.  Z  I)  I  ST.  *  ##.###  M. 

;  d,  zc:  :1 

LOCATE  1 , 1 

FOR  I  =■•  1  TO  5000 

NEXT  I 

IF  I»LQT*  -  "M"  GOTO  1300 

STORE  INTENSITY  PLOT  ON  DISK 

DEF  3ES-=  &H38Q0 
BSAVE  F I  L2$ .  0 ,  5414000 

PRINT:  INPUT  "DO  YOU  WISH  A  LISTING  OF  THE  INTENSITIES (Y/N) ";D1* 

IF  Di$  -  "Y"  THEN  GOTO  1820  ELSE  GOTO  1940 

LPF:  I  NT :  LPR  t  NT  .*  LPR I  NT :  LPR I  NT :  LPR I  NT :  LPR I  NT :  LPR I  NT :  LPR I  NT  :  LPR I  NT 
LPRINT  :  LPR  I  NT  :  LPR  I  NT  :  LPR  I  NT  .*  LPR  I  NT  :  LPR  I  NT  :  LPR  I  NT 
LPR I NT  "RADIAL  DISTANCE (MM)  INTENSITY  EXTREMUM  TYPE 

IF  10  -  INTEN  (1)  >  0  THEN  EXT*  =  "MAXIMUM"  ELSE  EXT*  «  "MINIMUM" 

LPR I NT  US  I NG  "  # . ###  # . ### " ; 0 , 1 0  5 

LPRINT  "  "JEXT* 

S I GN  =  SGN ( I 0- I NTEN ( 1 ) ) 

FOR  I  :=  1  TO  439 

IF  SQN(INTEN(I+1) -INTEN <I) )  =  SIGN  GOTO  1930 

IF  SIGN  >  0  THEN  EXT*  «  "MAXIMUM"  ELSE  EXT*  =  "MINIMUM" 

S I  GN  =■  SGN  ( I  NTEN  ( I  + 1 )  - 1  NTEN  ( I )  ) 

LPR I NT  US  I NG  "  # . ###  # . #### " 5  I *R* 1 000 , I NTEN ( I ) ; 

LPRINT  "  ";EXT* 

NEXT  I 
GOTO  970 
CLS 

INPUT  "DO  YOU  WISH  TO  MAKE  ANOTHER  RUN  (Y/N) " ?  D* 

IF  D*  =  "N"  THEN  END  ELSE  GOTO  130 

LABEL  CONTOUR  PLOT 

LOCATE  24, 1 

PRINT  "  LEGEND:  " ” CHR* (221 ) ; "  I  <  .002  la  " ; CHR* (173) ; "  I  <  .01  la  " ; CHR* 
-  < "  I  <  .05  la  CHR* <176) 5 "  I  <  .2  I  a  " ‘ CHR* (0) ; "  I  >  .2  la"; 

9PTI  IQt\| 


FIND  SHADING  FOR  CONTOUR  PLOT  POINTS 


2035  IF  I NTEN ( R2 ) 
2040  DLJM  =  5 
2050  RETURN 
2060  IF  I NTEN  < R2 ' 
2070  DLJM  «  a 
2080  RETURN 
2090  IF  INTEN (R2) 
2100  DUM  ~  3 
2110  RETURN 


002  GOTO  2060 


01  GOTO  2090 


.05  GOTO  2120 


2 1 20 
2130 
2140 
2 1 50 
2 1 60 


IF  I NTEN ( R2 )  .2  GOTO  2150 

rjr-i  -  2 

RETURN 
r<  ir-t  --  1 
RETURN 


3001  ’  SUBROUTINE  TO  FIND  BESSEL.  FUNCTIONS  OF  THE  FIRST  KIND 

3002  ’ 

00 1 0  T1  -  9.99999919—06 
8020  PI  =  3.14159 
3030  X  2  "  X  -r  X 

0040  IF  (  \  :<-=0 '  AMD  (D9=l)  )  THEN  3240 

3050  IF  (X  .  15)  THEN  S260 

3060  IF  (D9  :>  0)  THEN  3090 

3070  S6  =  1 

9030  SOTO  Pi -0 

3090  X3  =  X 

n  1 00  X  -  D9  *•  L 

8110  GOSUB  8290 

••>120  X  -  X3 


3 1  30 

36 

-  ( X 

3140 

T3 

=  S6 

0 1 50 

17. 

=  0 

•  160 

IF 

ABSC 

8170 

17. 

--  17. 

1130 

T  4 

—  T"' 

3190 

T3 

=  —  T  ■ 

2200 

36 

=  36 

32 1 0 

GOTO  31, 

3220 

J  -7 

iJ  / 

•=  S6 

3230 

RETURN 

3240 

3  7 

(I) 

-'250  RETURN 

8260  ’  ASYMPTOTIC  EXPANSION 

8270  J7  »  SQR  (2/  <PI  *  X)  )  *  COS  ( X--PI  /4-~D9*PI  /2) 

3230  RETURN 

3290  ’ GAMMA  FUNCTION 

3300  IF  <X<0)  THEM  GOTO  8350 

9310  Y  *  X  +  2 

8320  32  =  SQR  (2*PI/Y)  *  EXP  < Y*LOG  ( Y )  +  ( 1  -- 1  /  ( 30*Y*Y ))/(]. 2-*Y )  - Y ) 
3330  05  -  G2/<X*<X+1>> 

3340  RETURN 

2350  ’NEGATIVE  ARGUMENT 


■  ?  rs'..' 

1  4  O  0 
24  10 
34  20 
'432' 
14  40 
'4rrO 


K57.  ==  0 

IF  (X  >-  0;  THEN  GOTO  8420 

H57.  «  K:57.  -i-  1 

X  *  X  1 

GOTO  3300 

GOSUB  3290 

FOR  7  i  TO  K57. 

X  ..=  y  ...  t 
03  -  03/ X 
NEXT  J 


>0  KEY  OFF 

is  dim  iNTEi'i  ■: -i-k>' 

20  CL 2 
“0  PRINT 

TO  PRINT  "AR-ASO  PLOTTING  ROUTINE: 

50  PRINT 

SO  PRINT  "THIS  PROGRAM  WILL.  GIVE  A  CONTOUR  PLOT  OF  THE  SPOT  OF  ARAGO  AT  A 
70  PRINT  "CHOSEN  DISTANCE  BEHIND  AN  ANNULAR  TARGET.  THE  APERTURE  AND  OBSTACLE 
PR  HIT  "C  AM  l?:I  D 1 3FLACED 
30  PRINT 

?0  PRINT  "THIS  PROGRAM  ASSUMES  A  WAVELENGTH  OF  6320  ANGSTROMS.  INPUTS  ARE 

100  PRINT  "THE  INNER  AND  OUTER  DIAMETERS,  THE  TOTAL  RADIAL  DISTANCE  OFF  THE 

110  PRINT  "OPTICAL  AXIS  FnR  COMPUTATION,  THE  ABERRATION  COEFFICIENTS  FOR  DEFOCUS 

120  PRINT  "AND  SPHERICAL  ABERRATION,  THE  GAUSSIAN  CONSTANT  AND  THE  DISTANCES  Z1 
125  PRINT  "AND  Z2  (MEASURED  FROM  THE  APERTURE  AND  OBSTACLE  RESPECTIVELY) . 

130  PRINT 
.140  PRINT 

'50  INPUT  "INPUT  THE  TOTAL  RADIAL  DISTANCE  TO  BE  CALCULATED  (IN  METERS) "?R1 


360  PRINT:  INPUT  "WHAT  IS  THE  NAME  OF  THE  FILE  TO  STORE  THE  PLOT  UNDER" : F2* 
370  FILES  •■=  "8:  "  +  F2.S 
ISO  ON  ERROR  GOTO  390 
390  RESUME  NEXT 

IOO  PRINT:  INPUT  "DO  YOU  WISH  TO  BYPASS  THE  CONTOUR  PLOT  (Y/N) " ; PBYS 
401  ’ 

102  '  CALCULATE  ABERRATION  AND  GAUSSIAN  FACTORS 


-  r  •  --  • •  it: -- 


'3  17  AB  1  i  '"t;  '  i  ,  '  l.'J-lsU 

413  A 32  -  ( t-B  4,4020 '*•  3  2  ' 

41?  AB3  -  i - 1  c:  -C  ^ 7 !  /  D B  D C  • 

420  AB4  -  •'it. -i  ■=■  (DD*D0-;  E  OP' 

•422  (3 “2  -  E  <P  \  -3 0>D»D  ■-*CD*  ED  '  ■ 

•423  2  =-  203  (  12  :  r»  ■  nr  nr  ••  i  *  .  -  - 

D  >  D  *■  D  f  •  0  /  ■  D  D 0  Z:  D  D  *  0  D 1  ;  1  > 

24  C  0 E F 1  =  A  0  :  ■*  G  -  1  •+  G  F 1  '  ■  A  Dl  »  A  3 1  ' 


rn  >n  .  i.  ~nn..»  nn  ■  , 


D>  0.'  :DD*0C  •  5  -!••  G  .WO-4 


426  C0EF3  =  2H3F1  *GF2  <■  ■:  1  --B>*W040'*  (  .  7  •  D  /  (DD-i-D  D'  )  )  -*C/  (AB1-AB2) 

■429  10  «  00FF1  30EF2  -  CGEF3 

430  PRINT  :  PRINT  "COMPUTING  ...” 

500  D?  =»  0 
510  R  =  R  1/3:2  :• 

520  FOR  I  «  1  TO  440 

525  FOR  J  =  1  TO  2 

526  IF  J  =  1  GOTO  530 

527  2  22 

523  TEMP  =  EPS 
52?  GOTO  535 

530  ? 

531  *  CALCULATE  ARGUMENT  FOR  BESSEL  FUNCTION  AND  EVALUATE  RADIAL  INTENSITY 

FUNCT I ON 

>rr  -*•  '“•>  y 
-J  ^ 

335  X  ■=  (DD  *  TEMP  ■*  I  *  R  *  PI-/  (6. 32BE-07  *  2) 

540  GO CUB  8000 

541  IF  J  =  1  THEN  TEMP 2  -  J7  ELSE  GOTO  550 

542  GOTO  55? 

550  INTENCI)  -  COEFi  *TB1P7-7EMP2  +  C0EF2*J7*J7  -  C  0  E  F  3  T  E  M  P  2  *  J  7 
55?  NEXT  J 
560  NEXT  I 
570  CL 5 

390  * 

39 1  ’  TRANSFER  TO  COLOR  DISPLAY 

e?2  ’ 


TRANSFER  TO 


900 

PR  I NT :  I NPU 

T  "MAKE  S 

LIRE 

i j 

V 

I  DEC) 

SCREEN 

IS  C 

M  AN 

D  HIT  RETURN 

WHEN  READY" 

?  1 0 

DEF  BEG  -  0 

920 

POKE  &H410, 

<  PEEK  ( .EH*- 

•  t  »  * 

H 

mo  k: 

iCF)  OR 

?<h  i  <: 

930 

SCREEN  1,0, 

0,0  :  BCR 

c  rr  h  ( 

!  j 

:  W I 

JTH  80 

940 

LOCATE  ,  ,  1  , 

6,3  ;  nEr.; 

950 

GOTO  ?P5 

960 

96  1. 

’  TRANSFER 

BACK  TO 

A  i'“> 

VOMIT 

"FO 

DEF  GEG  --  0 

900 

POKE  ?■:!■!  4  1 0  , 

(PEEK (OH 

l  i 

OR  7<! 

■!  30 ) 

390 

SCREEN  0:  W 

I DTH  40: 

U  I  n 

TH 

30 : 

LOCATE 

■i  i  1  n 

12.,  1 

3:  DEF  '3 EG:  G 

□TO  1950 

995 

SCREEN  2:  W 

TDTI4  30 

1 

IF  PBY5  =* 

"Y"  GOTO 

1 20 

I) 

1.001 

1* 

1 002 

’  CONTOUR 

plotting 

F:o 

JT 

INE. 

FIND 

INTEL 

S  t  TV 

AT  RADIAL  01 

STANCE  FROM 

SCREEN 

CENTER  ( X 

t*W'» 

-  ■  i 

0)  „ 

]  007 

1 004 

FOR  I  =  1 

TO  30 

1 005 

FOR  J  =  1 

TO  24 

1  00  7 

■n  -■  j*  i .  7 

1 0 1 0 

R  2  -=  C I  NT  ( 

(3QR  (  ( 1 70 . 

1  -*:■ 

3) 

v  (  *  7f 

>-  JI  G! 

(7 

20— T 

*0) * (320-1 >0) 

)  ) 

PAUSE  FOP  PRINTING  AND  STORING  OF  P 


1012  IF  F 
1015  GQGL'3  2070 

1 020  ON  CNN  ‘5 JVC  1  ■  ;  ;>  '•*?•>,  10/' 

1030  AT  ~r\ ;r  *  •  o  ; 

1040  ;-3ra  i  ■  .*•.< 

LOGO  A "5  -  oir-:--  ; 

1 060  GOTO  1  L  3 

10  70  AS  s  C!  IR  ‘  ■  1  77* 

-‘'S'"  7  C  7  ’  '  " ; 

1090  AS  -  c;  -F  A  173; 

1100  GST3  I"'.;. 

1150  ltutz  j,  : 

1160  PRINT  ,7: 

1170  NEXT  2 
■  1  T  NEXT  I 
1.1 7"  30  SUB  2000 
1  1  GO  ’ 

1181  ’  PAUSE  FOR  PRINTING  AND 

•t  i  G  7  ■’ 

L 1 33  FOR  1  -  1  TO  5000:  NEXT  K 

1  IBS  IF  CPLOT$  -  " N "  GOTO  1200 
1190  DEF  3 EG-  7HPS00 
1 1 93  B3AVE  F I L.  1  $ »  0 , 57-14000 
1 200  ? 

1201  ’  DRAW  AXES 

1 202  •' 

1203  CL3 

1210  FOR  I  «  1  TO  20 
!  220  LOCATE  1,12 
1230  PRINT  CHRf.  (179); 
i  2  TO  NEXT  I 
1750  LOCATE  21,12 
1260  PRINT  CHRiS  (192? 

12"'..  FOR  I  ■  1  TO  56 
z  1  :*  LOCATE  21 , 12-*- 1 
1  PRINT  CHR$  <  196)  1 
1  .’:)■)  • 

1301  LABEL  Y  AXIS 

1302  ’ 

1 305  NEXT  I 

'7-0  LOCATE  16.,  9 

1 320  PR I NT  1 1 CHR* ( 1 97 ) 

1370  LOCATE  11,9 


10/0, 1090, 1 1 


DRAW  AXES 


j  2  40 

•  050 
t  *•/ . » 


0  4  ') 

PRINT  2 'I  CHR 

3  SO 

7  ,jO 

L..OCATE  6  ,,9 

370 

L..GCA  1 12  1 , 9 

330 

PRINT  44 CHR 

17?  0  FOR  I  ;=  1  TO  9 
1  400  LAB IT  «  "NORMAL I  21 
1  4 1 0  FOR  I  ~  1  TO  9 
1420  LOCATE  /  +■ »;  ,  " 

1  130  PRINT  M  IDS  1 1  A  R  IS. 


LAB2T  ~  "INTENSIT 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


:  u.->  i- 


-  >  L  J  -,.-1  .  L  a  ,  -■ 

id'j  fr in-  " r " ■ 


.<  AXIS 


!  '.j4;.>  Lu.-.i  I  .  -  ,4,  -fc- o 

•  550  PRINT  ”  iADIAL  DISTANCE  (.001  M) "5 

1  360 

>361  ’  PLOT  RADIAL  INTENSITY  FUNCTION 

•  cr  .*  “t  ? 

’565  lJD  =  1X30  *  F:1 
570  F  OW  “  I  NT  <  LOG  ( DDD )  /  LOG  (10)) 

I  3 so  SCALE  -  DDD/lO-'-POW 

. 590  IF  GCA  E  <  4  THEN  FACT  *  2  ELSE  FACT  »  1 
1 600  UNIT  =  30/ (SCALE  *  FACT) 

. 1  •  .<  FOR  I  ---  1  TO  I  NT  (SCALE  *  FACT) 

!  A7  0  LOCATE  21,12  +  t*UNIT 
i  *,30  PRINT  CHR*  ( 197) 

3  40  LOCATE  22,10  +  I*UNIT 
■  ■:  30  PRINT  1*1 O"  ROW / FACT 
1  -.60  NEXT  I 

•  -:V3  YO  =  164  -  CINT(  10*40) 

.  06  70  =  93 

•  FOR  I  -  1  TO  440 

6  CO  V  =  164  -  Cl  NT  (  ( INTEN  ( I )  *  160/4) 

:  :  33  X  »  93  +  I 

.',95  LINE  ( X  ,  Y ) - ( XO,  YO ) 

I  .’00  XO  a  X 
' 705  YO  «  Y 
1710  NEXT  I 
: ’70  LOCATE  25,6 

i.  730  PRINT  USING  "O.D.=  .  tttt#  M.  I.D.  =  .  ####  H.  Z  DIST.  =  #.###  M.  GAUSS 

:mst.  •-=  #.#)!" ;dd.,d.,  zi,gc; 

i  733  LOCATE  1 , 1 

1740  FOR  I  =  1  TO  5000 

!  7  45  NEXT  I 

:  750  IF  I PLOTS  a  "N”  GOTO  1300 

1760  ’ 

1761  ■’  STORE  INTENSITY  PLOT  ON  DISK 

1762  ’ 

1765  DEF  3EG*  &HB800 
!  770  D3AVE  FIL2*, 0, &H4000 

<300  PRINT:  INPUT  "DO  YOU  WISH  A  LISTING  OF  THE  INTENSITIES ( Y/N) "  J  Dl* 

1 3 10  IF  Dl$  *  "Y"  THEN  GOTO  1320  ELSE  GOTO  1940 

1 820  LF'R  I  NT :  LF'R  I  NT :  LPR I  NT :  LPR I  NT :  LPR  I  NT :  LF'R  I  NT :  LPR I  NT :  LPR I  NT  :  LPR I  NT 
•330  LPR I NT  :  LPR I NT  :  LPR I NT  :  LPR I NT  :  LPR I NT  :  LPR I NT  :  LPR I NT 
1340  LF'R  I  NT  "RADIAL  DISTANCE  (MM)  INTENSITY  EXTREMUM  TYPE 

1350  LPR  I  NT  " -  -  - 

1155  I'"  10  -  INTEN(l)  >  0  THEN  EXT*  »  "MAXIMUM"  ELSE  EXT*  =  "MINIMUM" 

»  360  LT’R  I  NT  US  I NG  "  # .  ###  # .  tttt#  "  ;  O,  10; 

1363  I. -•PINT  "  M;EXTS 

1 269  S I GN  a  GGN ( 1 0- 1 NTEN ( 1 ) ) 

<370  FDR  I  =  1  TO  439 

•  300  IF  SGN  (  I  NTEN  ( I  + 1 )  -  INTEN  ( I )  )  =«  SIGN  GOTO  1930 

< 090  IF  SIGN  >  0  THEN  EXT$  a  "MAXIMUM"  ELSE  EXT*  =  "MINIMUM" 

1900  SIGN  =  SGN <  INTEN ( 1  + 1 )  —  INTEN  ( I )  ) 

1910  LPR I NT  USING  "  #.###  #.####" 5 I*R*1000, INTEN ( I ) J 

1720  LPR I NT  "  " ;E- 


REPRODUCED  AT  GOVERfMfeNT  EXPENSE 


1  ? '3  1  CLS 

I960  I i  IF'UT  "DO  YOU  WISH  TO  MAKE  ANOTHER  RUN  (Y/N)  "  5  D$ 

1970  IF  D*  =  “N"  THEN  END  ELSE  GOTO  130 

2"0'>  ’  | 
3  ■.'•’01  •  LAE  EL  CONTOUR  PLOT 

O'!  02  ’ 

O'.N'E  LOCATE  24,1 

2010  PRINT  "LEGEND:  " :  CHR*  <221  >  5  11  I  <  .01  la  "  5  CHR*  <  170)  ;  "  I  <  .05  la  "5CHR4U1 
-J-'  ;  "  I  '  .  1  la  ■’  ?CHR$<176)  ;  ”  I  <  .5  la  " :  CHF$  <0)  ;  '*  I  >  .5  la"; 


j  SHADING  FOR  CONTOUR  PLOT  POINTS 


IF  INI  EM (RE)  >  .002  GOTO  2060 

DUN  ■---  5 

RETURN 

IF  IN TEN (R2>  >  .05  GOTO  2090 

D'jM  ••=  4 

RETURN 

IF  INTEN (R2)  >  .1  GOTO  2120 

BUM  ---  3 

RETURN 

IF  INTEN <R2)  >  .5  GOTO  2150 

BUM  ~  2 

RETURN 

DUM  =  1 

RETURN 


COai  '  SUBROUTINE  TO  FIND  BESSEL  FUNCTIONS  OF  THE  FIRST  KIND 


Cl  00 
C 1 1 0 
0120 
1 30 
A 1  40 


c  1  70 
1 00 
II O.;, 


I -  9.  999999E-06 
rC  -  3.  14.159 
v2  =  X  *  X 

if  (  < X=0)  AND  <D9~1 ) )  THEN  8240 
IF  -X  '  .15)  THEN  8260 

IF  CD?  <>  0)  THEN  8090 

GOTO  8140 
X  3  •=  X 
X  =  D9  +  1 
COSUB  8290 
X  =  y  3 

86  =  ( X/2)  '"D9/G5 
T  _■  --  S6 
17.  =  0 

IF  ABS <T3)  ABSCS6*T1)  THEN  8220 
n:  “  17.  +  1 
T4  =  T3 

T3  =  -T4  *  X2  *  .25/(17.  *  <D9  +  17.)) 
B6  =  S6  +  T3 
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